# Stellar activity with LAMOST – II. Chromospheric activity in open clusters

Stellar activity with LAMOST – II. Chromospheric activity in open clusters Abstract We use the LAMOST spectra of member stars in Pleiades, M34, Praesepe, and Hyades to study how chromospheric activity varies as a function of mass and rotation at different age. We measured excess equivalent widths of H α, H β, and Ca ii K based on estimated chromospheric contributions from old and inactive field dwarfs, and excess luminosities are obtained by normalizing bolometric luminosity, for more than 700 late-type stars in these open clusters. Results indicate two activity sequences in cool spot coverage and H α excess emission among GK dwarfs in Pleiades and M dwarfs in Praesepe and Hyades, paralleling with well-known rotation sequences. A weak dependence of chromospheric emission on rotation exists among ultrafast rotators in saturated regime with Rossby number Ro ≲ 0.1. In the unsaturated regime, chromospheric and coronal emission show similar dependence on Ro, but with a shift towards larger Ro, indicating chromospheric emission gets easily saturated than coronal emission, and/or convective turnover time-scales based on X-ray data do not work well with chromospheric emission. More interestingly, our analysis shows fully convective slow rotators obey the rotation–chromospheric activity relation similar to hotter stars, confirming the previous finding. We found correlations among H α, H β, and Ca ii K emissions, in which H α losses are more important than Ca ii K for cooler and more active stars. In addition, a weak correlation is seen between chromospheric emission and photospheric activity that shows dependence on stellar spectral type and activity level, which provides some clues on how spot configuration varies as a function of mass and activity level. stars: activity, stars: chromospheres, stars: late-type, stars: rotation 1 INTRODUCTION Stellar magnetic activity, triggered by magnetic fields and maintained by a magnetic dynamo, is an ubiquitous phenomenon in late-type stars. Such activity includes a variety of phenomena, e.g. appearance of cool star-spots that are thought to be the fingerprints of emergent magnetic lines on the photosphere; the chromospheric and coronal heating that are responsible for excess chromospheric emissions and thermal X-rays, respectively; impulsive events like flares; and other solar phenomena. Stellar mass and rotation are thought to be two main basic parameters that affect the dynamo configuration and thus the properties of stellar magnetic activity. On the other hand, the loss of stellar angular momentum through the magnetized stellar wind depends on the magnetic field configuration on stellar surface (Gallet & Bouvier 2013). As the rotation of the star is slowed down due to angular momentum loss, its magnetic field weakens (or the efficiency of the dynamo decreases), and thus activity decays. The correlations between rotation, activity, and age are evident from observations of various tracers of activity, e.g. the rotation–activity connection is found based on stellar X-ray emissions (e.g. Pizzolato et al. 2003; Mamajek & Hillenbrand 2008; Wright et al. 2011; Reiners, Schüssler & Passegger 2014), chromospheric Ca ii HK emissions (e.g. Noyes et al. 1984; Mamajek & Hillenbrand 2008) and H α emissions (e.g. Stauffer et al. 1997; Douglas, Agüeros & Covey 2014; Newton et al. 2017), and photospheric activity indicator such as cool spot fractional coverage (Fang et al. 2016, hereafter Paper I). The connection between rotation and activity among cool stars directly allows us to understand the nature of the magnetic dynamo in these stars, however, such connection appears to be sensitive to tracers of activity, e.g. chromospheric emission shows different dependence on rotation compared to coronal X-ray emission (Douglas et al. 2014). Chromospheric activity (hereafter CA) compasses a wide range of phenomena that produce excess emission (e.g. Ca ii and Balmer lines) with respect to a radiative equilibrium atmosphere, stemming largely from changes in their magnetic field (see a comprehensive review on this topic by Hall 2008). Thanks to classic work of Mount Wilson Observatory Ca ii HK long-term observations that are operated from 1966 through 2003, for the initial purpose of searching stellar analogues of the solar cycle (Wilson 1968), wherein variety of CA properties including chromospheric variations are explored (Wilson 1978; Baliunas et al. 1995, 1998), e.g. the well-known CA long-term variation patterns (cyclic variation, irregular variation, little variation). What followed was a series of attempts trying to gain the general understanding (statistical characters) of CA among nearby late-type stars, e.g. statistical studies based on chromospheric emissions at the core of Ca ii H and K lines among nearby FGK field stars (Henry et al. 1996; Gray 2003, 2006), confirmed the bi-modal distribution in stellar activity (Vaughan & Preston 1980); West et al. (2004, 2008) investigate the fraction of active stars as a function of spectral type and the vertical distance from the Galactic plane, based on a large sample of nearby M-type stars using H α emission line as an activity indicator. Compared to field stars, open cluster members have an advantage for studying CA properties in light of rotation–activity–age connection due to their homogeneity in age and chemical composition. Therefore, CA in several open clusters got special attention, e.g. the young open cluster Pleiades (Soderblom et al. 1993) and NGC 2516 (Jackson & Jeffries 2010), the median-age open cluster Hyades (Stauffer et al. 1991, 1997) and Praesepe (Douglas et al. 2014). Based on these investigations a general picture is drawn that CA also varies as a function of both stellar mass and rotation. However, due to the limitation of sample stars used in previous studies in terms of mass coverage and available rotation periods (e.g. Pleiades), the behaviours of CA with respect to stellar mass changes and rotation is still not well understood, which demands further investigation on this issue. LAMOST (Cui et al. 2012) observed millions of stellar spectra for cool stars in our Milky Way (Liu, Zhao & Hou 2015), and more and more stars’ rotation periods become available thanks to recent and on-going photometric surveys, e.g. for ∼1000 Pleiades members (Hartman et al. 2010; Rebull et al. 2016), both of which allow us to initiate systematic study of CA among low-mass stars in different population. In Paper I, we found that almost all Pleiades members show excess emission in various indicators (e.g. Ca ii HK, H α, H β, Ca i λ4226, Mg i b triplet, Na i D lines, and Ca ii IRT lines, see figs 7 and 8 in Paper I) relative to old and inactive field stars, indicating they are chromospherically active. In this paper, we continue to investigate their CA properties in detail among members of nearby open clusters including Pleiades, M34, Praesepe, and Hyades by focusing on three chromospheric indicators: H α, H β, and Ca ii K lines. We examined their incidence and strength of CA as a function of effective temperature, and checked the correlation between activity and rotation. We also investigated the relationship amongst these three chromospheric indicators, and the relations between chromospheric emission and several manifestations of magnetic activity in different stellar atmosphere layers (cool spot coverage, light variation amplitude, and coronal emission), which holds key clues not only to the configuration of the magnetically active region on the stellar surface but also to the magnetically heating processes in stellar outer atmosphere. 2 DATA AND SAMPLE 2.1 LAMOST DR3 The LAMOST, characterized by both wide field of view (5 deg in diameter) and large effective aperture of ∼4 m, is a reflecting Schmidt telescope located at the Xinglong Observatory, China. LAMOST spectra are obtained in low resolution R ≈ 1800 (e.g. ∼3.6 Å around 6500 Å) (see Zhao et al. 2006, 2012; Luo et al. 2015, for more details). The LAMOST Date Release 3 (DR3) was made available to the Chinese astronomical community and international partners during 2015 December, containing nearly 5.3 million stellar spectra that were collected from pilot survey through 3 yr regular survey (autumn 2011 to summer 2015). LAMOST stellar parameter pipeline (LASP) provides the stellar parameters like effective temperature (Teff), surface gravity (log g), metallicity ([Fe/H]), and radial velocity (RV) for ∼3.2 million stellar spectra of AFGK stars whose spectra meet the signal-to-noise ratio criterion (e.g. g-band SNR of SNRg ≥ 15 and SNRg ≥ 6 for the bright and dark nights, respectively) (Wu et al. 2011; Luo et al. 2015). In addition, DR3 also includes a catalogue of ∼0.3 million M-type stars. 2.2 Sample stars We have selected four nearby open clusters with a range of ages from 100 to 700 Myr (Pleiades, M34, Praesepe, and Hyades), and retrieved the available spectra of good quality from LAMOST DR3 archive. The Pleiades (seven sisters) is a young open cluster with an age of ∼125 Myr (Stauffer, Schultz & Kirkpatrick 1998) and metallicity of [Fe/H] ∼ 0.03 (Soderblom et al. 2009), at a distance of around 130 pc (e.g. ∼136 pc; Melis et al. 2014). In addition to LAMOST DR2 archive used in Paper I, DR3 provides spectra for more than 70 targets from the Pleiades potential member catalogues provided by Stauffer et al. (2007) and Bouy et al. (2015) (whose probability is more than 75 per cent). In total, we selected 378 probable Pleiades members with 458 LAMOST spectra having signal-to-noise ratio in r-band SNR r > 10. M34 (NGC 1039) has an age of about 220 Myr (see Meibom et al. 2011, and references therein), located at a distance of ∼470 pc towards the Galactic anticentre (Jones & Prosser 1996). M34 has metallicity close to solar, [Fe/H] ∼ 0.07 (Schuler et al. 2003). Jones & Prosser (1996) provided proper motion membership probability analysis for more than 600 stars with brightness down to V ∼ 16.2 in the M34 field, wherein, the cooler members tend to have lower membership probability due to the faintness, e.g. the typical probability of faint stars with V > 15 (∼K2 and cooler) is ∼40 per cent, while the brighter members have probabilities close to 100 per cent. We selected candidate M34 members with membership probability greater than 40 per cent. However, real cool members might be excluded in this criterion. To collect more member stars, in particular at cooler end, we chose periodic variables detected in M34 (Irwin et al. 2006; James et al. 2010; Meibom et al. 2011), considering them as younger and more active. Unfortunately, only a small number of cool stars in the vicinity of M34 field have been observed by LAMOST (DR3), partly due to their faintness. In total, we collected 45 candidates that have LAMOST spectra with SNR r > 10, most of which are FG-type stars. The Hyades is the nearest open cluster at a distance of about 46 pc (Perryman et al. 1998), and has solar metallicity [Fe/H] ∼ 0.1 (Taylor & Joner 2005; Carrera & Pancino 2011). The Praesepe (NGC 2632; the Beehive Cluster) has a distance of 180–190 pc (e.g. An et al. 2007; van Leeuwen 2009) and a metallicity of [Fe/H] ∼ 0.16 (Carrera & Pancino 2011; Yang, Chen & Zhao 2016). The Hyades and Praesepe are both part of the Hyades supercluster, and seem to share an intermediate-age around 600–800 Myr (Perryman et al. 1998; Fossati et al. 2008; Brandt & Huang 2015). We adopt the membership catalogue provided by Douglas et al. (2014) for these two clusters, and retrieved 197 spectra for 142 Hyades candidates, and 227 spectra for 176 candidates of Praesepe from LAMOST DR3 archive with SNR r > 10. Fig. 1 and Table 1 show the summary of total sample (741 stars), and it is easy to note that more than half (∼60 per cent) of them are late K- and M-type stars. Figure 1. View largeDownload slide Distribution of the sample stars as a function of temperature. The dashed black line shows the distribution for total sample stars, while the grey line represents these stars for which rotation period is now known. The solid lines with different colours correspond to cumulative fraction profiles (right y-axis). Figure 1. View largeDownload slide Distribution of the sample stars as a function of temperature. The dashed black line shows the distribution for total sample stars, while the grey line represents these stars for which rotation period is now known. The solid lines with different colours correspond to cumulative fraction profiles (right y-axis). Table 1. Basic information about the sample stars in current work. Cluster  d  Age  [Fe/H]  E(B − V)  Samplea    (pc)  (Myr)    (mag)    Pleiades  136  ∼125  +0.03  0.03b  378 (458), 231  M34  470  ∼220  +0.07  0.07  45 (49), 21  Praesepe  184  ∼650  +0.16  0.027  176 (227), 76  Hyades  46  ∼650  +0.11  0.001  142 (197), 24  Cluster  d  Age  [Fe/H]  E(B − V)  Samplea    (pc)  (Myr)    (mag)    Pleiades  136  ∼125  +0.03  0.03b  378 (458), 231  M34  470  ∼220  +0.07  0.07  45 (49), 21  Praesepe  184  ∼650  +0.16  0.027  176 (227), 76  Hyades  46  ∼650  +0.11  0.001  142 (197), 24  aTotal number of stars (spectra), and those with rotation period. bThere is significant differential extinction in the Pleiades, see Soderblom et al. (1993). View Large 3 QUANTIFYING CHROMOSPHERIC ACTIVITY It is well known that enhanced emission of H α, produced by collisional excitation in the relatively dense chromospheres of low-mass dwarf stars, is the strongest and widely used activity indicator (e.g. Stauffer & Hartmann 1986; Soderblom et al. 1993; West et al. 2004, 2015; Newton et al. 2017). In this work, H α line was used as main CA indicator, considering that more than 50 per cent of our sample stars are late K- and M-type stars (see Fig. 1). However, we also investigated the behaviours of other two activity indicators, Ca ii K and H β lines. 3.1 Excess equivalent width The equivalent widths were measured from RV-corrected spectra of all sample stars. For FGK dwarfs, we adopted RV values derived by LASP, while we measured RVs for M dwarfs following methodology described in Paper I. The equivalent widths of H α line (hereafter EWH α) were measured using the formula (equation 4) in Paper I, where the line was centred at 6563 Å over 12 Å bandpasses and continuum flux was taken to be the average flux between 6547–6557 and 6570–6580 Å. The measurements of EWH α for entire sample are plotted against effective (quiescent) temperatures (see Appendix D for details) in Fig. 2, wherein the mean values of EWH α for inactive reference stars (a large sample of very less active field FGKM dwarfs with solar metallicities over the temperature range 3000–6500 K, see Paper I for more details) are treated as the minimum or basal value of chromospheric emissions, and shown as black solid lines. It is clear from the figure that the members of Pleiades show higher EWH α compared to inactive counterparts. In fact, the correlation between EWH α and CA level is complex, due to presence of non-monotonic behaviour (Cram & Mullan 1979), e.g. first the H α absorption increases with increase in activity, then fill-in, and eventually appear as emission line. Therefore, among the less active stars at the same temperature, the stronger Hα absorption line indicates its higher activity. In this work, we adapted simple approach by ignoring such non-monotonic effect, and consider excess EWH α above the basal value as the tracer of pure chromospheric emission. To determine the excess H α emission, we calculated the excess equivalent widths (the ΔEWH α in Paper I, hereafter denoted as EW$$^{\prime }_{\text{H}\,\alpha}$$) by subtracting the basal value at the same temperature, EWH α, basal, i.e.   $$\text{EW}^{\prime }_{\text{H}\,\alpha } = \text{EW}_{\text{H}\,\alpha } - \text{EW}_{\text{H}\,\alpha ,\text{basal}}.$$ (1)Similarly, we obtained the excess equivalent widths, EW$$^{\prime }_{\text{Ca K}}$$ and EW$$^{\prime }_{\text{H}\,\beta}$$, for Ca ii K and H β lines, respectively (see Appendix A for details). Figure 2. View largeDownload slide Top: The equivalent widths of H α line for stars in Pleiades (left), and M 34, Praesepe, and Hyades (right). Black solid lines show the mean relations for inactive reference stars, upper and lower 3σ to the mean value are shown as grey lines. Filled symbols represent probable single members and open symbols denote probable binary members or non-members. For stars having LAMOST multi-epoch observations, we displayed their average values rather than individual observational values. Middle: The excess equivalent width of H α, EW$$^{\prime }_{\text{H}\,\alpha}$$. Bottom: The excess fractional luminosity of H α, $$R^{\prime }_{\text{H}\,\alpha}$$, where the grey solid (dashed) lines denote the stars locating in the upper 3σ (σ) regime. Figure 2. View largeDownload slide Top: The equivalent widths of H α line for stars in Pleiades (left), and M 34, Praesepe, and Hyades (right). Black solid lines show the mean relations for inactive reference stars, upper and lower 3σ to the mean value are shown as grey lines. Filled symbols represent probable single members and open symbols denote probable binary members or non-members. For stars having LAMOST multi-epoch observations, we displayed their average values rather than individual observational values. Middle: The excess equivalent width of H α, EW$$^{\prime }_{\text{H}\,\alpha}$$. Bottom: The excess fractional luminosity of H α, $$R^{\prime }_{\text{H}\,\alpha}$$, where the grey solid (dashed) lines denote the stars locating in the upper 3σ (σ) regime. 3.2 Excess fractional luminosity As shown in Fig. 2, the strong change of EWH α (and EW$$^{\prime }_{\text{H}\,\alpha}$$) with temperature does reflect not only the intrinsic change of CA, but also the effects of the drop in the nearby continuum level. To remove the dependence of the measured EWH α on the surrounding continuum level, and to compare the strength of CA in different spectral types, we used the fractional luminosity of H α (normalized with the bolometric luminosity) (e.g. Reid, Hawley & Mateo 1995; West et al. 2004). However, quantifying fractional luminosity for a given line requires a flux-calibrated spectrum and the distance. For all activity indicators in this work, the excess fractional luminosities were derived from the excess equivalent widths using a distance-independent value, χ, the ratio between the continuum flux near the line of interest and the bolometric flux, following commonly used χ-methods (e.g. West et al. 2004). Using the χ ratios of H α line (χH α) obtained from the model spectra (see Appendix C), we obtained the H α excess fractional luminosity (denote as $$R^{\prime }_{\text{H}\,\alpha}$$), the ratio of excess H α emission luminosity ($$L^{\prime }_{\text{H}\,\alpha}$$) to the bolometric luminosity (Lbol), as follows:   $$R^{\prime }_{\text{H}\,\alpha } \equiv \frac{L^{\prime }_{\text{H}\,\alpha }}{L_{\text{bol}}}= {\chi }_{\text{H}\,\alpha } \times \text{EW}^{\prime }_{\text{H}\,\alpha }.$$ (2)In similar way, we estimated Ca ii K and H β excess fractional luminosities, $$R^{\prime }_{\text{Ca K}}$$ and $$R^{\prime }_{\text{H}\,\beta}$$, respectively (see Appendix A). These values provide the fraction of respective bolometric luminosity emitted from the chromosphere in the lines, and are analogous to the chromospheric Ca ii emission ratio $$R^{\prime }_{\text{HK}}$$ (e.g. Noyes et al. 1984) and $$R^{\prime }_{\text{IRT}}$$ (e.g. Marsden, Carter & Donati 2009; Jackson & Jeffries 2010). 3.3 Comparison with literature measurements We compared our EWH α measurements with Douglas et al. (2014), wherein they measured EWH α from medium resolution (about two times higher than LAMOST) spectra, for common stars in Hyades and Praesepe (see Fig. 3), and found to be in good agreement. Furthermore, Soderblom et al. (1993) and Soderblom, Jones & Fischer (2001) obtained excess H α emission from high-resolution echelle spectra for about 100 FGK dwarfs in the Pleiades and about 50 FGK dwarfs in the M34, respectively, using a synthetic subtraction technique, i.e. the H α profile of each target star is subtracted by the rotationally broadened H α of inactive counterpart. In Fig. 4 we compared their measurements of excess H α emission with our measurements for common stars in the Pleiades and M34, and found they are in close agreement. However, it seems there still exist systematic offset in fractional luminosity for less active hotter stars with $$\log R^{\prime }_{\text{H}\,\alpha } \lesssim -4.5$$, e.g. our measured $$R^{\prime }_{\text{H}\,\alpha}$$ are ∼0.2 dex larger. In fact, as discussed in Appendix B, a small error in excess equivalent width would result in larger error in fractional luminosity (e.g. see Table B1), particularly when the excess emission is very weak. In a careful check on the excess equivalent width, we indeed found there is a small offset (∼0.1 Å, see the upper panel of Fig. 4) that is comparable with the uncertainties due to 100 K uncertainty in temperature (see Table B1). Figure 3. View largeDownload slide Comparison of our measured EWH α with those of Douglas et al. (2014) for common stars in Hyades and Praesepe. Figure 3. View largeDownload slide Comparison of our measured EWH α with those of Douglas et al. (2014) for common stars in Hyades and Praesepe. Figure 4. View largeDownload slide Comparison of H α excess emissions with those of Soderblom et al. (1993) (Pleiades members) and Soderblom et al. (2001) (M34 members) for common stars. Figure 4. View largeDownload slide Comparison of H α excess emissions with those of Soderblom et al. (1993) (Pleiades members) and Soderblom et al. (2001) (M34 members) for common stars. 3.4 Variability in chromospheric emission The larger scatter in Figs 3 and 4 are of course partly from the uncertainties of our measurements. However, we believe that they are mainly related to intrinsic variation in activity, which is well-known for Pleiades members (e.g. Soderblom et al. 1993). To further check the variation in H α emission, we have plotted EWH α against the temperature in Fig. 5, in which we show minimum and maximum values of EWH α for stars having multi-epoch spectra. From the figure it is clear that very cool members show high variation in H α emission, which is evidently larger than H α equivalent width measurement errors, indicating intrinsic variability of CA. Also, it shows that more active stars have large equivalent width variations, indicating stars with high level activity suffer more often large variations, which may be due to short-time-scale evolution of active regions or flare-like events. Figure 5. View largeDownload slide Variation of EWH α for stars having multi-epoch observations. Also shown by grey-colour solid line are the basal values. Figure 5. View largeDownload slide Variation of EWH α for stars having multi-epoch observations. Also shown by grey-colour solid line are the basal values. 4 RESULTS AND DISCUSSION 4.1 Chromospheric activity as a function of mass As shown in top-left panel of Fig. 2, the EWH α is generally a strong function of stellar effective temperature, i.e. EWH α increases with decreasing temperature, from H α absorption at Teff ∼ 6500 K (EWH α ∼ −2 Å) to H α emission at Teff ∼ 3000 K (EWH α ∼ 10 Å). Such trend in the Pleiades was also observed by Stauffer & Hartmann (1987). The upper envelop of EWH α versus Teff profile indicates the activity levels of most active stars at a given temperature, and provides the clues on chromospherically cooling effects for stars with different mass. Also, it shows the dependence of the upper envelop on spectral types among Pleiades members, e.g. the EWH α increases moderately as temperature decreases for late F- and G-type stars (Teff ≳ 5000 K), and for K-type stars it goes into a plateau where exists a slight increment of EWHα, and then it increases sharply when Teff ≲ 4000 K. However, at the same time, the spread in EW H α also increases. Such feature is more overt in the middle-left panel of Fig. 2. The lower EWH α of Pleiades members follows the mean relations of inactive reference stars, which varies with temperature, i.e. from EWHα ∼ −2.5 Å at Teff ∼ 6500 K to EWH α ∼ 0 Å at Teff ∼ 3800 K. We noticed that some Pleiades G-type members have larger EWH α values compared to inactive counterparts, indicating they have shallower H α absorption line profile (i.e. the core of the line is filled in). The excess H α emission for all clusters are shown in middle panels of Fig. 2. It is clear from the figure that most of Pleiades candidates show excess emission compared to inactive counterparts, i.e. EW$$^{\prime }_{\text{H}\,\alpha }>0$$, indicating they are in chromospherically active phase, which is consistent with its young age. However, a small fraction of Pleiades stars have weaker excess H α emission compared to other stars in the same spectral type, e.g. about 40 very low-mass stars with Teff < 4000 K have EW$$^{\prime }_{\text{H}\,\alpha }<0.5$$ Å, while other members in this temperature regime have evident larger EWs (e.g. EW$$^{\prime }_{\text{H}\,\alpha }>2$$ Å), making a gap around Teff ∼ 3500 K in the diagram. One general interpretation is that these stars may be field dwarfs rather than Pleiades members. We found that about one -third of these stars are photometrically single candidate members according to their location on CMDs. To further check their membership, we plotted their RVs distribution in Fig. 6, and found about half of them have RVs from −8 to 20 km s−1 suggesting they are likely Pleiades members within 3σ-RV regime by adopting the mean RV value of μ ∼6 km s−1 and a scatter of σ ∼5 km s−1 for Pleiades (see Fig. 6). Several stars (e.g. about one-sixth of them) are probably Pleiades members which pass both the CMD and RV criterion, such as BPL167, DH 668, and DH 908. On the other hand, we noticed that all these less active stars have no rotation periods available till the date. In other words, stars having detected rotation periods show larger activity in the cooler end, which makes the issue complex. Therefore, their memberships and potential light variation through monitoring indeed need to be further investigated. Figure 6. View largeDownload slide RV distributions of Pleiades candidate members with Teff < 4000 K in inactive regime, shown by the grey line. The red line is a Gaussian fitting, which gives μ ≈ 6.2 km s−1 and σ ≈ 4.6 km s−1. The expected Pleiades RV at 5.9 km s−1 provided by Mermilliod, Mayor & Udry (2009) is marked as a black dashed line. Figure 6. View largeDownload slide RV distributions of Pleiades candidate members with Teff < 4000 K in inactive regime, shown by the grey line. The red line is a Gaussian fitting, which gives μ ≈ 6.2 km s−1 and σ ≈ 4.6 km s−1. The expected Pleiades RV at 5.9 km s−1 provided by Mermilliod, Mayor & Udry (2009) is marked as a black dashed line. The large contrast of continuum flux at H α line between M-type and hotter stars dilutes their larger differences in equivalent widths, which makes the upper envelops of $$R^{\prime }_{\text{H}\,\alpha}$$ among these stars have similar values, e.g. $$\log R^{\prime }_{\text{H}\,\alpha }\sim -3.7$$ (see bottom panels of Fig. 2). Still, the upper envelop of $$R^{\prime }_{\text{H}\,\alpha}$$ in Pleiades M-type stars (e.g. Teff ≲ 3900 K) show a slight increment with decreasing effective temperature, confirms the previous finding that the low-mass stars in young clusters with age less than 150 Myr show an increase in activity strength (LH α/Lbol) towards lower mass (Hawley et al. 1999). Also, previous studies show that later spectral type stars have lower activity level among late M-type stars (e.g. Hodgkin et al. 1995; West et al. 2004). Though such trends are evident from $$R^{\prime }_{\text{H}\,\alpha}$$ when Teff ≲3300 K, and $$R^{\prime }_{\text{H}\,\beta}$$ (see Fig. A1), we cannot make a solid conclusion about this issue based on current data. Another remarkable feature in Fig. 2 is a large scatter of $$R^{\prime }_{\text{H}\alpha}$$ among hotter Pleiades members. In particular, there evidently are two sequences among GK-type Pleiades candidates, i.e. one sequence of stars have similar activity strengths, and other sequence being with lower activity strength showing an increment towards lower temperature. The largest difference of activity level between these two sequences, e.g. 0.5–0.6 dex in $$\log R^{\prime }_{\text{H}\alpha}$$ (a factor of about four in $$R^{\prime }_{\text{H}\alpha}$$), occurs at about 5000 K. Such a difference in activity levels is unlikely due to measurement uncertainties, indicating bi-sequence behaviour among these stars. Therefore, apart from the stellar mass, other factors may play a role in producing the observed dispersion pattern. In fact, we have noticed these two sequences correspond to rotational sequences (see Section 4.2.1). It is known that the CA depends not only on stellar mass but also on age, e.g. the H α emission lasts longer in cooler stars (Stauffer et al. 1991; Hawley et al. 1999). Indeed, Fig. 2 shows a clear age effect on CA, i.e. compared to Pleiades, Hyades and Praesepe are typically less active, and M34 with intermediate age between Pleiades and Praesepe/Hyades shows activity levels in-between them. Additionally, the onset of H α direct emission occurs at about Teff ∼5000–5500, 4500–5000, 4000–4500 K for Pleiades, M34, and Hyades/Praesepe members, respectively. Furthermore, two activity sequences appear among cooler members (e.g. late K- and M-type stars) in Praesepe and Hyades, similar to the feature observed for GK-type members in Pleiades. Interestingly, the members of Praesepe and Hyades with Teff less than 3600 K (spectral type later than M1–M2) have activity levels similar to those of Pleiades active members; it indicates that the activity lifetime for M1–M2-type stars is around 600–700 Myr, a value agreed with that reported by West et al. (2008). It is possible to get such kind of activity-time-scale for earlier spectral type stars by comparing M34 with Pleiades, however, we could not able to provide any solid conclusion due to limited sample of M34 candidates in our data. 4.2 Rotation–activity connection 4.2.1 Rotation and activity sequences As discussed above, there exist two sequences among GK-type stars in Pleiades, and M-type stars in Praesepe and Hyades in $$R^{\prime }_{\text{H}\alpha}$$–Teff diagram, i.e. H α emissions in Pleiades begin to bifurcate at ∼6000 K and merge again at ∼4000 K; similar sequences in Praesepe and Hyades appear at ∼4000 K and disappear at about 3100 K. In order to understand further on this scenario, we provide $$R^{\prime }_{\text{H}\alpha}$$ as a function of temperature again in Fig. 7, wherein shown only members with known rotation periods. For comparison purpose, we also plot the corresponding rotation period (Prot) distribution in the upper panels, where the Prot values were collected from the literature (see Appendix E). The colour contrast in Fig. 7 denotes the Rossby Number Ro (Ro = Prot/τ, normalized rotation period by the convective turnover time τ). The τ was estimated from stellar mass (mass were converted from the quiescent photosphere temperature using PARSEC models, see Appendix D) using the correlation between convective turnover time and mass, derived by Wright et al. (2011). It is clear that stars in the upper sequence rotate faster than those in the lower activity sequence. Figure 7. View largeDownload slide Rotation period (top), $$R^{\prime }_{\text{H},\alpha}$$ (middle), and spot filling factor (bottom) as a function of Teff for stars with rotation period available. Left: Pleiades. Right: M34, Praesepe, and Hyades. The colour gradient indicates Rossby numbers. Overplotted in top-right panel are contours of constant Rossby numbers (Ro = 0.01, 0.1, and 1.0; grey dashed lines) based on the empirically calibrated mass–convective overturn time relation by Wright et al. (2011). The canonical rotational I/C sequences are illustrated by black lines in top panels. Figure 7. View largeDownload slide Rotation period (top), $$R^{\prime }_{\text{H},\alpha}$$ (middle), and spot filling factor (bottom) as a function of Teff for stars with rotation period available. Left: Pleiades. Right: M34, Praesepe, and Hyades. The colour gradient indicates Rossby numbers. Overplotted in top-right panel are contours of constant Rossby numbers (Ro = 0.01, 0.1, and 1.0; grey dashed lines) based on the empirically calibrated mass–convective overturn time relation by Wright et al. (2011). The canonical rotational I/C sequences are illustrated by black lines in top panels. We noticed that there exist similar features in cool star-spot coverages (spot filling factors), as shown in the lower panels of Fig. 7, where the spot filling factors were estimated by modelling their TiO band near 7050 Å (TiO2) using the standard values obtained based on references stars with solar metallicity (see Paper I for more details). It is clear that there exists a trend for fast rotators among K-type stars in Pleiades having larger spot coverages, and the spot coverages belong to two branches among GK-type Pleiades members. Further, we collected X-ray emission data from literature (Wright et al. 2011; Gondoin 2012) for stars in these open clusters, and found the existence of similar feature with dual sequences, as shown in Fig. 8. Figure 8. View largeDownload slide Fractional X-ray emission luminosities RX as a function of Teff for stars in Pleiades (left), and M34, Praesepe, and Hyades (right). These data were collected from literature, as labelled in the panels, where W11: Wright et al. (2011), G12: Gondoin (2012). Figure 8. View largeDownload slide Fractional X-ray emission luminosities RX as a function of Teff for stars in Pleiades (left), and M34, Praesepe, and Hyades (right). These data were collected from literature, as labelled in the panels, where W11: Wright et al. (2011), G12: Gondoin (2012). It is known that the open cluster members locate in three well-known regions in the rotation–colour diagram (Barnes 2003a), i.e. I (Interface) and C (Convective) sequence stars and gap stars between I and C sequence (e.g. see the top-left panel of Fig. 7). From Fig. 7, we can get a general conclusion that the activity level sequences correspond to the rotation sequences, namely, the more active sequence stars correspond to the C sequence stars, and the stars in less active sequence map to I sequence stars. The dynamo configuration in I sequences may be different with that in C sequences, e.g. an interface dynamo dominates in I sequence stars, and a turbulent or convective dynamo works in C sequence stars, correspondingly causing a large-scale interface magnetic field and a small-scale turbulent field (Barnes 2003a,b). Therefore, such a bifurcation of CA of Pleiades GK-type members is probably due to a change in the dynamo configuration. For older open clusters such as Hyades with an age of 600–700 Myr their GK-type member stars rotate slowly and should be in I sequence in the period–colour diagram, but the M members display clearly both I and C sequence (see the top-right panel of Fig. 7). As mentioned above, indeed there is a presence of two H α emission sequences for early and Mid-M stars in Praesepe and Hyades. Moreover, it is clear that the two activity sequences correspond to I- and C-rotation sequences, as shown in Fig. 7. For M34, an open cluster with an age of about 220 Myr, the two activity sequences should be observed among K-type stars. However, we could not notice such feature due to lack of observational data of the M34 K-type stars. 4.2.2 Rotation–activity correlation We have investigated the rotation–activity connection in terms of Rossby number Ro, however, it is difficult to judge whether Rossby number is superior to rotation period when modernizing this rotation–activity relation (e.g. see a comprehensive study on this topic by Reiners et al. 2014). In Fig. 9 the $$R^{\prime }_{\text{H}\,\alpha}$$ against Ro is displayed for our sample stars whose periods are available. As expected, it shows two different regimes, saturated and unsaturated, with different dependence of activity on rotation around a turnoff point (e.g. Ro ∼ 0.1), usually shown by X-ray emissions (e.g. Wright et al. 2011), chromospheric Ca ii and H α emissions (e.g. Noyes et al. 1984; Jackson & Jeffries 2010; Douglas et al. 2014; Newton et al. 2017), and cool spot coverages (Fang et al. 2016). Figure 9. View largeDownload slide $$R^{\prime }_{\text{H}\,\alpha}$$, $$R^{\prime }_{\text{H}\,\beta}$$, and $$R^{\prime }_{\text{Ca K}}$$ versus Ro for stars in Pleiades (left) and M 34, Praesepe, and Hyades (right). The colour gradient represents different effective temperatures. The cross symbols show slowly rotating Kepler M dwarfs (Prot > 10  d, Mass <0.6 M⊙). The fully convective candidates are marked by large circles. The grey solid and dashed line in top-right panel denotes the relation of Newton et al. (2017) and Douglas et al. (2014), respectively. Figure 9. View largeDownload slide $$R^{\prime }_{\text{H}\,\alpha}$$, $$R^{\prime }_{\text{H}\,\beta}$$, and $$R^{\prime }_{\text{Ca K}}$$ versus Ro for stars in Pleiades (left) and M 34, Praesepe, and Hyades (right). The colour gradient represents different effective temperatures. The cross symbols show slowly rotating Kepler M dwarfs (Prot > 10  d, Mass <0.6 M⊙). The fully convective candidates are marked by large circles. The grey solid and dashed line in top-right panel denotes the relation of Newton et al. (2017) and Douglas et al. (2014), respectively. The rotation–activity relationship is often parametrized by a flat region connected to a power-law decay, i.e. the activity has constant (saturated) value below a turnoff of Ro (namely Rosat), and above Rosat it declines as a power law of Ro, indicates activity is proportional to Roβ. Wright et al. (2011) found that the decline in coronal emission follow a power law with β = −2.18 for Ro > 0.13 based on large sample of 824 late-type stars. Furthermore, they derived β = −2.7 for a set of solar-type stars. However, the chromospheric emissions show a much shallower power-law decay for slow rotators, e.g. the first two Ca ii IRT excess emissions in early M-type dwarfs in NGC 2516 decrease as a power law of β ∼ −1 for Ro ≳ 0.1 (Jackson & Jeffries 2010). Recently, Douglas et al. (2014) derived much shallower power law with β = −0.73 for Ro > 0.11 in H α emissions of late K and M dwarfs of Praesepe and Hyades. More recently, Newton et al. (2017) detected a power law index β = −1.7 with Rosat = 0.21 in excess H α emission based on large sample of nearby M dwarfs. The rapid rotators with Ro ≲ 0.1 in Fig. 9 indicates that the chromospheric emissions have similar values shaping a plateau feature (e.g. Delfosse et al. 1998; Reiners, Joshi & Goldman 2012; Douglas, Agüeros & Covey 2014; Newton et al. 2017). Approximately, H α emission saturates at $$\log R^{\prime }_{\text{H}\,\alpha }\approx -3.7$$ dex, while H β and Ca ii K emission saturate respectively at $$\log R^{\prime }_{\text{H}\,\beta }\approx -4.1$$ dex and $$\log R^{\prime }_{\text{Ca K}}\approx -4.0$$ dex. The activity level is usually thought to be a constant value in the saturated regime. However, our data show that the chromospheric emissions still increase slightly as Ro decreases, a phenomenon exists in other activity indicators, e.g. Ca ii IR emissions of early M-type members of NGC 2516 (Jackson & Jeffries 2010). We found that it can be roughly fitted by some power of Ro, e.g. $$R^{\prime }_{\text{H}\,\alpha }\propto \text{Ro}^{-0.2}$$ for Pleiades stars with Ro ≲ 0.1, as shown by black dashed line in the top-left panel of Fig. 9, a relationship consistent with that (LX/Lbol ∝ Ro−0.16) derived by Reiners et al. (2014) based on X-ray emission. A similar increment could be seen in the H β emission, e.g. $$R^{\prime }_{\text{H}\,\beta }\propto \text{Ro}^{-0.28}$$ (the black dashed line in middle-left panel of Fig. 9). We could not see such trend in Ca ii K emission in this regime for total sample stars, which may be due to large uncertainties in $$R^{\prime }_{\text{Ca K}}$$ for cooler stars. However, one can see that there also exist a slight increase relation in GK-type stars, e.g. $$R^{\prime }_{\text{Ca K}}\propto \text{Ro}^{-0.2}$$ (the dashed line in bottom-left panel of Fig. 9). Such a slight slope of the activity–rotation relationship in the saturated regime indicates some unclear dependence of the activity on rotation is remaining even when saturation is reached, as pointed out by Reiners et al. (2014). The activity levels of ultrafast rotating G/K stars are decreasing below the saturation level, a phenomenon known as supersaturation which has been detected in the coronal emission (e.g. Randich et al. 1996; Wright et al. 2011). we found no clear evidence in our sample for any systematic fall in CA levels when the Rossby numbers become very small. However, we noticed that the average CA of M-type stars follows a decreasing trend as rotation rate increases at Ro ≲ 0.01 (particularly in terms of the excess Ca ii K emission), showing a weak evidence of damped cooling efficiency in ultrafast rotators. The regime with Ro ≳ 0.1 is usually thought to be unsaturated regime. At first sight, the activity decreases with increasing Ro, which is expected from qualitative arguments based on the αΩ-type shell dynamo theory (e.g. Noyes et al. 1984; Charbonneau 2014, and references therein). Moreover, the stars with Ro in the regime of 0.1–0.4 have large scatter in activity with upper envelop near to saturated activity level, which makes it difficult to explain the location of turn-off point. Even worse is the strong dependence of power law with onset of saturation, Rosat, i.e. β anticorrelates with Rosat, as seen in Douglas et al. (2014) and Newton et al. (2017). Anyway, we sought to represent the potential trends using the canonical rotation–activity relation (a saturated regime and a power-law decay) following the previous studies, being fitted by eye rather than using any χ2-methods. For Pleiades candidate members (see left-hand panels), we found that the saturation may occurs at Ro near to 0.2 rather than the canonical value of 0.1, and a power law is likely in the range of −2.0 to − 2.5, e.g. it can be represented by a relation of $$R^{\prime }_{\text{H}\,\alpha }\propto \text{Ro}^{-2.2}$$ for Ro ≳ 0.19, as illustrated by the black solid line in top-left panel of Fig. 9. We found a similar power-law decay with Rosat near to 0.2 in H β emissions for the Pleiades sample stars, as illustrated by the black solid line ($$R^{\prime }_{\text{H}\,\beta }\propto \text{Ro}^{-2.3}$$ for Ro ≳ 0.19) in middle-left panel of Fig. 9. Compared to these two Balmer emissions, the Ca ii K emissions show a slight different trend, i.e. the unsaturated regime start at a larger Ro (e.g. Rosat ≳ 0.25), but appears to be a reasonably fit by a similar power law (an example relation with β = −2.3 and Rosat = 0.27 is shown by the black solid line in bottom-left panel). The results for M34, Hyades, and Praesepe are shown in the right-hand panels of Fig. 9. For comparison purpose, in each right-hand panel, we re-plotted the black lines in corresponding left-hand panel. Among these older stars, most hotter stars are slow rotators, in particular members of Praesepe and Hyades, which are in unsaturated regime as expected, and only some M dwarfs in the Praesepe and Hyades and K-type members in M34 are fast rotators (see top panels of Fig. 7) that are occupied in the saturated regime. In the saturated regime, their excess chromospheric emissions increase as Ro decreases similar to Pleiades stars, but systematically having lower level. In the unsaturated regime, they appear to still follow the power law in a similar way as the Pleiades members, but with larger dispersion. However, it seems that Ca ii K emission saturates at a larger Rosat value, e.g. Rosat ≳ 0.3 (alternatively, these stars follow a shallower power laws at the same onset saturation in Pleiades). In top-right panel of Fig. 9, also shown is the relation of β = −0.73 with Rosat = 0.11 in dashed grey line, which was derived by Douglas et al. (2014) for M dwarfs with H α in emission in the Praesepe and Hyades. However, our data are inconsistent with their results, our H α data can be fitted by a clearly deeper power law with a larger Rosat value. The strong degeneracy between β and Rosat may be a potential contributor. Moreover, Douglas et al. (2014) sample consists of stars with r″ − K ≳ 3 (roughly later than K7) and definite H α emission, neglected potential less active members with H α in absorption (but in fact having excess emissions such as these GK-type members), which may also accounted for the notable difference between their results and ours. The relation for β = −1.7 with Rosat = 0.21 derived by Newton et al. (2017) for nearby M dwarfs also shown in solid grey line in the top-right panel. The Rosat shown by our H α data is consistent with their value, however, the power-law index is slightly deeper than their value. The transformation between the solar-type stars (with radiative cores and convective envelopes) and fully convective stars takes place roughly in stars of spectral type M3–M4. The interface layer between radiative cores and convective envelopes (known as the tachocline) is believed to play an important role in the generation of the magnetic field in solar-type stars (Spiegel & Zahn 1992). Therefore, the dynamo mechanism in the fully convective stars such as middle/late M-type stars is expected to be very different because of the absence of a tachocline. However, it is evident that slowly rotating, fully convective stars still show rotation–activity (X-ray emission) connection in the same way similar to solar-type stars (Wright & Drake 2016). Recently, Newton et al. (2017) found a single relation between excess H α emission and Ro for stars including both early- and late-type M dwarfs, supporting the previous finding. Unfortunately, most M-type stars in our sample are fast rotators (see Fig. 7), and only two slow rotators in Hyades are fully convective candidate stars (marked as large open circles, and see Table 2). To take a simple look at this issue, we picked up LAMOST spectra for 47 slowly rotating M dwarfs (Prot >10 d, mass < 0.6 M⊙) from the main-sequence Kepler targets with known rotation periods provided by McQuillan, Mazeh & Aigrain (2014), and measured their excess CA emissions, as shown by crosses in right panels of Fig. 9 (only 29 of them show excess H α emission). Note that two of them are the fully convective candidates, which were marked with large open circles (but also see Table 2). It is clear that these slowly rotating M dwarfs, including the fully convective candidates, still follow the power-law decay that show no evident difference from hotter stars, which is consistent with the finding of Newton et al. (2017), though our results show a slight different power-law decay. If this scenario is real, it then suggests that the magnetic dynamo may be similar to that in solar-type stars. Wright & Drake (2016) suggest a common magnetic dynamo in both partially and fully convective stars, supporting those models in which the dynamo originates throughout the stellar convection zone. Table 2. Slow rotating, potential fully convective stars with detected excess emissions in both H α and H β lines. Object name  SpT.  Teff  Mass  Prot  Ro  $$\log R^{\prime }_{\text{H}\,\alpha}$$  $$\log R^{\prime }_{\text{H}\,\beta}$$  $$\log R^{\prime }_{\text{Ca K}}$$  Note      (K)  (M⊙)  (d)            04202761+1853499  M3.6  3244  0.375  20.3105  0.4079  −4.3609  −4.7871  −4.5722  Hyades  04303385+1444532  M3.5  3257  0.382  18.4100  0.3771  −4.4220  −4.8098  −4.6805  Hyades  KIC 10731839  M3.3  3272  0.362  46.1150  0.8946  −5.4351  −5.7745  −5.2008  Kepler target  KIC 11495571  M3.4  3292  0.374  13.6400  0.2737  −4.0450  −4.4538  –  Kepler target  Object name  SpT.  Teff  Mass  Prot  Ro  $$\log R^{\prime }_{\text{H}\,\alpha}$$  $$\log R^{\prime }_{\text{H}\,\beta}$$  $$\log R^{\prime }_{\text{Ca K}}$$  Note      (K)  (M⊙)  (d)            04202761+1853499  M3.6  3244  0.375  20.3105  0.4079  −4.3609  −4.7871  −4.5722  Hyades  04303385+1444532  M3.5  3257  0.382  18.4100  0.3771  −4.4220  −4.8098  −4.6805  Hyades  KIC 10731839  M3.3  3272  0.362  46.1150  0.8946  −5.4351  −5.7745  −5.2008  Kepler target  KIC 11495571  M3.4  3292  0.374  13.6400  0.2737  −4.0450  −4.4538  –  Kepler target  View Large 4.3 Correlations between activity indicators The correlations among activity indicators of different temperature layers contain information of overall heating process and systematic deficiencies in these layers (Güdel 2004). Therefore, the correlations between different CA indicators and the relation of excess chromospheric emission with other activity indicators had gained lot of attention (e.g. Schrijver, Dobson & Radick 1992; Montes et al. 1995, 1996; Martínez-Arnáiz et al. 2011). In fact, this topic need further investigation, e.g. need complex physical models to interpret the correlations, which is beyond the scope of this paper. Our intent here is to merely investigate the basic observational correlations between different tracers of activity. 4.3.1 H α, H β, and Ca ii K Ca ii H&K lines are the most widely used CA indicators as their source function are collisionally dominated. H α and H β are the two strongest lines among Balmer series, and their excess emissions are formed at the middle chromosphere (Montes et al. 2004). Our data show that there exist good correlations of excess emissions of H α to that of H β and Ca ii K lines, as shown in Figs 10 and 11. Figure 10. View largeDownload slide The correlation of excess equivalent widths (upper panels) and fractional excess luminosities (lower panels) of H β and H α lines. Dashed lines in upper panels: power-law approximation (see the text). The dashed lines with grey colours in the lower panels denote constant ratios of $$R^{\prime }_{\text{H}\,\alpha }/R^{\prime }_{\text{H}\,\beta } = 0.5, 1, 2, 4,\ \text{and}\ 8$$, respectively. Figure 10. View largeDownload slide The correlation of excess equivalent widths (upper panels) and fractional excess luminosities (lower panels) of H β and H α lines. Dashed lines in upper panels: power-law approximation (see the text). The dashed lines with grey colours in the lower panels denote constant ratios of $$R^{\prime }_{\text{H}\,\alpha }/R^{\prime }_{\text{H}\,\beta } = 0.5, 1, 2, 4,\ \text{and}\ 8$$, respectively. Figure 11. View largeDownload slide Relationships between excess emissions of H α and Ca ii K lines (see the text). Figure 11. View largeDownload slide Relationships between excess emissions of H α and Ca ii K lines (see the text). Fig. 10 shows a power-law correlation (linear relation in log–log scale) between excess equivalent widths of the two Balmer lines, EW$$^{\prime }_{\text{H}\,\alpha}$$ and EW$$^{\prime }_{\text{H}\,\beta}$$, among the Pleiades stars. We found that it can be represented well by a power-law index in the range of 1.2–1.4, e.g. EW$$^{\prime }_{\text{H}\,\beta } \propto ({\text{EW}}^{\prime }_{\text{H}\,\alpha })^{1.3}$$, as illustrated by the black dashed line in top-left panel of Fig. 10. For members in M34, Praesepe, and Hyades, it also can be represented roughly by the same power-law relation, as shown by the black dashed line in top-right panel. It is noticeable that the main trend between $$R^{\prime }_{\text{H}\,\beta}$$ and $$R^{\prime }_{\text{H}\,\alpha}$$ also roughly linear in log–log scale. Moreover, we found that the power-law index is around 1.0, which indicates a linear relation roughly even in linear scale, e.g. $$R^{\prime }_{\text{H}\,\beta } \propto R^{\prime }_{\text{H}\,\alpha}$$, as demonstrated by dashed lines with constant ratios ($$R^{\prime }_{\text{H}\,\alpha }/R^{\prime }_{\text{H}\,\beta } = 0.5, 1, 2, 4, {\text{ and }} 8$$, from top to bottom, respectively) in lower panels of Fig. 10. A constant H α to H β ratio has been detected in previous studies, e.g. Reid et al. (1995) found a good linear correlation between the flux emitted at H α and H β for low-mass stars in the Hyades (FH α ∼ 4.6 × FH β). However, some late F-type or early G-type stars show departures from the general trend, e.g. their H β losses seem to be more important than H α losses. In fact, compared to H α, H β is preferred the higher chromospheric density, e.g. while H β enhanced during a flare-like event, the H α is not proportionally increased in strength (Huenemoerder & Ramsey 1987). Thus the ratio of energy emitted in the H α to H β ($${{E}}_{{\rm H}\,\,\alpha }/{{E}}_{\rm{H}\,\,\beta}$$) is widely used as indicator of the presence of flare-like events (Huenemoerder & Ramsey 1987), and also a diagnostic indicator for discriminating between stellar plages and prominences (Hall & Ramsey 1992). The typical value of $${{E}}_{\rm{H}\,\alpha }/{{E}}_{\rm{H}\,\beta}$$ is around 1 in flare-like events, it then become 4 or larger for RS CVn stars (Huenemoerder & Ramsey 1987). Therefore, the departure that some earlier type stars with lower H α to H β ratios ($$R^{\prime }_{\text{H}\,\alpha }/R^{\prime }_{\text{H}\,\beta }<1$$) were likely underlying flare/plage-like activities. If this is the real case, then it suggests that these stars preferred to suffer flare events than cooler stars. It is clear from Fig. 11 that there exist power-law relationship (linear relation in log–log scale) between excess equivalent widths of H α and Ca ii K lines. Here we did not attempt any χ2 fitting but merely sought to reproduce the correlation by eye as before. It can be represented approximately by a power-law index of 1.0–1.1, e.g. EW$$^{\prime }_{\text{Ca K}} \propto ({\text{EW}}^{\prime }_{\text{H}\,\alpha })^{1.05}$$ (see the black dashed lines in upper panels), for which the power-law exponent is slight larger than the value (0.97) derived by Montes et al. (1995) based on a sample of 51 chromospherically active binaries. The lower panels of Fig. 11 correspond to relations between excess fractional luminosities of Ca ii K and H α, which could be represented approximately by a relation of $$R^{\prime }_{\text{Ca K}}\propto ({R^{\prime }_{\text{H}\,\alpha }})^{0.5}$$ for GK-type Pleiades stars (see the black dashed lines in lower panels). However, cooler and more chromospherically active stars such as M-type stars do not obey the general relations, showing departure feature of much sharper slope compared to hotter stars. It is illustrated as an example relation $$R^{\prime }_{\text{Ca K}}\propto ({R^{\prime }_{\text{H}\,\alpha }})^{2.0}$$ by grey solid lines in the figure, which indicates that H α losses tend to be more important than Ca ii K losses for cooler stars that are more active. In fact, previous studies show such departures from the general trend in surface flux–flux relations between H α and other chromospheric indicators among some late-K and M dwarfs, and they are found to be young stars or flare stars being in the X-ray saturation regime (e.g. Martínez-Arnáiz et al. 2011). One can see that the departures hosted in M stars with non-universality trends appear both in young open cluster, Pleiades, and median-age open clusters, Hyades, and most of these stars are indeed in H α emission saturation regime. It has been suggested that the stars deviating from the general relation probably have a different magnetic structure (Martínez-Arnáiz et al. 2011), e.g. nanoflare heating. 4.3.2 H α and photospheric activity In Fig. 12, we showed the spot filling factor (or fractional spot coverage, fs) of stars with detected cool spots (fs > 0) in Pleiades as a function of $$R^{\prime }_{\text{H}\,\alpha}$$. The spot coverage of stars have large scatter, which makes harder to say about existence of any believable connections between these two activity indicators. However, the overall trend, in particular among more chromospherically active members (e.g. if we consider stars with $$\log R^{\prime }_{\text{H}\,\alpha }>-5.0$$ and keep aside the outliers at the end of weak H α emission), shows that more active stars tend to have larger spot coverages. We use a Pearson's product–moment correlation analysis to test the statistical significance of the potential connections between fs and $$\log R^{\prime }_{\text{H}\,\alpha}$$. The Pearson's correlation coefficient r ≈ 0.55 among stars with $$\log R^{\prime }_{\text{H}\alpha }>-5.0$$ indicates that there exist some correlation between spot coverage and chromospheric emission. Compared to hotter stars such as GK-type members, most M-type Pleiades members (see red symbols in Fig. 12) are in H α emission saturation regime, showing a different and slightly tighter trend in the plot, which also evidenced by a stronger positive correlation of r ≈ 0.63 (for stars with Teff < 3900 K and $$\log R^{\prime }_{\text{H}\,\alpha }>-4.0$$). It is clear that such non-universality feature is very similar to the relation between chromospheric emissions of H α and Ca ii K, as discussed above. Note that the uncertainties are larger in determination of spot coverages for M-type star as discussed in Paper I, and the correlation showed in the plot needs to be further investigated, which puts a caution for further discussion based on these features. Figure 12. View largeDownload slide Relations between $$R^{\prime }_{\text{H}\,\alpha}$$ and spot filling factor for the Pleiades sample stars with detected cool spots coverage. Colours as in Fig. 9. Figure 12. View largeDownload slide Relations between $$R^{\prime }_{\text{H}\,\alpha}$$ and spot filling factor for the Pleiades sample stars with detected cool spots coverage. Colours as in Fig. 9. The amplitude of light variation due to spot rotational modulation carries information of inhomogeneity in the stellar surface, thus acts as another indicator of photospheric activity, e.g. the light variation becomes larger with increase in CA among Sun-like stars (e.g. Lockwodd, Skiff & Radick 1997; Radick et al. 1998; Lockwodd et al. 2007) and M dwarfs (e.g. Newton et al. 2017). The amplitudes of periodic light variation versus excess H α emissions for Pleiades candidates with detected rotation period are shown in Fig. 13, where the r-band peak-to-peak photometric amplitudes Ar were derived by Hartman et al. (2010) based on the Hungarian-made Automated Telescope Network (HATNet) light curves. On the whole, there is a trend showing larger chromospheric emission corresponds to larger light variation, however, the dispersion is large as shown in the left-hand panel. The Pearson's correlation coefficient in log–log scale is r ∼ 0.71 (with a p-value of < 2.2 × 10−16). The right-hand panel of Fig. 13 is analogous to the left-hand panel, wherein the amplitude of variation (between 10 per cent and 90 per cent of the distribution of points) AK2 is derived by Rebull et al. (2016) based on K2 observations. We noticed some additional features from this panel. The correlation is weaker that is evidenced by a smaller value of r ∼ 0.47 (p ∼ 1.4 × 10−9). And there exist systematic offset between the AK2 and Ar for very cool stars (see red symbols). By neglecting the very cool stars with Teff < 3900 K, the correlation become tighter (r ∼ 0.62, p ∼ 1.9 × 10−10). Previous studies show that the short-term photometric variability (mainly due to rotational spot modulation) in less active Sun-like stars can be related to their average level of CA by a power law (e.g. Radick et al. 1998). The Pleiades GK-type members in the unsaturated regime are roughly linear on the log–log scale plot, indicating a power-law relation. On the other hand, M dwarfs still follow the general trend of more active stars showing higher levels of photometric variability (e.g. Newton et al. 2017). As most of the Pleiades M-type stars are in saturated CA regime, there is no believable trend among these very cool stars with Teff < 3900 K in both panels (left: r ∼ 0.21, p ∼ 0.12; right: r ∼ 0.08, p ∼ 0.50), which probably due to their very high activity levels. In fact, the CA levels vary in short term (as shown in Fig. 5), and long term such as cyclic variation (e.g. Baliunas et al. 1995), thus the panels of Fig. 13 actually represent a snapshot of stellar variability. Moreover, the pattern and evolution of spots, and stellar inclination (the orientation angle of the axis of rotation relative to our line of sight) could effect the light variation pattern and contaminate the potential relation between chromospheric emission and light variation, thus the feature shown in this figure is only instructive rather than conclusive. Figure 13. View largeDownload slide Amplitudes of light variation (left: HATNet, right: K2) as a function of $$\log R^{\prime }_{\text{H}\,\alpha}$$. Note that the amplitude is in log scale. Figure 13. View largeDownload slide Amplitudes of light variation (left: HATNet, right: K2) as a function of $$\log R^{\prime }_{\text{H}\,\alpha}$$. Note that the amplitude is in log scale. 4.3.3 H α and coronal X-ray emission To investigate the relation between chromospheric heating and coronal emission, we have taken RX (=LX/Lbol) from Wright et al. (2011) for common Pleiades candidate members, and plotted against $$R^{\prime }_{\text{H}\,\alpha}$$ in Fig. 14. A simple linear fit in log scale indicated that the power law index is around 1.44, e.g. $$R^{\prime }_{\text{H}\,\alpha }\propto (R^{\prime }_{\text{H}\,\alpha })^{1.44}$$ (see black dashed line in Fig. 14). A slight steeper correlation of coronal and chromospheric emission flux with a power of 1.5–1.7 based on magnetic induced excess flux in the Ca ii H&K lines was reported (Schrijver et al. 1992; Güdel 2004, and further references cited therein). Such a non-linearity in the chromospheric-coronal flux–flux correlation may be a manifestation of the chromospheric radiative losses (Güdel 2004). The correlation between H α emission and X-ray emission indicates that the formation process may be same or at least correlated. Figure 14. View largeDownload slide The relations between H α fractional luminosity $$R^{\prime }_{\text{H}\,\alpha}$$ and X-ray fractional luminosity RX. Figure 14. View largeDownload slide The relations between H α fractional luminosity $$R^{\prime }_{\text{H}\,\alpha}$$ and X-ray fractional luminosity RX. 4.4 What can we learn from the results The different onsets of saturation in various CA indicators suggest the existence of transition regime at Ro range of 0.1–0.4, which connects the unsaturated and saturated regime with steep and shallower power law, respectively. In fact, Gondoin (2012, 2013) showed the evidence for transition from saturated to non-saturated X-ray emission that occurs at Ro between ∼0.13 and 0.4 in open clusters M34 and M35. Indeed, this transition scenario could describe the scatter in Ro region of 0.1–0.4. However, such scenario fails to explain the whole observational features in this Ro regime, e.g. the Ca K emission that shows a clear rotation–activity relationship, suggesting the investigation towards other unknown parameters or the activity behaviours in this regime. In addition, the chromospheric emission follows the similar (slight shallower) power laws in the unsaturated regime compared to coronal X-ray emission, but appears to be shifted to a larger Rosat value, which indicates that chromospheric emission gets easily saturated at a given rotation rate than coronal emission. Alternatively, it may suggest that the convective overturn time-scale, empirically derived based on X-ray emission (Wright et al. 2011), is not suitable for chromospheric emission. A corresponding relation between the saturated/non-saturated regimes of X-ray emission and the C/I rotation sequences among late-type stars has been reported in literature (e.g. Barnes 2003b; Gondoin 2012, 2013). In our work, we detected a similar correspondence in H α emission among cool stars in Pleiades, Praesepe, and Hyades, as discussed in previous sections. Barnes (2003b) and Gondoin (2012) suggested such transition from the saturated to the non-saturated regime of activity is the manifestation of a dynamo transition, e.g. from a turbulent dynamo to an interface-type dynamo. Durney, De Young & Roxburgh (1993) showed a candidate turbulent dynamo, which can generate small-scale turbulent magnetic field in late-type stars, and have a weak dependence on rotation. Such a turbulent dynamo can naturally explain the weak dependence of activity on rotation among fast rotators in the saturation regime. Moreover, the generation of small-scale magnetic field supports some observational behaviour shown by very active stars, e.g. large spot coverage of stars in young open clusters may be due to many small and randomly located spots on the stellar photosphere (Jackson & Jeffries 2012, 2013; Fang et al. 2016), and the possibility of multiple small flares such as nanoflares that might be responsible for the departures of relation between various chromospheric emissions among very active late-K and M dwarfs reported by previous studies and discussed above in current work. We thus conclude that a turbulent-like dynamo dominates in fast rotating cool stars in the saturated regime. 5 CONCLUSION We have quantified CA (excess emissions in H α, H β, and Ca ii K lines) in more than 700 stars from various open clusters including Pleiades, M34, Praesepe, and Hyades, using LAMOST DR3 spectra. The CA was investigated as a function of stellar parameters: mass, age, and rotation. We also analysed the correlation between activity indicators of various atmosphere layers. CA was found to be overall increased as temperature decreases, showing a mass-dependence, which partially could be due to deeper convective zones in cooler stars. CA also shows a clear age-dependence, i.e. activity level of GK-type stars declines steadily from 100 Myr to about 700 Myr. Additionally, we found that the activity lifetime for M1–M2-type stars is around 600–700 Myr, which is consistent with previous studies. More interestingly, there appear two sequences among Pleiades GK-type stars in activity (spot coverage, chromospheric, and coronal emission)–Teff diagram, paralleling with well-known rotation sequences. We also found the existence of similar sequences in Praesepe and Hyades M-type members. We found CA is correlated with rotation in the unsaturated regime, following a power law with power of β between −2 and −2.5, which is marginally consistent with results from X-ray emissions, but with a larger turnoff point of Ro (e.g. Rosat ∼ 0.2–0.3). A weak dependence of CA on rotation is noticed in the saturated regime of activity (Ro ≲ 0.1), e.g. a slight increase as Ro decreases (a power law with β ≈ 0.2), peaks around Ro ∼ 0.01, and decrease when Ro decreases for Ro ≲ 0.01. In addition, our analysis confirms the previous finding that slowly rotating, fully convective stars also follow the activity–rotation trend like the hotter stars. We found a good linear correlation between the luminosity emitted at H α and H β. The relations from excess emission of H α to that Ca ii K lines confirm that cool stars follow two different power-law flux–flux relationships. We also found weak correlation between chromospheric emission and photospheric activity indicators (cool spot filling factor and light variation amplitude), however, it shows dependence on spectral type and activity level, which provides clues on how spot configuration varies as a function of mass and activity. A good correlation between chromospheric and coronal emissions was also detected that confirms previous findings. ACKNOWLEDGEMENTS We thank the referee for his/her observant comments and constructive suggestions that helped to improve the manuscript. This study is supported by the National Natural Science Foundation of China under grant No. 11390371, 11233004, 11573035, 11625313, 11550110492, the National Key Basic Research Program of China (973 program) 2014CB845701. This work is supported by the Astronomical Big Data Joint Research Center, co-founded by the National Astronomical Observatories, Chinese Academy of Sciences and the Alibaba Cloud. YBK is thankful to the Chinese Academy of Sciences for supporting through CAS PIFI Fellowship. This work has made use of LAMOST data. The Guo Shou Jing Telescope (the Large sky Area Multi-Object fibre Spectroscopic Telescope, LAMOST) is a National Major Scientific Project built by the Chinese Academy of Sciences. Funding for the project has been provided by the National Development and Reform Commission. LAMOST is operated and managed by National Astronomical Observatories, Chinese Academy of Sciences. This research has made use of the VizieR catalogue access tool and the cross-match service provided by CDS, Strasbourg, France. This research has made use of the WEBDA data base, operated at the Department of Theoretical Physics and Astrophysics of the Masaryk University. REFERENCES Agüeros M. A. et al.  , 2011, ApJ , 740, 110 https://doi.org/10.1088/0004-637X/740/2/110 CrossRef Search ADS   An D., Terndrup D. M., Pinsonneault M. H., Paulson D. B., Hanson R. B., Stauffer J. R., 2007, ApJ , 655, 233 https://doi.org/10.1086/509653 CrossRef Search ADS   Baliunas S. L. et al.  , 1995, ApJ , 438, 269 https://doi.org/10.1086/175072 CrossRef Search ADS   Baliunas S. L., Donahue R. A., Soon W., Henry G. W., 1998, in Donahue R. A., Bookbinder J. A., eds, ASP Conf. Ser. Vol. 154, Cool Stars, Stellar Systems, and the Sun . Astron. Soc. Pac., San Francisco, p. 153 Barnes S. A., 2003a, ApJ , 586, 464 https://doi.org/10.1086/367639 CrossRef Search ADS   Barnes S. A., 2003b, ApJ , 586, L145 https://doi.org/10.1086/374681 CrossRef Search ADS   Bouy H. et al.  , 2015, A&A , 577, 148 CrossRef Search ADS   Brandt T. D., Huang C. X., 2015, ApJ , 807, 24 https://doi.org/10.1088/0004-637X/807/1/24 CrossRef Search ADS   Canterna R., Crawford D. L., Perry C. L., 1979, PASP , 91, 263 https://doi.org/10.1086/130482 CrossRef Search ADS   Carrera R., Pancino E., 2011, A&A , 535, 30 CrossRef Search ADS   Charbonneau P., 2014, ARA&A , 52, 251 CrossRef Search ADS   Chen Y., Girardi L., Bressan A., Marigo P., Barbieri M., Kong X., 2014, MNRAS , 444, 2525 https://doi.org/10.1093/mnras/stu1605 CrossRef Search ADS   Copenhagen University Institute of Astronomy Cambridge, UK Real Instituto y Observatorio de la Armada en San Fernando, 2011, yCat , 1327, 0 Covey K. R. et al.  , 2016, ApJ , 822, 81 https://doi.org/10.3847/0004-637X/822/2/81 CrossRef Search ADS   Cram L. E., Mullan D. J., 1979, ApJ , 234, 579 https://doi.org/10.1086/157532 CrossRef Search ADS   Cui X. -Q. et al.  , 2012, RA&A , 12, 1197 https://doi.org/10.1088/1674-4527/12/9/003 Delfosse X., Forveille T., Perrier C., Mayor M., 1998, A&A , 1998, 331 , 581 Delorme P., Collier Cameron A., Hebb L., Rostron J., Lister T. A., Norton A. J., Pollacco D., West R. G., 2011, MNRAS , 413, 2218 https://doi.org/10.1111/j.1365-2966.2011.18299.x CrossRef Search ADS   Douglas S. T., Agüeros M. A., Covey K. R., 2014, ApJ , 795, 161 https://doi.org/10.1088/0004-637X/795/2/161 CrossRef Search ADS   Douglas S. T., Agüeros M. A., Covey K. R., Cargile P. A., Barclay T., Cody A., Howell S. B., Kopytova T., 2016, ApJ , 822, 47 https://doi.org/10.3847/0004-637X/822/1/47 CrossRef Search ADS   Durney B. R., De Young D. S., Roxburgh I. W., 1993, Solar Phys. , 145, 207 https://doi.org/10.1007/BF00690652 CrossRef Search ADS   Fang X.-S., Zhao G., Zhao J.-K., Chen Y.-Q., Bharat Kumar Y., 2016, MNRAS , 463, 2494 (Paper I) https://doi.org/10.1093/mnras/stw1923 (Paper I) CrossRef Search ADS   Fossati L., Bagnulo S., Landstreet J., Wade G., Kochukhov O., Monier R., Weiss W., Gebran M., 2008, A&A , 483, 891 CrossRef Search ADS   Gallet F., Bouvier J., 2013, A&A , 556, A36 CrossRef Search ADS   Geller A. M., Latham D. W., Mathieu R. D., 2015, AJ , 150, 97 https://doi.org/10.1088/0004-6256/150/3/97 CrossRef Search ADS   Gondoin P., 2012, A&A , 546, A117 CrossRef Search ADS   Gondoin P., 2013, A&A , 556, A14 CrossRef Search ADS   Gray R. O., Corbally C. J., Garrison R. F., McFadden M. T., Robinson P. E., 2003, AJ , 126, 2048 https://doi.org/10.1086/378365 CrossRef Search ADS   Gray R. O., Corbally C. J., Garrison R. F., McFadden M. T., Bubar E. J., McGahee C. E., O'Donoghue A. A., Knox E. R., 2006, AJ , 132, 161 https://doi.org/10.1086/504637 CrossRef Search ADS   Güdel M., 2004, A&AR , 12, 71 Hall J. C., 2008, Living Rev. Solar Phys ., 5, 2 https://doi.org/10.12942/lrsp-2008-2 CrossRef Search ADS   Hall J. C., Ramsey L. W., 1992, AJ , 104, 1942 https://doi.org/10.1086/116370 CrossRef Search ADS   Hartman J. D., Bakos G. Á., Kovács G., Noyes R. W., 2010, MNRAS , 408, 475 https://doi.org/10.1111/j.1365-2966.2010.17147.x CrossRef Search ADS   Hartman J. D., Bakos G. Á., Noyes R. W., Sipocz B., Kovács G., Mazeh T., Shporer A., Pál A., 2011, AJ , 141, 166 https://doi.org/10.1088/0004-6256/141/5/166 CrossRef Search ADS   Hawley S. L., Reid I. N., Gizis J. E., Byrne P. B., 1999, in Butler C. J., Doyle J. G., eds, ASP Conf. Ser. Vol. 158, Solar and Stellar Activity: Similarities and Differences . Astron. Soc. Pac., San Francisco, p. 63 Henry T. J., Soderblom D. R., Donahue R. A., Baliunas S. L., 1996, AJ , 111, 439 https://doi.org/10.1086/117796 CrossRef Search ADS   Hodgkin S. T., Jameson R. F., Steele I. A., 1995, MNRAS , 274, 869 https://doi.org/10.1093/mnras/274.3.869 CrossRef Search ADS   Huenemoerder D. P., Ramsey L. W., 1987, ApJ , 319, 392 https://doi.org/10.1086/165463 CrossRef Search ADS   Husser T.-O., Wende-von Berg S., Dreizler S., Homeier D., Reiners A., Barman T., Hauschildt P. H., 2013, A&A , 553, A6 CrossRef Search ADS   Irwin J., Aigrain S., Hodgkin S., Irwin M., Bouvier J., Clarke C., Hebb L., Moraux E., 2006, MNRAS , 370, 954 https://doi.org/10.1111/j.1365-2966.2006.10521.x CrossRef Search ADS   Jackson R. J., Jeffries R. D., 2010, MNRAS , 407, 465 https://doi.org/10.1111/j.1365-2966.2010.16917.x CrossRef Search ADS   Jackson R. J., Jeffries R. D., 2012, MNRAS , 423, 2966 https://doi.org/10.1111/j.1365-2966.2012.21119.x CrossRef Search ADS   Jackson R. J., Jeffries R. D., 2013, MNRAS , 431, 1883 https://doi.org/10.1093/mnras/stt304 CrossRef Search ADS   Jacobson H. R., Pilachowski C. A., Friel E. D., 2011, AJ , 142, 59 https://doi.org/10.1088/0004-6256/142/2/59 CrossRef Search ADS   James D. J. et al.  , 2010, A&A , 515, 100 CrossRef Search ADS   Jones B. F., Prosser C. F., 1996, AJ , 111, 1193 https://doi.org/10.1086/117865 CrossRef Search ADS   Kovács G. et al.  , 2014, MNRAS , 442, 2081 https://doi.org/10.1093/mnras/stu946 CrossRef Search ADS   Lee Y. S. et al.  , 2008, AJ , 136, 2050 https://doi.org/10.1088/0004-6256/136/5/2050 CrossRef Search ADS   Liu X.-W., Zhao G., Hou J.-L., 2015, RA&A , 15, 1089 https://doi.org/10.1088/1674-4527/15/8/001 Lockwood G. W., Skiff B. A., Radick R. R., 1997, ApJ , 485, 789 https://doi.org/10.1086/304453 CrossRef Search ADS   Lockwood G. W., Skiff B. A., Henry G. W., Henry S., Radick R. R., Baliunas S. L., Donahue R. A., Soon W., 2007, ApJS , 171, 260 https://doi.org/10.1086/516752 CrossRef Search ADS   Luo A. -L. et al.  , 2015, RA&A , 15, 1095 https://doi.org/10.1088/1674-4527/15/8/002 Mamajek E. E., Hillenbrand L. A., 2008, ApJ , 687, 1264 https://doi.org/10.1086/591785 CrossRef Search ADS   Marsden S. C., Carter B. D., Donati J.-F., 2009, MNRAS , 399, 888 https://doi.org/10.1111/j.1365-2966.2009.15319.x CrossRef Search ADS   Martínez-Arnáiz R., López-Santiago J., Crespo-Chacón I., Montes D., 2011, MNRAS , 414, 2629 https://doi.org/10.1111/j.1365-2966.2011.18584.x CrossRef Search ADS   McQuillan A., Mazeh T., Aigrain S., 2014, ApJS , 211, 24 https://doi.org/10.1088/0067-0049/211/2/24 CrossRef Search ADS   Meibom S., Mathieu R. D., Stassun K. G., Liebesny P., Saar S. H., 2011, ApJ , 733, 115 https://doi.org/10.1088/0004-637X/733/2/115 CrossRef Search ADS   Melis C., Reid M. J., Mioduszewski A. J., Stauffer J. R., Bower G. C., 2014, Science , 345, 1029 https://doi.org/10.1126/science.1256101 CrossRef Search ADS PubMed  Mermilliod J.-C., Mayor M., Udry S., 2009, A&A , 498, 949 CrossRef Search ADS   Montes D., Fernández-Figueroa M. J., de Castro E., Cornide M., 1995, A&A , 294, 165 Montes D., Fernández-Figueroa M. J., Cornide M., de Castro E., 1996, A&A , 312, 221 Montes D., Crespo-Chacón I., Gálvez M. C., Fernández-Figueroa M. J., López-Santiago J., de Castro E., Cornide M., Hernán-Obispo M., 2004, Lect. Notes Essays Astrophys. , 1, 119 Newton E. R., Irwin J., Charbonneau D., Berlind P., Calkins M. L., Mink J., 2017, ApJ , 834, 85 https://doi.org/10.3847/1538-4357/834/1/85 CrossRef Search ADS   Noyes R. W., Hartmann L. W., Baliunas S. L., Duncan D. K., Vaughan A. H., 1984, ApJ , 279, 763 https://doi.org/10.1086/161945 CrossRef Search ADS   Önehag A., Korn A., Gustafsson B., Stempels E., Vandenberg D. A., 2011, A&A , 528, A85 CrossRef Search ADS   Perryman M. A. C. et al.  , 1998, A&A , 331, 81 Pizzolato N., Maggio A., Micela G., Sciortino S., Ventura P., 2003, A&A , 397, 147 CrossRef Search ADS   Radick R. R., Thompson D. T., Lockwood G. W., Duncan D. K., Baggett W. E., 1987, ApJ , 321, 459 https://doi.org/10.1086/165645 CrossRef Search ADS   Radick R. R., Lockwood G. W., Skiff B. A., Baliunas S. L., 1998, ApJS , 118, 239 https://doi.org/10.1086/313135 CrossRef Search ADS   Randich S., Schmitt J. H. M. M., Prosser C. F., Stauffer J. R., 1996, A&A , 305, 785 Rebull L. M. et al.  , 2016, AJ , 152, 113 https://doi.org/10.3847/0004-6256/152/5/113 CrossRef Search ADS   Reid N., Hawley S. L., Mateo M., 1995, MNRAS , 272, 828 Reiners A., Joshi N., Goldman B., 2012, AJ , 143, 93 https://doi.org/10.1088/0004-6256/143/4/93 CrossRef Search ADS   Reiners A., Schüssler M., Passegger V. M., 2014, ApJ , 794, 144 https://doi.org/10.1088/0004-637X/794/2/144 CrossRef Search ADS   Salaris M., Weiss A., Percival S. M., 2004, A&A , 414, 163 CrossRef Search ADS   Scholz A., Eislöffel J., 2007, MNRAS , 381, 1638 https://doi.org/10.1111/j.1365-2966.2007.12333.x CrossRef Search ADS   Scholz A., Irwin J., Bouvier J., SipH ocz B. M., Hodgkin S., Eislöffel J., 2011, MNRAS , 413, 2595 https://doi.org/10.1111/j.1365-2966.2011.18328.x CrossRef Search ADS   Schrijver C. J., Dobson A. K., Radick R. R., 1992, A&A , 258, 432 Schuler S. C., King J. R., Fischer D. A., Soderblom D. R., Jones B. F., 2003, AJ , 125, 2085 https://doi.org/10.1086/373927 CrossRef Search ADS   Skrutskie M. F. et al.  , 2006, AJ , 131, 1163 CrossRef Search ADS   Smolinski J. P. et al.  , 2011, AJ , 141, 89 https://doi.org/10.1088/0004-6256/141/3/89 CrossRef Search ADS   Soderblom D. R., Stauffer J. R., Hudon J. D., Jones B. F., 1993, ApJS , 85, 315 https://doi.org/10.1086/191767 CrossRef Search ADS   Soderblom D. R., Jones B. F., Fischer D., 2001, ApJ , 563, 334 https://doi.org/10.1086/323829 CrossRef Search ADS   Soderblom D. R., Laskar T., Valenti J. A., Stauffer J. R., Rebull L. M., 2009, AJ , 138, 1292 https://doi.org/10.1088/0004-6256/138/5/1292 CrossRef Search ADS   Spiegel E. A., Zahn J.-P., 1992, A&A , 265, 106 Stauffer J. R., Hartmann L. W., 1986, ApJS , 61, 531 https://doi.org/10.1086/191123 CrossRef Search ADS   Stauffer J. R., Hartmann L. W., 1987, ApJ , 318, 337 https://doi.org/10.1086/165371 CrossRef Search ADS   Stauffer J. R., Giampapa M. S., Herbst W., Vincent J. M., Hartmann L. W., Stern R. A., 1991, ApJ , 374, 142 https://doi.org/10.1086/170104 CrossRef Search ADS   Stauffer J. R., Balachandran S. C., Krishnamurthi A., Pinsonneault M., Terndrup D. M., Stern R. A., 1997, ApJ , 475, 604 https://doi.org/10.1086/303567 CrossRef Search ADS   Stauffer J. R., Schultz G., Kirkpatrick J. D., 1998, ApJ , 499, L199 https://doi.org/10.1086/311379 CrossRef Search ADS   Stauffer J. R. et al.  , 2007, ApJS , 172, 663 https://doi.org/10.1086/518961 CrossRef Search ADS   Taylor B. J., 2006, AJ , 132, 2453 https://doi.org/10.1086/508610 CrossRef Search ADS   Taylor B. J., Joner M. D., 2005, ApJS , 159, 100 https://doi.org/10.1086/430422 CrossRef Search ADS   van Leeuwen F., 2009, A&A , 497, 209 CrossRef Search ADS   Vaughan A. H., Preston G. W., 1980, PASP , 92, 385 https://doi.org/10.1086/130683 CrossRef Search ADS   Walkowicz L. M., Hawley S. L., West A. A., 2004, PASP , 116, 1105 https://doi.org/10.1086/426792 CrossRef Search ADS   West A. A., Hawley S. L., 2008, PASP , 120, 1161 https://doi.org/10.1086/593024 CrossRef Search ADS   West A. A. et al.  , 2004, AJ , 128, 426 https://doi.org/10.1086/421364 CrossRef Search ADS   West A. A., Hawley S. L., Bochanski J. J., Covey K. R., Reid I. N., Dhital S., Hilton E. J., Masuda M., 2008, AJ , 135, 785 https://doi.org/10.1088/0004-6256/135/3/785 CrossRef Search ADS   West A. A., Weisenburger K. L., Irwin J., Berta-Thompson Z. K., Charbonneau D., Dittmann J., Pineda J. S., 2015, ApJ , 812, 3 https://doi.org/10.1088/0004-637X/812/1/3 CrossRef Search ADS   Wilson O. C., 1968, ApJ , 153, 221 https://doi.org/10.1086/149652 CrossRef Search ADS   Wilson O. C., 1978, ApJ , 226, 379 https://doi.org/10.1086/156618 CrossRef Search ADS   Wright N. J., Drake J. J., 2016, Nature , 535, 526 https://doi.org/10.1038/nature18638 CrossRef Search ADS PubMed  Wright N. J., Drake J. J., Mamajek E. E., Henry G. W., 2011, ApJ , 743, 48 https://doi.org/10.1088/0004-637X/743/1/48 CrossRef Search ADS   Wu Y. et al.  , 2011, RA&A , 11, 924 https://doi.org/10.1088/1674-4527/11/8/006 Yang X. L., Chen Y. Q., Zhao G., 2015, AJ , 150, 158 https://doi.org/10.1088/0004-6256/150/5/158 CrossRef Search ADS   Zacharias N., Finch C. T., Girard T. M., Henden A., Bartlett J. L., Monet D. G., Zacharias M. I., 2013, AJ , 145, 44 CrossRef Search ADS   Zhao G., Chen Y.-Q., Shi J.-R., Liang Y.-C., Hou J.-L., Chen L., Zhang H.-W., Li A.-G., 2006, Chin. J. Astron. Astrophys. , 6, 265 CrossRef Search ADS   Zhao G., Zhao Y.-H., Chu Y.-C., Jing Y.-P., Deng L.-C., 2012, RA&A , 12, 723 https://doi.org/10.1088/1674-4527/12/7/002 SUPPORTING INFORMATION Supplementary data are available at MNRAS online. online_tables.tar.gz Please note: Oxford University Press is not responsible for the content or functionality of any supporting materials supplied by the authors. Any queries (other than missing material) should be directed to the corresponding author for the article. APPENDIX A: CHROMOSPHERIC INDICATORS: Ca ii K, H β We measured equivalent widths of Ca ii K and H β lines, EWCaK, and EWHβ, and wavelength regimes for these measurements are listed in Table A1. To get reliable measurements for these two blue lines, we merely measured EWCaK for stars with g-band signal-to-noise ratio SNRg >1.0, and measured EWH β for those with SNRg >3.0. Such SNRg trim criteria cut off about 18 per cent and 3 per cent of total sample stars, for Ca ii K and H β, respectively. Fig. A1 shows the measurements of EWCaK and EWH β as a function of effective temperature. Note that the corresponding basal values are also plotted, which were derived based on a large sample of inactive reference dwarf stars with solar metallicity, following the similar procedure that used to derive mean basal values for H α (see Paper I for more details). Similar to EW$$^{\prime }_{\text{H}\,\alpha}$$, we obtained the excess equivalent widths of the Ca ii K and H β lines, EW$$^{\prime }_{\text{Ca K}}$$ and EW$$^{\prime }_{\text{H}\,\beta}$$, respectively. Further, we derived the excess fractional luminosities, $$R^{\prime }_{\text{Ca K}}$$($$\equiv L^{\prime }_{\text{Ca K}}/L_{\text{bol}}$$) and $$R^{\prime }_{\text{H}\,\beta}$$ ($$\equiv L^{\prime }_{\text{H}\,\beta }/L_{\text{bol}}$$), by using corresponding χ ratios (see Appendix C). These measurements are shown in bottom panels of Fig. A1. Figure A1. View largeDownload slide Emissions of H β and Ca ii K lines for all sample stars. Figure A1. View largeDownload slide Emissions of H β and Ca ii K lines for all sample stars. Table A1. Equivalent width measurements of H α, H β, and Ca ii K lines. Line  Line bandpass (Å)  Pseudo-continua (Å)  H α  6557–6569  6547–6557, 6570–6580  H β  4855–4867  4842–4852, 4873–4883  Ca ii K  3930–3937  3910–3915, 3950–3955  Line  Line bandpass (Å)  Pseudo-continua (Å)  H α  6557–6569  6547–6557, 6570–6580  H β  4855–4867  4842–4852, 4873–4883  Ca ii K  3930–3937  3910–3915, 3950–3955  View Large Note that the derived mean value of inactive reference stars cut off at 6500 K. However, a very small fraction of sample stars are hotter than 6500 K (see Fig. 1). We simply linearly extrapolated the reference value when necessary. In order to get reliable measurements for sample stars with Teff > 6500 K, we merely derived the excess equivalent widths for those stars with Teff less than 6900 K in this work. APPENDIX B: MEASUREMENT UNCERTAINTIES The uncertainties in the equivalent width measurements mainly result from the errors in the RV-corrected LAMOST spectra. The uncertainty in the LAMOST spectra at each wavelength (Poisson noise, and any noises during data reduction) is available and was directly propagated to the equivalent width measurements. Fig. B1 shows the errors (σEWH α) in EWH α for stars in each open cluster, and the increment in errors with decreasing temperature is seen as expected. The typical error of EWH α is less than 0.1 Å for GK-type sample stars, e.g. ∼0.05 Å for Pleiades members with temperatures around 5000 K; however, it become large for M-type stars, e.g. typically greater than 0.1 Å and up to 0.5 Å. Also, the lower SNR in the blue band of spectra increases the uncertainties in the measurements of Ca ii K and H β lines, especially for those M-type stars, e.g. the typical errors of EWH β are 0.05, 0.15, 1.0 Å, for Pleiades members with temperatures around 5500, 4500, 3500 K, respectively. Figure B1. View largeDownload slide Errors in measurements of EWH α. Note Hyades members have smaller errors because of the proximity. Figure B1. View largeDownload slide Errors in measurements of EWH α. Note Hyades members have smaller errors because of the proximity. The spectral resolution of LAMOST spectra (R ∼ 1800) limits the accuracy of RV measurements. To check how RV uncertainties affect EWH α measurements, we have taken PHOENIX model spectra and degraded the resolution to 1800 by convolving a Gaussian profile, found that the differences of EWH α due to RV shifts of ±20 km s−1 are less than 0.02 Å (a relative difference of less than 3 per cent), as shown in Fig. B2. In addition, many cool Pleiades members are fast rotators (Rebull et al. 2016), we thus evaluated the effects of rotational broadening on measurements of EWH α, and found the differences of EWH α due to a vsin i of ∼100 km s−1 are typically less than 0.04 Å (see Fig. B2), which is slightly larger than the errors due to RV uncertainties, but typically less than average errors of EWH α measurements propagated from the uncertainties of the observed flux. So the uncertainties from RV-correction and rotational broadening are marginal compared to those due to uncertainty of observed spectra. Figure B2. View largeDownload slide Top: Estimated differential EWH α values due to wavelength shifts, RV = ±20 km s−1, and rotational broadening, vsin i = 50 and 100 km s−1. Bottom: Relative differences compared to EWH α scales, note that the large values around 4300 K result from the fact that EWH α is near to zero for late K-type stars. Figure B2. View largeDownload slide Top: Estimated differential EWH α values due to wavelength shifts, RV = ±20 km s−1, and rotational broadening, vsin i = 50 and 100 km s−1. Bottom: Relative differences compared to EWH α scales, note that the large values around 4300 K result from the fact that EWH α is near to zero for late K-type stars. The uncertainties in EW$$^{\prime }_{\text{H}\,\alpha}$$ are mainly propagates from the errors of EWH α measurement and corresponding basal value used at a given temperature. The basal value of equivalent width for a given line was derived based on a large sample of reference dwarfs with solar metallicity (e.g. for H α line, see Paper I). The scatter of EWH α of reference stars in each 50 K bin is typically smaller than the errors of measured EWH α, as demonstrated in the upper panels of Fig. 2, where three times of the scatter (±3σ) are shown in grey lines. To estimate the uncertainties on the derived EW$$^{\prime }_{\text{H}\,\alpha}$$ resulting from the flux errors of spectra and uncertainties of temperature, we performed a Monte Carlo simulation 1000 times for six representative EWH α values with Teff from 6000 to 3500 K, using the typical measurement error (see Fig. B1) as the standard deviation of EWH α, and setting 100 K as the standard deviation of temperature, as shown in Fig. B3 and listed in Table B1. It is clear that the EW$$^{\prime }_{\text{H}\,\alpha}$$ uncertainties are dominated by the uncertainties of temperatures for GK-type stars. It is not the same case with M-type stars since the reference EWH α value is not so sensitive to temperature compared to hotter stars. Figure B3. View largeDownload slide Top: Examples of typical values of EWH α with measured error bars σ(EWH α). Middle and bottom: Resulted excess emissions with estimated error bars from different uncertainty sources by using Monte Carlo simulations. The data points at the same temperature are separated by 50 K for display purpose: uncertainties from EWH α (left), from Teff (centre), and from both of them (right). Figure B3. View largeDownload slide Top: Examples of typical values of EWH α with measured error bars σ(EWH α). Middle and bottom: Resulted excess emissions with estimated error bars from different uncertainty sources by using Monte Carlo simulations. The data points at the same temperature are separated by 50 K for display purpose: uncertainties from EWH α (left), from Teff (centre), and from both of them (right). Table B1. Estimated uncertainties from Monte Carlo simulations. Teff (σ)  EWH α(σ)  σ1  σ2  σ3  σ4  σ5  (K)  (Å)  (Å)  (Å)  (dex)  (dex)  (dex)  6000 (100)  −2.0 (0.03)  0.09  0.10  0.08  0.19  0.20  5500 (100)  −1.0 (0.04)  0.12  0.12  0.03  0.08  0.08  5000 (100)  0.0 (0.05)  0.09  0.10  0.02  0.04  0.04  4000 (100)  2.0 (0.10)  0.09  0.13  0.02  0.04  0.05  3500 (100)  4.0 (0.20)  0.06  0.22  0.02  0.06  0.07  4500 (100)  0.0 (0.07)  0.08  0.11  0.04  0.06  0.08  4500 (100)  1.0 (0.07)  0.08  0.11  0.02  0.03  0.04  4500 (100)  2.0 (0.07)  0.08  0.11  0.01  0.03  0.03  Teff (σ)  EWH α(σ)  σ1  σ2  σ3  σ4  σ5  (K)  (Å)  (Å)  (Å)  (dex)  (dex)  (dex)  6000 (100)  −2.0 (0.03)  0.09  0.10  0.08  0.19  0.20  5500 (100)  −1.0 (0.04)  0.12  0.12  0.03  0.08  0.08  5000 (100)  0.0 (0.05)  0.09  0.10  0.02  0.04  0.04  4000 (100)  2.0 (0.10)  0.09  0.13  0.02  0.04  0.05  3500 (100)  4.0 (0.20)  0.06  0.22  0.02  0.06  0.07  4500 (100)  0.0 (0.07)  0.08  0.11  0.04  0.06  0.08  4500 (100)  1.0 (0.07)  0.08  0.11  0.02  0.03  0.04  4500 (100)  2.0 (0.07)  0.08  0.11  0.01  0.03  0.03  σ1: $$\sigma (\text{EW}^{\prime }_{\text{H}\,\alpha}$$) due to σ(Teff). σ2: $$\sigma (\text{EW}^{\prime }_{\text{H}\,\alpha}$$) due to both σ(EWH α) and σ(Teff). σ3: $$\sigma (\log {{R}}^{\prime }_{\rm{H}\,\alpha })$$ due to σ(EWH α). σ4: $$\sigma (\log {{R}}^{\prime }_{\rm{H}\,\alpha })$$ due to σ(Teff). σ5: $$\sigma (\log {{R}}^{\prime }_{\rm{H}\,\alpha })$$ due to both σ(EWH α) and σ(Teff). View Large The excess fractional luminosity $$R^{\prime }_{\text{H}\,\alpha}$$ depends on both derived EW$$^{\prime }_{\text{H}\,\alpha}$$ and estimated value of χ at a given temperature. To check the temperature selection effect on derived $$R^{\prime }_{\text{H}\,\alpha}$$, we performed a simple Monte Carlo simulations (see Table B1), which show that temperature uncertainty dominates the derived $$R^{\prime }_{\text{H}\,\alpha}$$ uncertainty, as shown in Fig. B3. APPENDIX C: THE RATIO OF THE CONTINUUM FLUX TO THE BOLOMETRIC FLUX To quantify the strength of CA as the ratio between H α emission luminosity and bolometric luminosity, Walkowicz, Hawley & West (2004) introduced a distance-independent value χ, the ratio between the apparent continuum flux near H α, f(λ6560), and the apparent bolometric flux, fbol, namely, χ = f(λ6560)/fbol. Then the measured equivalent width of H α emission line could be converted to the fractional luminosity by using the χ value, i.e. LH α/Lbol = χ × EWH α. χ-method is widely used where flux-calibrated spectrum and distance information is not available for quantifying CA level (e.g. Douglas et al. 2014; West et al. 2015). χ values could be estimated based on observed spectra (e.g. West & Hawley 2008) or model spectra (e.g. Reiners et al. 2012). In this work, we estimated the χ values as the ratio of the surface continuum flux, F(λc), to the stellar surface bolometric flux, Fbol, using the following formula:   $$\chi = \frac{F(\lambda _{c})}{F_{\text{bol}}} = \frac{F(\lambda _{c})}{\sigma T_{\text{eff}}^4},$$ (C1)where σ is Stefan–Boltzmann constant, σ ≈ 5.6704 × 10−5 erg s−1 cm−2 K−4. We evaluated F(λc) for the lines of interest using PHOENIX ACES model spectra with solar metallicity (Husser et al. 2013). We then obtained χ values for H α, Ca ii K, and H β lines, as shown in Fig. C1. Figure C1. View largeDownload slide χ values of H α, H β, and Ca ii K lines, χH α, χH β, and χCaK, respectively, derived based on SDSS spectra, and PHOENIX model spectra with log g = 4.5 (black solid lines) and log g = 5.0 (black dashed lines). The red long dashed line denotes the relation derived from empirically determined H α line surface flux provided by Soderblom et al. (1993). Figure C1. View largeDownload slide χ values of H α, H β, and Ca ii K lines, χH α, χH β, and χCaK, respectively, derived based on SDSS spectra, and PHOENIX model spectra with log g = 4.5 (black solid lines) and log g = 5.0 (black dashed lines). The red long dashed line denotes the relation derived from empirically determined H α line surface flux provided by Soderblom et al. (1993). To check the agreement of these χ values derived from models with those from the observations, we collected the flux calibrated SDSS spectra of stars in NGC 2420, M 67, and Praesepe. 140 members of NGC 2420 with log g > 4.0 and 52 dwarfs in M 67 were selected from the catalogue of Lee et al. (2008). We estimated their radii using the temperatures from Lee et al. (2008), adopting the PARSEC model (Chen et al. 2014) with a metallicity of [Fe/H] ∼ −0.2 (Jacobson, Pilachowski & Friel 2011) and an age of ∼2 Gyr (Salaris, Weiss & Percival 2004), and the PARSEC model with a solar metallicity (Jacobson et al. 2011) and an age of ∼4 Gyr (Önehag et al. 2011; Geller, Latham & Mathieu 2015), for NGC 2420 and M 67, respectively. We collected SDSS spectra with good quality (SNR > 10) for 36 Praesepe stars of Douglas et al. (2014) based on the SDSS DR12 SSPP output catalogue (Lee et al. 2008; Smolinski et al. 2011). We estimated their effective temperatures using the colours (V − Ks and/or r − Ks), and derived corresponding radii adopting the PARSEC model with a metallicity of [Fe/H] ∼ +0.16 (Carrera & Pancino 2011) and an age of about 650 Myr. The distances adopted from WEBDA (http://www.univie.ac.at/webda/) for NGC 2420, M 67, and Praesepe are 3085, 908, and 187 pc, respectively. Based on the observed continuum flux, f(λc), and their corresponding radii (R) and distances (D), the χ ratios could be derived using the following relation:   $$\chi = \frac{F(\lambda _{c})}{F_{\text{bol}}} = \frac{(D/R)^{2}f(\lambda _{c})}{\sigma T_{\text{eff}}^4},$$ (C2)as shown in Fig. C1, wherein the error bar is equivalent to its uncertainty due to 10 per cent error of distance. There exists a slight offset between members in NGC 2420 and M 67 at the same effective temperatures, which may be due to the uncertainties in distance, and/or the difference in basic physical parameters between these two clusters. Also shown in Fig. C1 by the red dashed line is a relation for H α line, which was derived from empirically determined surface flux of H α line provided by Soderblom et al. (1993). As shown by the figure, the χ values from model spectra are overall in good agreement with those from observations, there however exist some small discrepancies between them, e.g. ∼0.1 dex in H α line around 4500 K. APPENDIX D: ESTIMATE OF QUIESCENT PHOTOSPHERE TEMPERATURE AND MASS The cool spots on the stellar surface contribute light that would affect spectral features and broad-band colours, for which the derived effective temperature becomes the average surface temperature (flux-weighted mean temperature) in spotted stars. Therefore, no direct way to derive a quiescent temperature for active stars, particularly those Pleiades members with large spot coverages (see Paper I for more discussion on this topic). For Pleiades members, we followed the procedures in Paper I to estimate the quiescent temperature (see Paper I for more details). For M34 members, we collected V-band magnitudes from literature (Jones & Prosser 1996; Irwin et al. 2006; Meibom et al. 2011), and got their Ks-band magnitudes from 2MASS archive (Skrutskie et al. 2006) via Vizier x-match service. We estimated their quiescent temperatures based on V − Ks colours using PARSEC isochrones (Chen et al. 2014) with an age of 220 Myr and solar metallicity, wherein we adopted the reddening of E(B − V) = 0.07 mag (Canterna, Crawford & Perry 1979). For Praesepe and Hyades candidates, first we derived their effective temperatures based on V − Ks colours, where V magnitudes taken from the fourth U.S. Naval Observatory CCD Astrograph Catalogue (Zacharias et al. 2013). A small fraction of sample stars have no V magnitudes, we then derived their temperatures based on r − Ks, where r magnitudes from the Carlsberg Meridian Catalogue 15 (Copenhagen University et al. 2011). The PARSEC models with an age of 650 Myr and a metallicity of [Fe/H] ∼ 0.1 were used to convert colours to temperatures, where the reddening $$E(B-V) = 0{^{\rm m}_{.}}027$$ and $$E(B-V)=0{^{\rm m}_{.}}001$$ (Taylor 2006) were adopted for Praesepe and Hyades, respectively. We then compared the colour-based temperatures with the spectra-based temperatures (Tsp, for FGK-type stars, we adopt the values provided by LASP, for M-type stars, we derived it via the molecular band CaH using the method in Paper I), and found that there exist complex differences between colour-based temperatures and spectra-based temperatures, as shown in Fig. D1. We finally adopted the effective temperatures from spectral features as the quiescent temperatures for stars in these two clusters. Figure D1. View largeDownload slide The differences between V − Ks-based temperature (TVK), r − Ks-based temperature (TrK), and the temperature from spectral features (Tsp) among stars in Praesepe and Hyades. Figure D1. View largeDownload slide The differences between V − Ks-based temperature (TVK), r − Ks-based temperature (TrK), and the temperature from spectral features (Tsp) among stars in Praesepe and Hyades. For the M-type slow rotators in Kepler field, first we estimated temperature for each target based on the molecular bands of CaH and TiO following the method in Paper I, and then took the average of our measured value and the value provided by McQuillan et al. (2014), as finally adopted quiescent temperature. To identify potential binaries and non-members for the member candidates in open clusters, we followed the procedures used in Paper I by using available CMDs, e.g. we identified these stars as photometrically binaries/non-members that have larger deviation compared with empirical single-star loci in V versus V − Ks diagram, namely, brighter stars with $$-1{^{\rm m}_{.}}0 < \Delta V < -0{^{\rm m}_{.}}5$$ were classified as probable binary members, while the fainter stars with $$\Delta V >0{^{\rm m}_{.}}3$$ or very much brighter stars with $$\Delta V < -1{^{\rm m}_{.}}0$$ were identified as probable non-member stars. Once the effective temperatures are obtained for each star, we estimated the mass from temperature based on mass–Teff relation from PARSEC isochrones (Chen et al. 2014). Note that, for these M-stars in Kepler field, we used the isochrones with a solar metallicity and an age of 2 Gyr. APPENDIX E: ROTATION PERIODS COLLECTION For Pleiades, Hartman et al. (2010) detected rotation periods for nearly 400 members using the HATNet transit survey data. Recent studies based on dedicated photometric surveys, e.g. the Palomar Transient Factory (PTF) (Covey et al. 2016), and the K2 (Rebull et al. 2016), increased the measurements of rotation periods for more than 500 new members of this open cluster, vastly expanding the number of Pleiades members with rotation periods. For these members with multirotation measurements, we adopted values from the K2 (Rebull et al. 2016), and for others, we adopted their rotation periods provided by Covey et al. (2016) and Hartman et al. (2010). Totally, 231 Pleiades members in our sample have available rotation periods. For M34, we collected rotation periods from literature, Irwin et al. (2006), James et al. (2010), and Meibom et al. (2011). The latest measurements were adopted in this paper. Totally, we found 21 members have rotation periods. For stars in Praesepe, their rotation periods were collected from Scholz & Eislöffel (2007), Delorme et al. (2011), Scholz et al. (2011), Agüeros et al. (2011), and Kovács et al. (2014). Totally, we got their rotation periods for 76 stars in Praesepe sample. For Hyades, we collected the rotation periods from Radick et al. (1987), Delorme et al. (2011), Hartman et al. (2011), and Douglas et al. (2016). Totally, 24 stars in Hyades sample have rotation periods. APPENDIX F: ON-LINE TABLES We compiled the measurements for stars in Pleiades, M34, Praesepe, and Hyades in on-line Table F1–F4, respectively, including effective (quiescent) temperatures, equivalent widths of H α, H β, and Ca ii K lines, their excess equivalent widths, excess fractional luminosities, and spot filling factors. © 2018 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Monthly Notices of the Royal Astronomical Society Oxford University Press

# Stellar activity with LAMOST – II. Chromospheric activity in open clusters

, Volume 476 (1) – May 1, 2018
19 pages

/lp/ou_press/stellar-activity-with-lamost-ii-chromospheric-activity-in-open-6a4vVD0BG1
Publisher
The Royal Astronomical Society
ISSN
0035-8711
eISSN
1365-2966
D.O.I.
10.1093/mnras/sty212
Publisher site
See Article on Publisher Site

### Abstract

Abstract We use the LAMOST spectra of member stars in Pleiades, M34, Praesepe, and Hyades to study how chromospheric activity varies as a function of mass and rotation at different age. We measured excess equivalent widths of H α, H β, and Ca ii K based on estimated chromospheric contributions from old and inactive field dwarfs, and excess luminosities are obtained by normalizing bolometric luminosity, for more than 700 late-type stars in these open clusters. Results indicate two activity sequences in cool spot coverage and H α excess emission among GK dwarfs in Pleiades and M dwarfs in Praesepe and Hyades, paralleling with well-known rotation sequences. A weak dependence of chromospheric emission on rotation exists among ultrafast rotators in saturated regime with Rossby number Ro ≲ 0.1. In the unsaturated regime, chromospheric and coronal emission show similar dependence on Ro, but with a shift towards larger Ro, indicating chromospheric emission gets easily saturated than coronal emission, and/or convective turnover time-scales based on X-ray data do not work well with chromospheric emission. More interestingly, our analysis shows fully convective slow rotators obey the rotation–chromospheric activity relation similar to hotter stars, confirming the previous finding. We found correlations among H α, H β, and Ca ii K emissions, in which H α losses are more important than Ca ii K for cooler and more active stars. In addition, a weak correlation is seen between chromospheric emission and photospheric activity that shows dependence on stellar spectral type and activity level, which provides some clues on how spot configuration varies as a function of mass and activity level. stars: activity, stars: chromospheres, stars: late-type, stars: rotation 1 INTRODUCTION Stellar magnetic activity, triggered by magnetic fields and maintained by a magnetic dynamo, is an ubiquitous phenomenon in late-type stars. Such activity includes a variety of phenomena, e.g. appearance of cool star-spots that are thought to be the fingerprints of emergent magnetic lines on the photosphere; the chromospheric and coronal heating that are responsible for excess chromospheric emissions and thermal X-rays, respectively; impulsive events like flares; and other solar phenomena. Stellar mass and rotation are thought to be two main basic parameters that affect the dynamo configuration and thus the properties of stellar magnetic activity. On the other hand, the loss of stellar angular momentum through the magnetized stellar wind depends on the magnetic field configuration on stellar surface (Gallet & Bouvier 2013). As the rotation of the star is slowed down due to angular momentum loss, its magnetic field weakens (or the efficiency of the dynamo decreases), and thus activity decays. The correlations between rotation, activity, and age are evident from observations of various tracers of activity, e.g. the rotation–activity connection is found based on stellar X-ray emissions (e.g. Pizzolato et al. 2003; Mamajek & Hillenbrand 2008; Wright et al. 2011; Reiners, Schüssler & Passegger 2014), chromospheric Ca ii HK emissions (e.g. Noyes et al. 1984; Mamajek & Hillenbrand 2008) and H α emissions (e.g. Stauffer et al. 1997; Douglas, Agüeros & Covey 2014; Newton et al. 2017), and photospheric activity indicator such as cool spot fractional coverage (Fang et al. 2016, hereafter Paper I). The connection between rotation and activity among cool stars directly allows us to understand the nature of the magnetic dynamo in these stars, however, such connection appears to be sensitive to tracers of activity, e.g. chromospheric emission shows different dependence on rotation compared to coronal X-ray emission (Douglas et al. 2014). Chromospheric activity (hereafter CA) compasses a wide range of phenomena that produce excess emission (e.g. Ca ii and Balmer lines) with respect to a radiative equilibrium atmosphere, stemming largely from changes in their magnetic field (see a comprehensive review on this topic by Hall 2008). Thanks to classic work of Mount Wilson Observatory Ca ii HK long-term observations that are operated from 1966 through 2003, for the initial purpose of searching stellar analogues of the solar cycle (Wilson 1968), wherein variety of CA properties including chromospheric variations are explored (Wilson 1978; Baliunas et al. 1995, 1998), e.g. the well-known CA long-term variation patterns (cyclic variation, irregular variation, little variation). What followed was a series of attempts trying to gain the general understanding (statistical characters) of CA among nearby late-type stars, e.g. statistical studies based on chromospheric emissions at the core of Ca ii H and K lines among nearby FGK field stars (Henry et al. 1996; Gray 2003, 2006), confirmed the bi-modal distribution in stellar activity (Vaughan & Preston 1980); West et al. (2004, 2008) investigate the fraction of active stars as a function of spectral type and the vertical distance from the Galactic plane, based on a large sample of nearby M-type stars using H α emission line as an activity indicator. Compared to field stars, open cluster members have an advantage for studying CA properties in light of rotation–activity–age connection due to their homogeneity in age and chemical composition. Therefore, CA in several open clusters got special attention, e.g. the young open cluster Pleiades (Soderblom et al. 1993) and NGC 2516 (Jackson & Jeffries 2010), the median-age open cluster Hyades (Stauffer et al. 1991, 1997) and Praesepe (Douglas et al. 2014). Based on these investigations a general picture is drawn that CA also varies as a function of both stellar mass and rotation. However, due to the limitation of sample stars used in previous studies in terms of mass coverage and available rotation periods (e.g. Pleiades), the behaviours of CA with respect to stellar mass changes and rotation is still not well understood, which demands further investigation on this issue. LAMOST (Cui et al. 2012) observed millions of stellar spectra for cool stars in our Milky Way (Liu, Zhao & Hou 2015), and more and more stars’ rotation periods become available thanks to recent and on-going photometric surveys, e.g. for ∼1000 Pleiades members (Hartman et al. 2010; Rebull et al. 2016), both of which allow us to initiate systematic study of CA among low-mass stars in different population. In Paper I, we found that almost all Pleiades members show excess emission in various indicators (e.g. Ca ii HK, H α, H β, Ca i λ4226, Mg i b triplet, Na i D lines, and Ca ii IRT lines, see figs 7 and 8 in Paper I) relative to old and inactive field stars, indicating they are chromospherically active. In this paper, we continue to investigate their CA properties in detail among members of nearby open clusters including Pleiades, M34, Praesepe, and Hyades by focusing on three chromospheric indicators: H α, H β, and Ca ii K lines. We examined their incidence and strength of CA as a function of effective temperature, and checked the correlation between activity and rotation. We also investigated the relationship amongst these three chromospheric indicators, and the relations between chromospheric emission and several manifestations of magnetic activity in different stellar atmosphere layers (cool spot coverage, light variation amplitude, and coronal emission), which holds key clues not only to the configuration of the magnetically active region on the stellar surface but also to the magnetically heating processes in stellar outer atmosphere. 2 DATA AND SAMPLE 2.1 LAMOST DR3 The LAMOST, characterized by both wide field of view (5 deg in diameter) and large effective aperture of ∼4 m, is a reflecting Schmidt telescope located at the Xinglong Observatory, China. LAMOST spectra are obtained in low resolution R ≈ 1800 (e.g. ∼3.6 Å around 6500 Å) (see Zhao et al. 2006, 2012; Luo et al. 2015, for more details). The LAMOST Date Release 3 (DR3) was made available to the Chinese astronomical community and international partners during 2015 December, containing nearly 5.3 million stellar spectra that were collected from pilot survey through 3 yr regular survey (autumn 2011 to summer 2015). LAMOST stellar parameter pipeline (LASP) provides the stellar parameters like effective temperature (Teff), surface gravity (log g), metallicity ([Fe/H]), and radial velocity (RV) for ∼3.2 million stellar spectra of AFGK stars whose spectra meet the signal-to-noise ratio criterion (e.g. g-band SNR of SNRg ≥ 15 and SNRg ≥ 6 for the bright and dark nights, respectively) (Wu et al. 2011; Luo et al. 2015). In addition, DR3 also includes a catalogue of ∼0.3 million M-type stars. 2.2 Sample stars We have selected four nearby open clusters with a range of ages from 100 to 700 Myr (Pleiades, M34, Praesepe, and Hyades), and retrieved the available spectra of good quality from LAMOST DR3 archive. The Pleiades (seven sisters) is a young open cluster with an age of ∼125 Myr (Stauffer, Schultz & Kirkpatrick 1998) and metallicity of [Fe/H] ∼ 0.03 (Soderblom et al. 2009), at a distance of around 130 pc (e.g. ∼136 pc; Melis et al. 2014). In addition to LAMOST DR2 archive used in Paper I, DR3 provides spectra for more than 70 targets from the Pleiades potential member catalogues provided by Stauffer et al. (2007) and Bouy et al. (2015) (whose probability is more than 75 per cent). In total, we selected 378 probable Pleiades members with 458 LAMOST spectra having signal-to-noise ratio in r-band SNR r > 10. M34 (NGC 1039) has an age of about 220 Myr (see Meibom et al. 2011, and references therein), located at a distance of ∼470 pc towards the Galactic anticentre (Jones & Prosser 1996). M34 has metallicity close to solar, [Fe/H] ∼ 0.07 (Schuler et al. 2003). Jones & Prosser (1996) provided proper motion membership probability analysis for more than 600 stars with brightness down to V ∼ 16.2 in the M34 field, wherein, the cooler members tend to have lower membership probability due to the faintness, e.g. the typical probability of faint stars with V > 15 (∼K2 and cooler) is ∼40 per cent, while the brighter members have probabilities close to 100 per cent. We selected candidate M34 members with membership probability greater than 40 per cent. However, real cool members might be excluded in this criterion. To collect more member stars, in particular at cooler end, we chose periodic variables detected in M34 (Irwin et al. 2006; James et al. 2010; Meibom et al. 2011), considering them as younger and more active. Unfortunately, only a small number of cool stars in the vicinity of M34 field have been observed by LAMOST (DR3), partly due to their faintness. In total, we collected 45 candidates that have LAMOST spectra with SNR r > 10, most of which are FG-type stars. The Hyades is the nearest open cluster at a distance of about 46 pc (Perryman et al. 1998), and has solar metallicity [Fe/H] ∼ 0.1 (Taylor & Joner 2005; Carrera & Pancino 2011). The Praesepe (NGC 2632; the Beehive Cluster) has a distance of 180–190 pc (e.g. An et al. 2007; van Leeuwen 2009) and a metallicity of [Fe/H] ∼ 0.16 (Carrera & Pancino 2011; Yang, Chen & Zhao 2016). The Hyades and Praesepe are both part of the Hyades supercluster, and seem to share an intermediate-age around 600–800 Myr (Perryman et al. 1998; Fossati et al. 2008; Brandt & Huang 2015). We adopt the membership catalogue provided by Douglas et al. (2014) for these two clusters, and retrieved 197 spectra for 142 Hyades candidates, and 227 spectra for 176 candidates of Praesepe from LAMOST DR3 archive with SNR r > 10. Fig. 1 and Table 1 show the summary of total sample (741 stars), and it is easy to note that more than half (∼60 per cent) of them are late K- and M-type stars. Figure 1. View largeDownload slide Distribution of the sample stars as a function of temperature. The dashed black line shows the distribution for total sample stars, while the grey line represents these stars for which rotation period is now known. The solid lines with different colours correspond to cumulative fraction profiles (right y-axis). Figure 1. View largeDownload slide Distribution of the sample stars as a function of temperature. The dashed black line shows the distribution for total sample stars, while the grey line represents these stars for which rotation period is now known. The solid lines with different colours correspond to cumulative fraction profiles (right y-axis). Table 1. Basic information about the sample stars in current work. Cluster  d  Age  [Fe/H]  E(B − V)  Samplea    (pc)  (Myr)    (mag)    Pleiades  136  ∼125  +0.03  0.03b  378 (458), 231  M34  470  ∼220  +0.07  0.07  45 (49), 21  Praesepe  184  ∼650  +0.16  0.027  176 (227), 76  Hyades  46  ∼650  +0.11  0.001  142 (197), 24  Cluster  d  Age  [Fe/H]  E(B − V)  Samplea    (pc)  (Myr)    (mag)    Pleiades  136  ∼125  +0.03  0.03b  378 (458), 231  M34  470  ∼220  +0.07  0.07  45 (49), 21  Praesepe  184  ∼650  +0.16  0.027  176 (227), 76  Hyades  46  ∼650  +0.11  0.001  142 (197), 24  aTotal number of stars (spectra), and those with rotation period. bThere is significant differential extinction in the Pleiades, see Soderblom et al. (1993). View Large 3 QUANTIFYING CHROMOSPHERIC ACTIVITY It is well known that enhanced emission of H α, produced by collisional excitation in the relatively dense chromospheres of low-mass dwarf stars, is the strongest and widely used activity indicator (e.g. Stauffer & Hartmann 1986; Soderblom et al. 1993; West et al. 2004, 2015; Newton et al. 2017). In this work, H α line was used as main CA indicator, considering that more than 50 per cent of our sample stars are late K- and M-type stars (see Fig. 1). However, we also investigated the behaviours of other two activity indicators, Ca ii K and H β lines. 3.1 Excess equivalent width The equivalent widths were measured from RV-corrected spectra of all sample stars. For FGK dwarfs, we adopted RV values derived by LASP, while we measured RVs for M dwarfs following methodology described in Paper I. The equivalent widths of H α line (hereafter EWH α) were measured using the formula (equation 4) in Paper I, where the line was centred at 6563 Å over 12 Å bandpasses and continuum flux was taken to be the average flux between 6547–6557 and 6570–6580 Å. The measurements of EWH α for entire sample are plotted against effective (quiescent) temperatures (see Appendix D for details) in Fig. 2, wherein the mean values of EWH α for inactive reference stars (a large sample of very less active field FGKM dwarfs with solar metallicities over the temperature range 3000–6500 K, see Paper I for more details) are treated as the minimum or basal value of chromospheric emissions, and shown as black solid lines. It is clear from the figure that the members of Pleiades show higher EWH α compared to inactive counterparts. In fact, the correlation between EWH α and CA level is complex, due to presence of non-monotonic behaviour (Cram & Mullan 1979), e.g. first the H α absorption increases with increase in activity, then fill-in, and eventually appear as emission line. Therefore, among the less active stars at the same temperature, the stronger Hα absorption line indicates its higher activity. In this work, we adapted simple approach by ignoring such non-monotonic effect, and consider excess EWH α above the basal value as the tracer of pure chromospheric emission. To determine the excess H α emission, we calculated the excess equivalent widths (the ΔEWH α in Paper I, hereafter denoted as EW$$^{\prime }_{\text{H}\,\alpha}$$) by subtracting the basal value at the same temperature, EWH α, basal, i.e.   $$\text{EW}^{\prime }_{\text{H}\,\alpha } = \text{EW}_{\text{H}\,\alpha } - \text{EW}_{\text{H}\,\alpha ,\text{basal}}.$$ (1)Similarly, we obtained the excess equivalent widths, EW$$^{\prime }_{\text{Ca K}}$$ and EW$$^{\prime }_{\text{H}\,\beta}$$, for Ca ii K and H β lines, respectively (see Appendix A for details). Figure 2. View largeDownload slide Top: The equivalent widths of H α line for stars in Pleiades (left), and M 34, Praesepe, and Hyades (right). Black solid lines show the mean relations for inactive reference stars, upper and lower 3σ to the mean value are shown as grey lines. Filled symbols represent probable single members and open symbols denote probable binary members or non-members. For stars having LAMOST multi-epoch observations, we displayed their average values rather than individual observational values. Middle: The excess equivalent width of H α, EW$$^{\prime }_{\text{H}\,\alpha}$$. Bottom: The excess fractional luminosity of H α, $$R^{\prime }_{\text{H}\,\alpha}$$, where the grey solid (dashed) lines denote the stars locating in the upper 3σ (σ) regime. Figure 2. View largeDownload slide Top: The equivalent widths of H α line for stars in Pleiades (left), and M 34, Praesepe, and Hyades (right). Black solid lines show the mean relations for inactive reference stars, upper and lower 3σ to the mean value are shown as grey lines. Filled symbols represent probable single members and open symbols denote probable binary members or non-members. For stars having LAMOST multi-epoch observations, we displayed their average values rather than individual observational values. Middle: The excess equivalent width of H α, EW$$^{\prime }_{\text{H}\,\alpha}$$. Bottom: The excess fractional luminosity of H α, $$R^{\prime }_{\text{H}\,\alpha}$$, where the grey solid (dashed) lines denote the stars locating in the upper 3σ (σ) regime. 3.2 Excess fractional luminosity As shown in Fig. 2, the strong change of EWH α (and EW$$^{\prime }_{\text{H}\,\alpha}$$) with temperature does reflect not only the intrinsic change of CA, but also the effects of the drop in the nearby continuum level. To remove the dependence of the measured EWH α on the surrounding continuum level, and to compare the strength of CA in different spectral types, we used the fractional luminosity of H α (normalized with the bolometric luminosity) (e.g. Reid, Hawley & Mateo 1995; West et al. 2004). However, quantifying fractional luminosity for a given line requires a flux-calibrated spectrum and the distance. For all activity indicators in this work, the excess fractional luminosities were derived from the excess equivalent widths using a distance-independent value, χ, the ratio between the continuum flux near the line of interest and the bolometric flux, following commonly used χ-methods (e.g. West et al. 2004). Using the χ ratios of H α line (χH α) obtained from the model spectra (see Appendix C), we obtained the H α excess fractional luminosity (denote as $$R^{\prime }_{\text{H}\,\alpha}$$), the ratio of excess H α emission luminosity ($$L^{\prime }_{\text{H}\,\alpha}$$) to the bolometric luminosity (Lbol), as follows:   $$R^{\prime }_{\text{H}\,\alpha } \equiv \frac{L^{\prime }_{\text{H}\,\alpha }}{L_{\text{bol}}}= {\chi }_{\text{H}\,\alpha } \times \text{EW}^{\prime }_{\text{H}\,\alpha }.$$ (2)In similar way, we estimated Ca ii K and H β excess fractional luminosities, $$R^{\prime }_{\text{Ca K}}$$ and $$R^{\prime }_{\text{H}\,\beta}$$, respectively (see Appendix A). These values provide the fraction of respective bolometric luminosity emitted from the chromosphere in the lines, and are analogous to the chromospheric Ca ii emission ratio $$R^{\prime }_{\text{HK}}$$ (e.g. Noyes et al. 1984) and $$R^{\prime }_{\text{IRT}}$$ (e.g. Marsden, Carter & Donati 2009; Jackson & Jeffries 2010). 3.3 Comparison with literature measurements We compared our EWH α measurements with Douglas et al. (2014), wherein they measured EWH α from medium resolution (about two times higher than LAMOST) spectra, for common stars in Hyades and Praesepe (see Fig. 3), and found to be in good agreement. Furthermore, Soderblom et al. (1993) and Soderblom, Jones & Fischer (2001) obtained excess H α emission from high-resolution echelle spectra for about 100 FGK dwarfs in the Pleiades and about 50 FGK dwarfs in the M34, respectively, using a synthetic subtraction technique, i.e. the H α profile of each target star is subtracted by the rotationally broadened H α of inactive counterpart. In Fig. 4 we compared their measurements of excess H α emission with our measurements for common stars in the Pleiades and M34, and found they are in close agreement. However, it seems there still exist systematic offset in fractional luminosity for less active hotter stars with $$\log R^{\prime }_{\text{H}\,\alpha } \lesssim -4.5$$, e.g. our measured $$R^{\prime }_{\text{H}\,\alpha}$$ are ∼0.2 dex larger. In fact, as discussed in Appendix B, a small error in excess equivalent width would result in larger error in fractional luminosity (e.g. see Table B1), particularly when the excess emission is very weak. In a careful check on the excess equivalent width, we indeed found there is a small offset (∼0.1 Å, see the upper panel of Fig. 4) that is comparable with the uncertainties due to 100 K uncertainty in temperature (see Table B1). Figure 3. View largeDownload slide Comparison of our measured EWH α with those of Douglas et al. (2014) for common stars in Hyades and Praesepe. Figure 3. View largeDownload slide Comparison of our measured EWH α with those of Douglas et al. (2014) for common stars in Hyades and Praesepe. Figure 4. View largeDownload slide Comparison of H α excess emissions with those of Soderblom et al. (1993) (Pleiades members) and Soderblom et al. (2001) (M34 members) for common stars. Figure 4. View largeDownload slide Comparison of H α excess emissions with those of Soderblom et al. (1993) (Pleiades members) and Soderblom et al. (2001) (M34 members) for common stars. 3.4 Variability in chromospheric emission The larger scatter in Figs 3 and 4 are of course partly from the uncertainties of our measurements. However, we believe that they are mainly related to intrinsic variation in activity, which is well-known for Pleiades members (e.g. Soderblom et al. 1993). To further check the variation in H α emission, we have plotted EWH α against the temperature in Fig. 5, in which we show minimum and maximum values of EWH α for stars having multi-epoch spectra. From the figure it is clear that very cool members show high variation in H α emission, which is evidently larger than H α equivalent width measurement errors, indicating intrinsic variability of CA. Also, it shows that more active stars have large equivalent width variations, indicating stars with high level activity suffer more often large variations, which may be due to short-time-scale evolution of active regions or flare-like events. Figure 5. View largeDownload slide Variation of EWH α for stars having multi-epoch observations. Also shown by grey-colour solid line are the basal values. Figure 5. View largeDownload slide Variation of EWH α for stars having multi-epoch observations. Also shown by grey-colour solid line are the basal values. 4 RESULTS AND DISCUSSION 4.1 Chromospheric activity as a function of mass As shown in top-left panel of Fig. 2, the EWH α is generally a strong function of stellar effective temperature, i.e. EWH α increases with decreasing temperature, from H α absorption at Teff ∼ 6500 K (EWH α ∼ −2 Å) to H α emission at Teff ∼ 3000 K (EWH α ∼ 10 Å). Such trend in the Pleiades was also observed by Stauffer & Hartmann (1987). The upper envelop of EWH α versus Teff profile indicates the activity levels of most active stars at a given temperature, and provides the clues on chromospherically cooling effects for stars with different mass. Also, it shows the dependence of the upper envelop on spectral types among Pleiades members, e.g. the EWH α increases moderately as temperature decreases for late F- and G-type stars (Teff ≳ 5000 K), and for K-type stars it goes into a plateau where exists a slight increment of EWHα, and then it increases sharply when Teff ≲ 4000 K. However, at the same time, the spread in EW H α also increases. Such feature is more overt in the middle-left panel of Fig. 2. The lower EWH α of Pleiades members follows the mean relations of inactive reference stars, which varies with temperature, i.e. from EWHα ∼ −2.5 Å at Teff ∼ 6500 K to EWH α ∼ 0 Å at Teff ∼ 3800 K. We noticed that some Pleiades G-type members have larger EWH α values compared to inactive counterparts, indicating they have shallower H α absorption line profile (i.e. the core of the line is filled in). The excess H α emission for all clusters are shown in middle panels of Fig. 2. It is clear from the figure that most of Pleiades candidates show excess emission compared to inactive counterparts, i.e. EW$$^{\prime }_{\text{H}\,\alpha }>0$$, indicating they are in chromospherically active phase, which is consistent with its young age. However, a small fraction of Pleiades stars have weaker excess H α emission compared to other stars in the same spectral type, e.g. about 40 very low-mass stars with Teff < 4000 K have EW$$^{\prime }_{\text{H}\,\alpha }<0.5$$ Å, while other members in this temperature regime have evident larger EWs (e.g. EW$$^{\prime }_{\text{H}\,\alpha }>2$$ Å), making a gap around Teff ∼ 3500 K in the diagram. One general interpretation is that these stars may be field dwarfs rather than Pleiades members. We found that about one -third of these stars are photometrically single candidate members according to their location on CMDs. To further check their membership, we plotted their RVs distribution in Fig. 6, and found about half of them have RVs from −8 to 20 km s−1 suggesting they are likely Pleiades members within 3σ-RV regime by adopting the mean RV value of μ ∼6 km s−1 and a scatter of σ ∼5 km s−1 for Pleiades (see Fig. 6). Several stars (e.g. about one-sixth of them) are probably Pleiades members which pass both the CMD and RV criterion, such as BPL167, DH 668, and DH 908. On the other hand, we noticed that all these less active stars have no rotation periods available till the date. In other words, stars having detected rotation periods show larger activity in the cooler end, which makes the issue complex. Therefore, their memberships and potential light variation through monitoring indeed need to be further investigated. Figure 6. View largeDownload slide RV distributions of Pleiades candidate members with Teff < 4000 K in inactive regime, shown by the grey line. The red line is a Gaussian fitting, which gives μ ≈ 6.2 km s−1 and σ ≈ 4.6 km s−1. The expected Pleiades RV at 5.9 km s−1 provided by Mermilliod, Mayor & Udry (2009) is marked as a black dashed line. Figure 6. View largeDownload slide RV distributions of Pleiades candidate members with Teff < 4000 K in inactive regime, shown by the grey line. The red line is a Gaussian fitting, which gives μ ≈ 6.2 km s−1 and σ ≈ 4.6 km s−1. The expected Pleiades RV at 5.9 km s−1 provided by Mermilliod, Mayor & Udry (2009) is marked as a black dashed line. The large contrast of continuum flux at H α line between M-type and hotter stars dilutes their larger differences in equivalent widths, which makes the upper envelops of $$R^{\prime }_{\text{H}\,\alpha}$$ among these stars have similar values, e.g. $$\log R^{\prime }_{\text{H}\,\alpha }\sim -3.7$$ (see bottom panels of Fig. 2). Still, the upper envelop of $$R^{\prime }_{\text{H}\,\alpha}$$ in Pleiades M-type stars (e.g. Teff ≲ 3900 K) show a slight increment with decreasing effective temperature, confirms the previous finding that the low-mass stars in young clusters with age less than 150 Myr show an increase in activity strength (LH α/Lbol) towards lower mass (Hawley et al. 1999). Also, previous studies show that later spectral type stars have lower activity level among late M-type stars (e.g. Hodgkin et al. 1995; West et al. 2004). Though such trends are evident from $$R^{\prime }_{\text{H}\,\alpha}$$ when Teff ≲3300 K, and $$R^{\prime }_{\text{H}\,\beta}$$ (see Fig. A1), we cannot make a solid conclusion about this issue based on current data. Another remarkable feature in Fig. 2 is a large scatter of $$R^{\prime }_{\text{H}\alpha}$$ among hotter Pleiades members. In particular, there evidently are two sequences among GK-type Pleiades candidates, i.e. one sequence of stars have similar activity strengths, and other sequence being with lower activity strength showing an increment towards lower temperature. The largest difference of activity level between these two sequences, e.g. 0.5–0.6 dex in $$\log R^{\prime }_{\text{H}\alpha}$$ (a factor of about four in $$R^{\prime }_{\text{H}\alpha}$$), occurs at about 5000 K. Such a difference in activity levels is unlikely due to measurement uncertainties, indicating bi-sequence behaviour among these stars. Therefore, apart from the stellar mass, other factors may play a role in producing the observed dispersion pattern. In fact, we have noticed these two sequences correspond to rotational sequences (see Section 4.2.1). It is known that the CA depends not only on stellar mass but also on age, e.g. the H α emission lasts longer in cooler stars (Stauffer et al. 1991; Hawley et al. 1999). Indeed, Fig. 2 shows a clear age effect on CA, i.e. compared to Pleiades, Hyades and Praesepe are typically less active, and M34 with intermediate age between Pleiades and Praesepe/Hyades shows activity levels in-between them. Additionally, the onset of H α direct emission occurs at about Teff ∼5000–5500, 4500–5000, 4000–4500 K for Pleiades, M34, and Hyades/Praesepe members, respectively. Furthermore, two activity sequences appear among cooler members (e.g. late K- and M-type stars) in Praesepe and Hyades, similar to the feature observed for GK-type members in Pleiades. Interestingly, the members of Praesepe and Hyades with Teff less than 3600 K (spectral type later than M1–M2) have activity levels similar to those of Pleiades active members; it indicates that the activity lifetime for M1–M2-type stars is around 600–700 Myr, a value agreed with that reported by West et al. (2008). It is possible to get such kind of activity-time-scale for earlier spectral type stars by comparing M34 with Pleiades, however, we could not able to provide any solid conclusion due to limited sample of M34 candidates in our data. 4.2 Rotation–activity connection 4.2.1 Rotation and activity sequences As discussed above, there exist two sequences among GK-type stars in Pleiades, and M-type stars in Praesepe and Hyades in $$R^{\prime }_{\text{H}\alpha}$$–Teff diagram, i.e. H α emissions in Pleiades begin to bifurcate at ∼6000 K and merge again at ∼4000 K; similar sequences in Praesepe and Hyades appear at ∼4000 K and disappear at about 3100 K. In order to understand further on this scenario, we provide $$R^{\prime }_{\text{H}\alpha}$$ as a function of temperature again in Fig. 7, wherein shown only members with known rotation periods. For comparison purpose, we also plot the corresponding rotation period (Prot) distribution in the upper panels, where the Prot values were collected from the literature (see Appendix E). The colour contrast in Fig. 7 denotes the Rossby Number Ro (Ro = Prot/τ, normalized rotation period by the convective turnover time τ). The τ was estimated from stellar mass (mass were converted from the quiescent photosphere temperature using PARSEC models, see Appendix D) using the correlation between convective turnover time and mass, derived by Wright et al. (2011). It is clear that stars in the upper sequence rotate faster than those in the lower activity sequence. Figure 7. View largeDownload slide Rotation period (top), $$R^{\prime }_{\text{H},\alpha}$$ (middle), and spot filling factor (bottom) as a function of Teff for stars with rotation period available. Left: Pleiades. Right: M34, Praesepe, and Hyades. The colour gradient indicates Rossby numbers. Overplotted in top-right panel are contours of constant Rossby numbers (Ro = 0.01, 0.1, and 1.0; grey dashed lines) based on the empirically calibrated mass–convective overturn time relation by Wright et al. (2011). The canonical rotational I/C sequences are illustrated by black lines in top panels. Figure 7. View largeDownload slide Rotation period (top), $$R^{\prime }_{\text{H},\alpha}$$ (middle), and spot filling factor (bottom) as a function of Teff for stars with rotation period available. Left: Pleiades. Right: M34, Praesepe, and Hyades. The colour gradient indicates Rossby numbers. Overplotted in top-right panel are contours of constant Rossby numbers (Ro = 0.01, 0.1, and 1.0; grey dashed lines) based on the empirically calibrated mass–convective overturn time relation by Wright et al. (2011). The canonical rotational I/C sequences are illustrated by black lines in top panels. We noticed that there exist similar features in cool star-spot coverages (spot filling factors), as shown in the lower panels of Fig. 7, where the spot filling factors were estimated by modelling their TiO band near 7050 Å (TiO2) using the standard values obtained based on references stars with solar metallicity (see Paper I for more details). It is clear that there exists a trend for fast rotators among K-type stars in Pleiades having larger spot coverages, and the spot coverages belong to two branches among GK-type Pleiades members. Further, we collected X-ray emission data from literature (Wright et al. 2011; Gondoin 2012) for stars in these open clusters, and found the existence of similar feature with dual sequences, as shown in Fig. 8. Figure 8. View largeDownload slide Fractional X-ray emission luminosities RX as a function of Teff for stars in Pleiades (left), and M34, Praesepe, and Hyades (right). These data were collected from literature, as labelled in the panels, where W11: Wright et al. (2011), G12: Gondoin (2012). Figure 8. View largeDownload slide Fractional X-ray emission luminosities RX as a function of Teff for stars in Pleiades (left), and M34, Praesepe, and Hyades (right). These data were collected from literature, as labelled in the panels, where W11: Wright et al. (2011), G12: Gondoin (2012). It is known that the open cluster members locate in three well-known regions in the rotation–colour diagram (Barnes 2003a), i.e. I (Interface) and C (Convective) sequence stars and gap stars between I and C sequence (e.g. see the top-left panel of Fig. 7). From Fig. 7, we can get a general conclusion that the activity level sequences correspond to the rotation sequences, namely, the more active sequence stars correspond to the C sequence stars, and the stars in less active sequence map to I sequence stars. The dynamo configuration in I sequences may be different with that in C sequences, e.g. an interface dynamo dominates in I sequence stars, and a turbulent or convective dynamo works in C sequence stars, correspondingly causing a large-scale interface magnetic field and a small-scale turbulent field (Barnes 2003a,b). Therefore, such a bifurcation of CA of Pleiades GK-type members is probably due to a change in the dynamo configuration. For older open clusters such as Hyades with an age of 600–700 Myr their GK-type member stars rotate slowly and should be in I sequence in the period–colour diagram, but the M members display clearly both I and C sequence (see the top-right panel of Fig. 7). As mentioned above, indeed there is a presence of two H α emission sequences for early and Mid-M stars in Praesepe and Hyades. Moreover, it is clear that the two activity sequences correspond to I- and C-rotation sequences, as shown in Fig. 7. For M34, an open cluster with an age of about 220 Myr, the two activity sequences should be observed among K-type stars. However, we could not notice such feature due to lack of observational data of the M34 K-type stars. 4.2.2 Rotation–activity correlation We have investigated the rotation–activity connection in terms of Rossby number Ro, however, it is difficult to judge whether Rossby number is superior to rotation period when modernizing this rotation–activity relation (e.g. see a comprehensive study on this topic by Reiners et al. 2014). In Fig. 9 the $$R^{\prime }_{\text{H}\,\alpha}$$ against Ro is displayed for our sample stars whose periods are available. As expected, it shows two different regimes, saturated and unsaturated, with different dependence of activity on rotation around a turnoff point (e.g. Ro ∼ 0.1), usually shown by X-ray emissions (e.g. Wright et al. 2011), chromospheric Ca ii and H α emissions (e.g. Noyes et al. 1984; Jackson & Jeffries 2010; Douglas et al. 2014; Newton et al. 2017), and cool spot coverages (Fang et al. 2016). Figure 9. View largeDownload slide $$R^{\prime }_{\text{H}\,\alpha}$$, $$R^{\prime }_{\text{H}\,\beta}$$, and $$R^{\prime }_{\text{Ca K}}$$ versus Ro for stars in Pleiades (left) and M 34, Praesepe, and Hyades (right). The colour gradient represents different effective temperatures. The cross symbols show slowly rotating Kepler M dwarfs (Prot > 10  d, Mass <0.6 M⊙). The fully convective candidates are marked by large circles. The grey solid and dashed line in top-right panel denotes the relation of Newton et al. (2017) and Douglas et al. (2014), respectively. Figure 9. View largeDownload slide $$R^{\prime }_{\text{H}\,\alpha}$$, $$R^{\prime }_{\text{H}\,\beta}$$, and $$R^{\prime }_{\text{Ca K}}$$ versus Ro for stars in Pleiades (left) and M 34, Praesepe, and Hyades (right). The colour gradient represents different effective temperatures. The cross symbols show slowly rotating Kepler M dwarfs (Prot > 10  d, Mass <0.6 M⊙). The fully convective candidates are marked by large circles. The grey solid and dashed line in top-right panel denotes the relation of Newton et al. (2017) and Douglas et al. (2014), respectively. The rotation–activity relationship is often parametrized by a flat region connected to a power-law decay, i.e. the activity has constant (saturated) value below a turnoff of Ro (namely Rosat), and above Rosat it declines as a power law of Ro, indicates activity is proportional to Roβ. Wright et al. (2011) found that the decline in coronal emission follow a power law with β = −2.18 for Ro > 0.13 based on large sample of 824 late-type stars. Furthermore, they derived β = −2.7 for a set of solar-type stars. However, the chromospheric emissions show a much shallower power-law decay for slow rotators, e.g. the first two Ca ii IRT excess emissions in early M-type dwarfs in NGC 2516 decrease as a power law of β ∼ −1 for Ro ≳ 0.1 (Jackson & Jeffries 2010). Recently, Douglas et al. (2014) derived much shallower power law with β = −0.73 for Ro > 0.11 in H α emissions of late K and M dwarfs of Praesepe and Hyades. More recently, Newton et al. (2017) detected a power law index β = −1.7 with Rosat = 0.21 in excess H α emission based on large sample of nearby M dwarfs. The rapid rotators with Ro ≲ 0.1 in Fig. 9 indicates that the chromospheric emissions have similar values shaping a plateau feature (e.g. Delfosse et al. 1998; Reiners, Joshi & Goldman 2012; Douglas, Agüeros & Covey 2014; Newton et al. 2017). Approximately, H α emission saturates at $$\log R^{\prime }_{\text{H}\,\alpha }\approx -3.7$$ dex, while H β and Ca ii K emission saturate respectively at $$\log R^{\prime }_{\text{H}\,\beta }\approx -4.1$$ dex and $$\log R^{\prime }_{\text{Ca K}}\approx -4.0$$ dex. The activity level is usually thought to be a constant value in the saturated regime. However, our data show that the chromospheric emissions still increase slightly as Ro decreases, a phenomenon exists in other activity indicators, e.g. Ca ii IR emissions of early M-type members of NGC 2516 (Jackson & Jeffries 2010). We found that it can be roughly fitted by some power of Ro, e.g. $$R^{\prime }_{\text{H}\,\alpha }\propto \text{Ro}^{-0.2}$$ for Pleiades stars with Ro ≲ 0.1, as shown by black dashed line in the top-left panel of Fig. 9, a relationship consistent with that (LX/Lbol ∝ Ro−0.16) derived by Reiners et al. (2014) based on X-ray emission. A similar increment could be seen in the H β emission, e.g. $$R^{\prime }_{\text{H}\,\beta }\propto \text{Ro}^{-0.28}$$ (the black dashed line in middle-left panel of Fig. 9). We could not see such trend in Ca ii K emission in this regime for total sample stars, which may be due to large uncertainties in $$R^{\prime }_{\text{Ca K}}$$ for cooler stars. However, one can see that there also exist a slight increase relation in GK-type stars, e.g. $$R^{\prime }_{\text{Ca K}}\propto \text{Ro}^{-0.2}$$ (the dashed line in bottom-left panel of Fig. 9). Such a slight slope of the activity–rotation relationship in the saturated regime indicates some unclear dependence of the activity on rotation is remaining even when saturation is reached, as pointed out by Reiners et al. (2014). The activity levels of ultrafast rotating G/K stars are decreasing below the saturation level, a phenomenon known as supersaturation which has been detected in the coronal emission (e.g. Randich et al. 1996; Wright et al. 2011). we found no clear evidence in our sample for any systematic fall in CA levels when the Rossby numbers become very small. However, we noticed that the average CA of M-type stars follows a decreasing trend as rotation rate increases at Ro ≲ 0.01 (particularly in terms of the excess Ca ii K emission), showing a weak evidence of damped cooling efficiency in ultrafast rotators. The regime with Ro ≳ 0.1 is usually thought to be unsaturated regime. At first sight, the activity decreases with increasing Ro, which is expected from qualitative arguments based on the αΩ-type shell dynamo theory (e.g. Noyes et al. 1984; Charbonneau 2014, and references therein). Moreover, the stars with Ro in the regime of 0.1–0.4 have large scatter in activity with upper envelop near to saturated activity level, which makes it difficult to explain the location of turn-off point. Even worse is the strong dependence of power law with onset of saturation, Rosat, i.e. β anticorrelates with Rosat, as seen in Douglas et al. (2014) and Newton et al. (2017). Anyway, we sought to represent the potential trends using the canonical rotation–activity relation (a saturated regime and a power-law decay) following the previous studies, being fitted by eye rather than using any χ2-methods. For Pleiades candidate members (see left-hand panels), we found that the saturation may occurs at Ro near to 0.2 rather than the canonical value of 0.1, and a power law is likely in the range of −2.0 to − 2.5, e.g. it can be represented by a relation of $$R^{\prime }_{\text{H}\,\alpha }\propto \text{Ro}^{-2.2}$$ for Ro ≳ 0.19, as illustrated by the black solid line in top-left panel of Fig. 9. We found a similar power-law decay with Rosat near to 0.2 in H β emissions for the Pleiades sample stars, as illustrated by the black solid line ($$R^{\prime }_{\text{H}\,\beta }\propto \text{Ro}^{-2.3}$$ for Ro ≳ 0.19) in middle-left panel of Fig. 9. Compared to these two Balmer emissions, the Ca ii K emissions show a slight different trend, i.e. the unsaturated regime start at a larger Ro (e.g. Rosat ≳ 0.25), but appears to be a reasonably fit by a similar power law (an example relation with β = −2.3 and Rosat = 0.27 is shown by the black solid line in bottom-left panel). The results for M34, Hyades, and Praesepe are shown in the right-hand panels of Fig. 9. For comparison purpose, in each right-hand panel, we re-plotted the black lines in corresponding left-hand panel. Among these older stars, most hotter stars are slow rotators, in particular members of Praesepe and Hyades, which are in unsaturated regime as expected, and only some M dwarfs in the Praesepe and Hyades and K-type members in M34 are fast rotators (see top panels of Fig. 7) that are occupied in the saturated regime. In the saturated regime, their excess chromospheric emissions increase as Ro decreases similar to Pleiades stars, but systematically having lower level. In the unsaturated regime, they appear to still follow the power law in a similar way as the Pleiades members, but with larger dispersion. However, it seems that Ca ii K emission saturates at a larger Rosat value, e.g. Rosat ≳ 0.3 (alternatively, these stars follow a shallower power laws at the same onset saturation in Pleiades). In top-right panel of Fig. 9, also shown is the relation of β = −0.73 with Rosat = 0.11 in dashed grey line, which was derived by Douglas et al. (2014) for M dwarfs with H α in emission in the Praesepe and Hyades. However, our data are inconsistent with their results, our H α data can be fitted by a clearly deeper power law with a larger Rosat value. The strong degeneracy between β and Rosat may be a potential contributor. Moreover, Douglas et al. (2014) sample consists of stars with r″ − K ≳ 3 (roughly later than K7) and definite H α emission, neglected potential less active members with H α in absorption (but in fact having excess emissions such as these GK-type members), which may also accounted for the notable difference between their results and ours. The relation for β = −1.7 with Rosat = 0.21 derived by Newton et al. (2017) for nearby M dwarfs also shown in solid grey line in the top-right panel. The Rosat shown by our H α data is consistent with their value, however, the power-law index is slightly deeper than their value. The transformation between the solar-type stars (with radiative cores and convective envelopes) and fully convective stars takes place roughly in stars of spectral type M3–M4. The interface layer between radiative cores and convective envelopes (known as the tachocline) is believed to play an important role in the generation of the magnetic field in solar-type stars (Spiegel & Zahn 1992). Therefore, the dynamo mechanism in the fully convective stars such as middle/late M-type stars is expected to be very different because of the absence of a tachocline. However, it is evident that slowly rotating, fully convective stars still show rotation–activity (X-ray emission) connection in the same way similar to solar-type stars (Wright & Drake 2016). Recently, Newton et al. (2017) found a single relation between excess H α emission and Ro for stars including both early- and late-type M dwarfs, supporting the previous finding. Unfortunately, most M-type stars in our sample are fast rotators (see Fig. 7), and only two slow rotators in Hyades are fully convective candidate stars (marked as large open circles, and see Table 2). To take a simple look at this issue, we picked up LAMOST spectra for 47 slowly rotating M dwarfs (Prot >10 d, mass < 0.6 M⊙) from the main-sequence Kepler targets with known rotation periods provided by McQuillan, Mazeh & Aigrain (2014), and measured their excess CA emissions, as shown by crosses in right panels of Fig. 9 (only 29 of them show excess H α emission). Note that two of them are the fully convective candidates, which were marked with large open circles (but also see Table 2). It is clear that these slowly rotating M dwarfs, including the fully convective candidates, still follow the power-law decay that show no evident difference from hotter stars, which is consistent with the finding of Newton et al. (2017), though our results show a slight different power-law decay. If this scenario is real, it then suggests that the magnetic dynamo may be similar to that in solar-type stars. Wright & Drake (2016) suggest a common magnetic dynamo in both partially and fully convective stars, supporting those models in which the dynamo originates throughout the stellar convection zone. Table 2. Slow rotating, potential fully convective stars with detected excess emissions in both H α and H β lines. Object name  SpT.  Teff  Mass  Prot  Ro  $$\log R^{\prime }_{\text{H}\,\alpha}$$  $$\log R^{\prime }_{\text{H}\,\beta}$$  $$\log R^{\prime }_{\text{Ca K}}$$  Note      (K)  (M⊙)  (d)            04202761+1853499  M3.6  3244  0.375  20.3105  0.4079  −4.3609  −4.7871  −4.5722  Hyades  04303385+1444532  M3.5  3257  0.382  18.4100  0.3771  −4.4220  −4.8098  −4.6805  Hyades  KIC 10731839  M3.3  3272  0.362  46.1150  0.8946  −5.4351  −5.7745  −5.2008  Kepler target  KIC 11495571  M3.4  3292  0.374  13.6400  0.2737  −4.0450  −4.4538  –  Kepler target  Object name  SpT.  Teff  Mass  Prot  Ro  $$\log R^{\prime }_{\text{H}\,\alpha}$$  $$\log R^{\prime }_{\text{H}\,\beta}$$  $$\log R^{\prime }_{\text{Ca K}}$$  Note      (K)  (M⊙)  (d)            04202761+1853499  M3.6  3244  0.375  20.3105  0.4079  −4.3609  −4.7871  −4.5722  Hyades  04303385+1444532  M3.5  3257  0.382  18.4100  0.3771  −4.4220  −4.8098  −4.6805  Hyades  KIC 10731839  M3.3  3272  0.362  46.1150  0.8946  −5.4351  −5.7745  −5.2008  Kepler target  KIC 11495571  M3.4  3292  0.374  13.6400  0.2737  −4.0450  −4.4538  –  Kepler target  View Large 4.3 Correlations between activity indicators The correlations among activity indicators of different temperature layers contain information of overall heating process and systematic deficiencies in these layers (Güdel 2004). Therefore, the correlations between different CA indicators and the relation of excess chromospheric emission with other activity indicators had gained lot of attention (e.g. Schrijver, Dobson & Radick 1992; Montes et al. 1995, 1996; Martínez-Arnáiz et al. 2011). In fact, this topic need further investigation, e.g. need complex physical models to interpret the correlations, which is beyond the scope of this paper. Our intent here is to merely investigate the basic observational correlations between different tracers of activity. 4.3.1 H α, H β, and Ca ii K Ca ii H&K lines are the most widely used CA indicators as their source function are collisionally dominated. H α and H β are the two strongest lines among Balmer series, and their excess emissions are formed at the middle chromosphere (Montes et al. 2004). Our data show that there exist good correlations of excess emissions of H α to that of H β and Ca ii K lines, as shown in Figs 10 and 11. Figure 10. View largeDownload slide The correlation of excess equivalent widths (upper panels) and fractional excess luminosities (lower panels) of H β and H α lines. Dashed lines in upper panels: power-law approximation (see the text). The dashed lines with grey colours in the lower panels denote constant ratios of $$R^{\prime }_{\text{H}\,\alpha }/R^{\prime }_{\text{H}\,\beta } = 0.5, 1, 2, 4,\ \text{and}\ 8$$, respectively. Figure 10. View largeDownload slide The correlation of excess equivalent widths (upper panels) and fractional excess luminosities (lower panels) of H β and H α lines. Dashed lines in upper panels: power-law approximation (see the text). The dashed lines with grey colours in the lower panels denote constant ratios of $$R^{\prime }_{\text{H}\,\alpha }/R^{\prime }_{\text{H}\,\beta } = 0.5, 1, 2, 4,\ \text{and}\ 8$$, respectively. Figure 11. View largeDownload slide Relationships between excess emissions of H α and Ca ii K lines (see the text). Figure 11. View largeDownload slide Relationships between excess emissions of H α and Ca ii K lines (see the text). Fig. 10 shows a power-law correlation (linear relation in log–log scale) between excess equivalent widths of the two Balmer lines, EW$$^{\prime }_{\text{H}\,\alpha}$$ and EW$$^{\prime }_{\text{H}\,\beta}$$, among the Pleiades stars. We found that it can be represented well by a power-law index in the range of 1.2–1.4, e.g. EW$$^{\prime }_{\text{H}\,\beta } \propto ({\text{EW}}^{\prime }_{\text{H}\,\alpha })^{1.3}$$, as illustrated by the black dashed line in top-left panel of Fig. 10. For members in M34, Praesepe, and Hyades, it also can be represented roughly by the same power-law relation, as shown by the black dashed line in top-right panel. It is noticeable that the main trend between $$R^{\prime }_{\text{H}\,\beta}$$ and $$R^{\prime }_{\text{H}\,\alpha}$$ also roughly linear in log–log scale. Moreover, we found that the power-law index is around 1.0, which indicates a linear relation roughly even in linear scale, e.g. $$R^{\prime }_{\text{H}\,\beta } \propto R^{\prime }_{\text{H}\,\alpha}$$, as demonstrated by dashed lines with constant ratios ($$R^{\prime }_{\text{H}\,\alpha }/R^{\prime }_{\text{H}\,\beta } = 0.5, 1, 2, 4, {\text{ and }} 8$$, from top to bottom, respectively) in lower panels of Fig. 10. A constant H α to H β ratio has been detected in previous studies, e.g. Reid et al. (1995) found a good linear correlation between the flux emitted at H α and H β for low-mass stars in the Hyades (FH α ∼ 4.6 × FH β). However, some late F-type or early G-type stars show departures from the general trend, e.g. their H β losses seem to be more important than H α losses. In fact, compared to H α, H β is preferred the higher chromospheric density, e.g. while H β enhanced during a flare-like event, the H α is not proportionally increased in strength (Huenemoerder & Ramsey 1987). Thus the ratio of energy emitted in the H α to H β ($${{E}}_{{\rm H}\,\,\alpha }/{{E}}_{\rm{H}\,\,\beta}$$) is widely used as indicator of the presence of flare-like events (Huenemoerder & Ramsey 1987), and also a diagnostic indicator for discriminating between stellar plages and prominences (Hall & Ramsey 1992). The typical value of $${{E}}_{\rm{H}\,\alpha }/{{E}}_{\rm{H}\,\beta}$$ is around 1 in flare-like events, it then become 4 or larger for RS CVn stars (Huenemoerder & Ramsey 1987). Therefore, the departure that some earlier type stars with lower H α to H β ratios ($$R^{\prime }_{\text{H}\,\alpha }/R^{\prime }_{\text{H}\,\beta }<1$$) were likely underlying flare/plage-like activities. If this is the real case, then it suggests that these stars preferred to suffer flare events than cooler stars. It is clear from Fig. 11 that there exist power-law relationship (linear relation in log–log scale) between excess equivalent widths of H α and Ca ii K lines. Here we did not attempt any χ2 fitting but merely sought to reproduce the correlation by eye as before. It can be represented approximately by a power-law index of 1.0–1.1, e.g. EW$$^{\prime }_{\text{Ca K}} \propto ({\text{EW}}^{\prime }_{\text{H}\,\alpha })^{1.05}$$ (see the black dashed lines in upper panels), for which the power-law exponent is slight larger than the value (0.97) derived by Montes et al. (1995) based on a sample of 51 chromospherically active binaries. The lower panels of Fig. 11 correspond to relations between excess fractional luminosities of Ca ii K and H α, which could be represented approximately by a relation of $$R^{\prime }_{\text{Ca K}}\propto ({R^{\prime }_{\text{H}\,\alpha }})^{0.5}$$ for GK-type Pleiades stars (see the black dashed lines in lower panels). However, cooler and more chromospherically active stars such as M-type stars do not obey the general relations, showing departure feature of much sharper slope compared to hotter stars. It is illustrated as an example relation $$R^{\prime }_{\text{Ca K}}\propto ({R^{\prime }_{\text{H}\,\alpha }})^{2.0}$$ by grey solid lines in the figure, which indicates that H α losses tend to be more important than Ca ii K losses for cooler stars that are more active. In fact, previous studies show such departures from the general trend in surface flux–flux relations between H α and other chromospheric indicators among some late-K and M dwarfs, and they are found to be young stars or flare stars being in the X-ray saturation regime (e.g. Martínez-Arnáiz et al. 2011). One can see that the departures hosted in M stars with non-universality trends appear both in young open cluster, Pleiades, and median-age open clusters, Hyades, and most of these stars are indeed in H α emission saturation regime. It has been suggested that the stars deviating from the general relation probably have a different magnetic structure (Martínez-Arnáiz et al. 2011), e.g. nanoflare heating. 4.3.2 H α and photospheric activity In Fig. 12, we showed the spot filling factor (or fractional spot coverage, fs) of stars with detected cool spots (fs > 0) in Pleiades as a function of $$R^{\prime }_{\text{H}\,\alpha}$$. The spot coverage of stars have large scatter, which makes harder to say about existence of any believable connections between these two activity indicators. However, the overall trend, in particular among more chromospherically active members (e.g. if we consider stars with $$\log R^{\prime }_{\text{H}\,\alpha }>-5.0$$ and keep aside the outliers at the end of weak H α emission), shows that more active stars tend to have larger spot coverages. We use a Pearson's product–moment correlation analysis to test the statistical significance of the potential connections between fs and $$\log R^{\prime }_{\text{H}\,\alpha}$$. The Pearson's correlation coefficient r ≈ 0.55 among stars with $$\log R^{\prime }_{\text{H}\alpha }>-5.0$$ indicates that there exist some correlation between spot coverage and chromospheric emission. Compared to hotter stars such as GK-type members, most M-type Pleiades members (see red symbols in Fig. 12) are in H α emission saturation regime, showing a different and slightly tighter trend in the plot, which also evidenced by a stronger positive correlation of r ≈ 0.63 (for stars with Teff < 3900 K and $$\log R^{\prime }_{\text{H}\,\alpha }>-4.0$$). It is clear that such non-universality feature is very similar to the relation between chromospheric emissions of H α and Ca ii K, as discussed above. Note that the uncertainties are larger in determination of spot coverages for M-type star as discussed in Paper I, and the correlation showed in the plot needs to be further investigated, which puts a caution for further discussion based on these features. Figure 12. View largeDownload slide Relations between $$R^{\prime }_{\text{H}\,\alpha}$$ and spot filling factor for the Pleiades sample stars with detected cool spots coverage. Colours as in Fig. 9. Figure 12. View largeDownload slide Relations between $$R^{\prime }_{\text{H}\,\alpha}$$ and spot filling factor for the Pleiades sample stars with detected cool spots coverage. Colours as in Fig. 9. The amplitude of light variation due to spot rotational modulation carries information of inhomogeneity in the stellar surface, thus acts as another indicator of photospheric activity, e.g. the light variation becomes larger with increase in CA among Sun-like stars (e.g. Lockwodd, Skiff & Radick 1997; Radick et al. 1998; Lockwodd et al. 2007) and M dwarfs (e.g. Newton et al. 2017). The amplitudes of periodic light variation versus excess H α emissions for Pleiades candidates with detected rotation period are shown in Fig. 13, where the r-band peak-to-peak photometric amplitudes Ar were derived by Hartman et al. (2010) based on the Hungarian-made Automated Telescope Network (HATNet) light curves. On the whole, there is a trend showing larger chromospheric emission corresponds to larger light variation, however, the dispersion is large as shown in the left-hand panel. The Pearson's correlation coefficient in log–log scale is r ∼ 0.71 (with a p-value of < 2.2 × 10−16). The right-hand panel of Fig. 13 is analogous to the left-hand panel, wherein the amplitude of variation (between 10 per cent and 90 per cent of the distribution of points) AK2 is derived by Rebull et al. (2016) based on K2 observations. We noticed some additional features from this panel. The correlation is weaker that is evidenced by a smaller value of r ∼ 0.47 (p ∼ 1.4 × 10−9). And there exist systematic offset between the AK2 and Ar for very cool stars (see red symbols). By neglecting the very cool stars with Teff < 3900 K, the correlation become tighter (r ∼ 0.62, p ∼ 1.9 × 10−10). Previous studies show that the short-term photometric variability (mainly due to rotational spot modulation) in less active Sun-like stars can be related to their average level of CA by a power law (e.g. Radick et al. 1998). The Pleiades GK-type members in the unsaturated regime are roughly linear on the log–log scale plot, indicating a power-law relation. On the other hand, M dwarfs still follow the general trend of more active stars showing higher levels of photometric variability (e.g. Newton et al. 2017). As most of the Pleiades M-type stars are in saturated CA regime, there is no believable trend among these very cool stars with Teff < 3900 K in both panels (left: r ∼ 0.21, p ∼ 0.12; right: r ∼ 0.08, p ∼ 0.50), which probably due to their very high activity levels. In fact, the CA levels vary in short term (as shown in Fig. 5), and long term such as cyclic variation (e.g. Baliunas et al. 1995), thus the panels of Fig. 13 actually represent a snapshot of stellar variability. Moreover, the pattern and evolution of spots, and stellar inclination (the orientation angle of the axis of rotation relative to our line of sight) could effect the light variation pattern and contaminate the potential relation between chromospheric emission and light variation, thus the feature shown in this figure is only instructive rather than conclusive. Figure 13. View largeDownload slide Amplitudes of light variation (left: HATNet, right: K2) as a function of $$\log R^{\prime }_{\text{H}\,\alpha}$$. Note that the amplitude is in log scale. Figure 13. View largeDownload slide Amplitudes of light variation (left: HATNet, right: K2) as a function of $$\log R^{\prime }_{\text{H}\,\alpha}$$. Note that the amplitude is in log scale. 4.3.3 H α and coronal X-ray emission To investigate the relation between chromospheric heating and coronal emission, we have taken RX (=LX/Lbol) from Wright et al. (2011) for common Pleiades candidate members, and plotted against $$R^{\prime }_{\text{H}\,\alpha}$$ in Fig. 14. A simple linear fit in log scale indicated that the power law index is around 1.44, e.g. $$R^{\prime }_{\text{H}\,\alpha }\propto (R^{\prime }_{\text{H}\,\alpha })^{1.44}$$ (see black dashed line in Fig. 14). A slight steeper correlation of coronal and chromospheric emission flux with a power of 1.5–1.7 based on magnetic induced excess flux in the Ca ii H&K lines was reported (Schrijver et al. 1992; Güdel 2004, and further references cited therein). Such a non-linearity in the chromospheric-coronal flux–flux correlation may be a manifestation of the chromospheric radiative losses (Güdel 2004). The correlation between H α emission and X-ray emission indicates that the formation process may be same or at least correlated. Figure 14. View largeDownload slide The relations between H α fractional luminosity $$R^{\prime }_{\text{H}\,\alpha}$$ and X-ray fractional luminosity RX. Figure 14. View largeDownload slide The relations between H α fractional luminosity $$R^{\prime }_{\text{H}\,\alpha}$$ and X-ray fractional luminosity RX. 4.4 What can we learn from the results The different onsets of saturation in various CA indicators suggest the existence of transition regime at Ro range of 0.1–0.4, which connects the unsaturated and saturated regime with steep and shallower power law, respectively. In fact, Gondoin (2012, 2013) showed the evidence for transition from saturated to non-saturated X-ray emission that occurs at Ro between ∼0.13 and 0.4 in open clusters M34 and M35. Indeed, this transition scenario could describe the scatter in Ro region of 0.1–0.4. However, such scenario fails to explain the whole observational features in this Ro regime, e.g. the Ca K emission that shows a clear rotation–activity relationship, suggesting the investigation towards other unknown parameters or the activity behaviours in this regime. In addition, the chromospheric emission follows the similar (slight shallower) power laws in the unsaturated regime compared to coronal X-ray emission, but appears to be shifted to a larger Rosat value, which indicates that chromospheric emission gets easily saturated at a given rotation rate than coronal emission. Alternatively, it may suggest that the convective overturn time-scale, empirically derived based on X-ray emission (Wright et al. 2011), is not suitable for chromospheric emission. A corresponding relation between the saturated/non-saturated regimes of X-ray emission and the C/I rotation sequences among late-type stars has been reported in literature (e.g. Barnes 2003b; Gondoin 2012, 2013). In our work, we detected a similar correspondence in H α emission among cool stars in Pleiades, Praesepe, and Hyades, as discussed in previous sections. Barnes (2003b) and Gondoin (2012) suggested such transition from the saturated to the non-saturated regime of activity is the manifestation of a dynamo transition, e.g. from a turbulent dynamo to an interface-type dynamo. Durney, De Young & Roxburgh (1993) showed a candidate turbulent dynamo, which can generate small-scale turbulent magnetic field in late-type stars, and have a weak dependence on rotation. Such a turbulent dynamo can naturally explain the weak dependence of activity on rotation among fast rotators in the saturation regime. Moreover, the generation of small-scale magnetic field supports some observational behaviour shown by very active stars, e.g. large spot coverage of stars in young open clusters may be due to many small and randomly located spots on the stellar photosphere (Jackson & Jeffries 2012, 2013; Fang et al. 2016), and the possibility of multiple small flares such as nanoflares that might be responsible for the departures of relation between various chromospheric emissions among very active late-K and M dwarfs reported by previous studies and discussed above in current work. We thus conclude that a turbulent-like dynamo dominates in fast rotating cool stars in the saturated regime. 5 CONCLUSION We have quantified CA (excess emissions in H α, H β, and Ca ii K lines) in more than 700 stars from various open clusters including Pleiades, M34, Praesepe, and Hyades, using LAMOST DR3 spectra. The CA was investigated as a function of stellar parameters: mass, age, and rotation. We also analysed the correlation between activity indicators of various atmosphere layers. CA was found to be overall increased as temperature decreases, showing a mass-dependence, which partially could be due to deeper convective zones in cooler stars. CA also shows a clear age-dependence, i.e. activity level of GK-type stars declines steadily from 100 Myr to about 700 Myr. Additionally, we found that the activity lifetime for M1–M2-type stars is around 600–700 Myr, which is consistent with previous studies. More interestingly, there appear two sequences among Pleiades GK-type stars in activity (spot coverage, chromospheric, and coronal emission)–Teff diagram, paralleling with well-known rotation sequences. We also found the existence of similar sequences in Praesepe and Hyades M-type members. We found CA is correlated with rotation in the unsaturated regime, following a power law with power of β between −2 and −2.5, which is marginally consistent with results from X-ray emissions, but with a larger turnoff point of Ro (e.g. Rosat ∼ 0.2–0.3). A weak dependence of CA on rotation is noticed in the saturated regime of activity (Ro ≲ 0.1), e.g. a slight increase as Ro decreases (a power law with β ≈ 0.2), peaks around Ro ∼ 0.01, and decrease when Ro decreases for Ro ≲ 0.01. In addition, our analysis confirms the previous finding that slowly rotating, fully convective stars also follow the activity–rotation trend like the hotter stars. We found a good linear correlation between the luminosity emitted at H α and H β. The relations from excess emission of H α to that Ca ii K lines confirm that cool stars follow two different power-law flux–flux relationships. We also found weak correlation between chromospheric emission and photospheric activity indicators (cool spot filling factor and light variation amplitude), however, it shows dependence on spectral type and activity level, which provides clues on how spot configuration varies as a function of mass and activity. A good correlation between chromospheric and coronal emissions was also detected that confirms previous findings. ACKNOWLEDGEMENTS We thank the referee for his/her observant comments and constructive suggestions that helped to improve the manuscript. This study is supported by the National Natural Science Foundation of China under grant No. 11390371, 11233004, 11573035, 11625313, 11550110492, the National Key Basic Research Program of China (973 program) 2014CB845701. This work is supported by the Astronomical Big Data Joint Research Center, co-founded by the National Astronomical Observatories, Chinese Academy of Sciences and the Alibaba Cloud. YBK is thankful to the Chinese Academy of Sciences for supporting through CAS PIFI Fellowship. This work has made use of LAMOST data. The Guo Shou Jing Telescope (the Large sky Area Multi-Object fibre Spectroscopic Telescope, LAMOST) is a National Major Scientific Project built by the Chinese Academy of Sciences. Funding for the project has been provided by the National Development and Reform Commission. LAMOST is operated and managed by National Astronomical Observatories, Chinese Academy of Sciences. This research has made use of the VizieR catalogue access tool and the cross-match service provided by CDS, Strasbourg, France. This research has made use of the WEBDA data base, operated at the Department of Theoretical Physics and Astrophysics of the Masaryk University. REFERENCES Agüeros M. A. et al.  , 2011, ApJ , 740, 110 https://doi.org/10.1088/0004-637X/740/2/110 CrossRef Search ADS   An D., Terndrup D. M., Pinsonneault M. H., Paulson D. B., Hanson R. B., Stauffer J. R., 2007, ApJ , 655, 233 https://doi.org/10.1086/509653 CrossRef Search ADS   Baliunas S. L. et al.  , 1995, ApJ , 438, 269 https://doi.org/10.1086/175072 CrossRef Search ADS   Baliunas S. L., Donahue R. A., Soon W., Henry G. W., 1998, in Donahue R. A., Bookbinder J. A., eds, ASP Conf. Ser. Vol. 154, Cool Stars, Stellar Systems, and the Sun . Astron. Soc. Pac., San Francisco, p. 153 Barnes S. A., 2003a, ApJ , 586, 464 https://doi.org/10.1086/367639 CrossRef Search ADS   Barnes S. A., 2003b, ApJ , 586, L145 https://doi.org/10.1086/374681 CrossRef Search ADS   Bouy H. et al.  , 2015, A&A , 577, 148 CrossRef Search ADS   Brandt T. D., Huang C. X., 2015, ApJ , 807, 24 https://doi.org/10.1088/0004-637X/807/1/24 CrossRef Search ADS   Canterna R., Crawford D. L., Perry C. L., 1979, PASP , 91, 263 https://doi.org/10.1086/130482 CrossRef Search ADS   Carrera R., Pancino E., 2011, A&A , 535, 30 CrossRef Search ADS   Charbonneau P., 2014, ARA&A , 52, 251 CrossRef Search ADS   Chen Y., Girardi L., Bressan A., Marigo P., Barbieri M., Kong X., 2014, MNRAS , 444, 2525 https://doi.org/10.1093/mnras/stu1605 CrossRef Search ADS   Copenhagen University Institute of Astronomy Cambridge, UK Real Instituto y Observatorio de la Armada en San Fernando, 2011, yCat , 1327, 0 Covey K. R. et al.  , 2016, ApJ , 822, 81 https://doi.org/10.3847/0004-637X/822/2/81 CrossRef Search ADS   Cram L. E., Mullan D. J., 1979, ApJ , 234, 579 https://doi.org/10.1086/157532 CrossRef Search ADS   Cui X. -Q. et al.  , 2012, RA&A , 12, 1197 https://doi.org/10.1088/1674-4527/12/9/003 Delfosse X., Forveille T., Perrier C., Mayor M., 1998, A&A , 1998, 331 , 581 Delorme P., Collier Cameron A., Hebb L., Rostron J., Lister T. A., Norton A. J., Pollacco D., West R. G., 2011, MNRAS , 413, 2218 https://doi.org/10.1111/j.1365-2966.2011.18299.x CrossRef Search ADS   Douglas S. T., Agüeros M. A., Covey K. R., 2014, ApJ , 795, 161 https://doi.org/10.1088/0004-637X/795/2/161 CrossRef Search ADS   Douglas S. T., Agüeros M. A., Covey K. R., Cargile P. A., Barclay T., Cody A., Howell S. B., Kopytova T., 2016, ApJ , 822, 47 https://doi.org/10.3847/0004-637X/822/1/47 CrossRef Search ADS   Durney B. R., De Young D. S., Roxburgh I. W., 1993, Solar Phys. , 145, 207 https://doi.org/10.1007/BF00690652 CrossRef Search ADS   Fang X.-S., Zhao G., Zhao J.-K., Chen Y.-Q., Bharat Kumar Y., 2016, MNRAS , 463, 2494 (Paper I) https://doi.org/10.1093/mnras/stw1923 (Paper I) CrossRef Search ADS   Fossati L., Bagnulo S., Landstreet J., Wade G., Kochukhov O., Monier R., Weiss W., Gebran M., 2008, A&A , 483, 891 CrossRef Search ADS   Gallet F., Bouvier J., 2013, A&A , 556, A36 CrossRef Search ADS   Geller A. M., Latham D. W., Mathieu R. D., 2015, AJ , 150, 97 https://doi.org/10.1088/0004-6256/150/3/97 CrossRef Search ADS   Gondoin P., 2012, A&A , 546, A117 CrossRef Search ADS   Gondoin P., 2013, A&A , 556, A14 CrossRef Search ADS   Gray R. O., Corbally C. J., Garrison R. F., McFadden M. T., Robinson P. E., 2003, AJ , 126, 2048 https://doi.org/10.1086/378365 CrossRef Search ADS   Gray R. O., Corbally C. J., Garrison R. F., McFadden M. T., Bubar E. J., McGahee C. E., O'Donoghue A. A., Knox E. R., 2006, AJ , 132, 161 https://doi.org/10.1086/504637 CrossRef Search ADS   Güdel M., 2004, A&AR , 12, 71 Hall J. C., 2008, Living Rev. Solar Phys ., 5, 2 https://doi.org/10.12942/lrsp-2008-2 CrossRef Search ADS   Hall J. C., Ramsey L. W., 1992, AJ , 104, 1942 https://doi.org/10.1086/116370 CrossRef Search ADS   Hartman J. D., Bakos G. Á., Kovács G., Noyes R. W., 2010, MNRAS , 408, 475 https://doi.org/10.1111/j.1365-2966.2010.17147.x CrossRef Search ADS   Hartman J. D., Bakos G. Á., Noyes R. W., Sipocz B., Kovács G., Mazeh T., Shporer A., Pál A., 2011, AJ , 141, 166 https://doi.org/10.1088/0004-6256/141/5/166 CrossRef Search ADS   Hawley S. L., Reid I. N., Gizis J. E., Byrne P. B., 1999, in Butler C. J., Doyle J. G., eds, ASP Conf. Ser. Vol. 158, Solar and Stellar Activity: Similarities and Differences . Astron. Soc. Pac., San Francisco, p. 63 Henry T. J., Soderblom D. R., Donahue R. A., Baliunas S. L., 1996, AJ , 111, 439 https://doi.org/10.1086/117796 CrossRef Search ADS   Hodgkin S. T., Jameson R. F., Steele I. A., 1995, MNRAS , 274, 869 https://doi.org/10.1093/mnras/274.3.869 CrossRef Search ADS   Huenemoerder D. P., Ramsey L. W., 1987, ApJ , 319, 392 https://doi.org/10.1086/165463 CrossRef Search ADS   Husser T.-O., Wende-von Berg S., Dreizler S., Homeier D., Reiners A., Barman T., Hauschildt P. H., 2013, A&A , 553, A6 CrossRef Search ADS   Irwin J., Aigrain S., Hodgkin S., Irwin M., Bouvier J., Clarke C., Hebb L., Moraux E., 2006, MNRAS , 370, 954 https://doi.org/10.1111/j.1365-2966.2006.10521.x CrossRef Search ADS   Jackson R. J., Jeffries R. D., 2010, MNRAS , 407, 465 https://doi.org/10.1111/j.1365-2966.2010.16917.x CrossRef Search ADS   Jackson R. J., Jeffries R. D., 2012, MNRAS , 423, 2966 https://doi.org/10.1111/j.1365-2966.2012.21119.x CrossRef Search ADS   Jackson R. J., Jeffries R. D., 2013, MNRAS , 431, 1883 https://doi.org/10.1093/mnras/stt304 CrossRef Search ADS   Jacobson H. R., Pilachowski C. A., Friel E. D., 2011, AJ , 142, 59 https://doi.org/10.1088/0004-6256/142/2/59 CrossRef Search ADS   James D. J. et al.  , 2010, A&A , 515, 100 CrossRef Search ADS   Jones B. F., Prosser C. F., 1996, AJ , 111, 1193 https://doi.org/10.1086/117865 CrossRef Search ADS   Kovács G. et al.  , 2014, MNRAS , 442, 2081 https://doi.org/10.1093/mnras/stu946 CrossRef Search ADS   Lee Y. S. et al.  , 2008, AJ , 136, 2050 https://doi.org/10.1088/0004-6256/136/5/2050 CrossRef Search ADS   Liu X.-W., Zhao G., Hou J.-L., 2015, RA&A , 15, 1089 https://doi.org/10.1088/1674-4527/15/8/001 Lockwood G. W., Skiff B. A., Radick R. R., 1997, ApJ , 485, 789 https://doi.org/10.1086/304453 CrossRef Search ADS   Lockwood G. W., Skiff B. A., Henry G. W., Henry S., Radick R. R., Baliunas S. L., Donahue R. A., Soon W., 2007, ApJS , 171, 260 https://doi.org/10.1086/516752 CrossRef Search ADS   Luo A. -L. et al.  , 2015, RA&A , 15, 1095 https://doi.org/10.1088/1674-4527/15/8/002 Mamajek E. E., Hillenbrand L. A., 2008, ApJ , 687, 1264 https://doi.org/10.1086/591785 CrossRef Search ADS   Marsden S. C., Carter B. D., Donati J.-F., 2009, MNRAS , 399, 888 https://doi.org/10.1111/j.1365-2966.2009.15319.x CrossRef Search ADS   Martínez-Arnáiz R., López-Santiago J., Crespo-Chacón I., Montes D., 2011, MNRAS , 414, 2629 https://doi.org/10.1111/j.1365-2966.2011.18584.x CrossRef Search ADS   McQuillan A., Mazeh T., Aigrain S., 2014, ApJS , 211, 24 https://doi.org/10.1088/0067-0049/211/2/24 CrossRef Search ADS   Meibom S., Mathieu R. D., Stassun K. G., Liebesny P., Saar S. H., 2011, ApJ , 733, 115 https://doi.org/10.1088/0004-637X/733/2/115 CrossRef Search ADS   Melis C., Reid M. J., Mioduszewski A. J., Stauffer J. R., Bower G. C., 2014, Science , 345, 1029 https://doi.org/10.1126/science.1256101 CrossRef Search ADS PubMed  Mermilliod J.-C., Mayor M., Udry S., 2009, A&A , 498, 949 CrossRef Search ADS   Montes D., Fernández-Figueroa M. J., de Castro E., Cornide M., 1995, A&A , 294, 165 Montes D., Fernández-Figueroa M. J., Cornide M., de Castro E., 1996, A&A , 312, 221 Montes D., Crespo-Chacón I., Gálvez M. C., Fernández-Figueroa M. J., López-Santiago J., de Castro E., Cornide M., Hernán-Obispo M., 2004, Lect. Notes Essays Astrophys. , 1, 119 Newton E. R., Irwin J., Charbonneau D., Berlind P., Calkins M. L., Mink J., 2017, ApJ , 834, 85 https://doi.org/10.3847/1538-4357/834/1/85 CrossRef Search ADS   Noyes R. W., Hartmann L. W., Baliunas S. L., Duncan D. K., Vaughan A. H., 1984, ApJ , 279, 763 https://doi.org/10.1086/161945 CrossRef Search ADS   Önehag A., Korn A., Gustafsson B., Stempels E., Vandenberg D. A., 2011, A&A , 528, A85 CrossRef Search ADS   Perryman M. A. C. et al.  , 1998, A&A , 331, 81 Pizzolato N., Maggio A., Micela G., Sciortino S., Ventura P., 2003, A&A , 397, 147 CrossRef Search ADS   Radick R. R., Thompson D. T., Lockwood G. W., Duncan D. K., Baggett W. E., 1987, ApJ , 321, 459 https://doi.org/10.1086/165645 CrossRef Search ADS   Radick R. R., Lockwood G. W., Skiff B. A., Baliunas S. L., 1998, ApJS , 118, 239 https://doi.org/10.1086/313135 CrossRef Search ADS   Randich S., Schmitt J. H. M. M., Prosser C. F., Stauffer J. R., 1996, A&A , 305, 785 Rebull L. M. et al.  , 2016, AJ , 152, 113 https://doi.org/10.3847/0004-6256/152/5/113 CrossRef Search ADS   Reid N., Hawley S. L., Mateo M., 1995, MNRAS , 272, 828 Reiners A., Joshi N., Goldman B., 2012, AJ , 143, 93 https://doi.org/10.1088/0004-6256/143/4/93 CrossRef Search ADS   Reiners A., Schüssler M., Passegger V. M., 2014, ApJ , 794, 144 https://doi.org/10.1088/0004-637X/794/2/144 CrossRef Search ADS   Salaris M., Weiss A., Percival S. M., 2004, A&A , 414, 163 CrossRef Search ADS   Scholz A., Eislöffel J., 2007, MNRAS , 381, 1638 https://doi.org/10.1111/j.1365-2966.2007.12333.x CrossRef Search ADS   Scholz A., Irwin J., Bouvier J., SipH ocz B. M., Hodgkin S., Eislöffel J., 2011, MNRAS , 413, 2595 https://doi.org/10.1111/j.1365-2966.2011.18328.x CrossRef Search ADS   Schrijver C. J., Dobson A. K., Radick R. R., 1992, A&A , 258, 432 Schuler S. C., King J. R., Fischer D. A., Soderblom D. R., Jones B. F., 2003, AJ , 125, 2085 https://doi.org/10.1086/373927 CrossRef Search ADS   Skrutskie M. F. et al.  , 2006, AJ , 131, 1163 CrossRef Search ADS   Smolinski J. P. et al.  , 2011, AJ , 141, 89 https://doi.org/10.1088/0004-6256/141/3/89 CrossRef Search ADS   Soderblom D. R., Stauffer J. R., Hudon J. D., Jones B. F., 1993, ApJS , 85, 315 https://doi.org/10.1086/191767 CrossRef Search ADS   Soderblom D. R., Jones B. F., Fischer D., 2001, ApJ , 563, 334 https://doi.org/10.1086/323829 CrossRef Search ADS   Soderblom D. R., Laskar T., Valenti J. A., Stauffer J. R., Rebull L. M., 2009, AJ , 138, 1292 https://doi.org/10.1088/0004-6256/138/5/1292 CrossRef Search ADS   Spiegel E. A., Zahn J.-P., 1992, A&A , 265, 106 Stauffer J. R., Hartmann L. W., 1986, ApJS , 61, 531 https://doi.org/10.1086/191123 CrossRef Search ADS   Stauffer J. R., Hartmann L. W., 1987, ApJ , 318, 337 https://doi.org/10.1086/165371 CrossRef Search ADS   Stauffer J. R., Giampapa M. S., Herbst W., Vincent J. M., Hartmann L. W., Stern R. A., 1991, ApJ , 374, 142 https://doi.org/10.1086/170104 CrossRef Search ADS   Stauffer J. R., Balachandran S. C., Krishnamurthi A., Pinsonneault M., Terndrup D. M., Stern R. A., 1997, ApJ , 475, 604 https://doi.org/10.1086/303567 CrossRef Search ADS   Stauffer J. R., Schultz G., Kirkpatrick J. D., 1998, ApJ , 499, L199 https://doi.org/10.1086/311379 CrossRef Search ADS   Stauffer J. R. et al.  , 2007, ApJS , 172, 663 https://doi.org/10.1086/518961 CrossRef Search ADS   Taylor B. J., 2006, AJ , 132, 2453 https://doi.org/10.1086/508610 CrossRef Search ADS   Taylor B. J., Joner M. D., 2005, ApJS , 159, 100 https://doi.org/10.1086/430422 CrossRef Search ADS   van Leeuwen F., 2009, A&A , 497, 209 CrossRef Search ADS   Vaughan A. H., Preston G. W., 1980, PASP , 92, 385 https://doi.org/10.1086/130683 CrossRef Search ADS   Walkowicz L. M., Hawley S. L., West A. A., 2004, PASP , 116, 1105 https://doi.org/10.1086/426792 CrossRef Search ADS   West A. A., Hawley S. L., 2008, PASP , 120, 1161 https://doi.org/10.1086/593024 CrossRef Search ADS   West A. A. et al.  , 2004, AJ , 128, 426 https://doi.org/10.1086/421364 CrossRef Search ADS   West A. A., Hawley S. L., Bochanski J. J., Covey K. R., Reid I. N., Dhital S., Hilton E. J., Masuda M., 2008, AJ , 135, 785 https://doi.org/10.1088/0004-6256/135/3/785 CrossRef Search ADS   West A. A., Weisenburger K. L., Irwin J., Berta-Thompson Z. K., Charbonneau D., Dittmann J., Pineda J. S., 2015, ApJ , 812, 3 https://doi.org/10.1088/0004-637X/812/1/3 CrossRef Search ADS   Wilson O. C., 1968, ApJ , 153, 221 https://doi.org/10.1086/149652 CrossRef Search ADS   Wilson O. C., 1978, ApJ , 226, 379 https://doi.org/10.1086/156618 CrossRef Search ADS   Wright N. J., Drake J. J., 2016, Nature , 535, 526 https://doi.org/10.1038/nature18638 CrossRef Search ADS PubMed  Wright N. J., Drake J. J., Mamajek E. E., Henry G. W., 2011, ApJ , 743, 48 https://doi.org/10.1088/0004-637X/743/1/48 CrossRef Search ADS   Wu Y. et al.  , 2011, RA&A , 11, 924 https://doi.org/10.1088/1674-4527/11/8/006 Yang X. L., Chen Y. Q., Zhao G., 2015, AJ , 150, 158 https://doi.org/10.1088/0004-6256/150/5/158 CrossRef Search ADS   Zacharias N., Finch C. T., Girard T. M., Henden A., Bartlett J. L., Monet D. G., Zacharias M. I., 2013, AJ , 145, 44 CrossRef Search ADS   Zhao G., Chen Y.-Q., Shi J.-R., Liang Y.-C., Hou J.-L., Chen L., Zhang H.-W., Li A.-G., 2006, Chin. J. Astron. Astrophys. , 6, 265 CrossRef Search ADS   Zhao G., Zhao Y.-H., Chu Y.-C., Jing Y.-P., Deng L.-C., 2012, RA&A , 12, 723 https://doi.org/10.1088/1674-4527/12/7/002 SUPPORTING INFORMATION Supplementary data are available at MNRAS online. online_tables.tar.gz Please note: Oxford University Press is not responsible for the content or functionality of any supporting materials supplied by the authors. Any queries (other than missing material) should be directed to the corresponding author for the article. APPENDIX A: CHROMOSPHERIC INDICATORS: Ca ii K, H β We measured equivalent widths of Ca ii K and H β lines, EWCaK, and EWHβ, and wavelength regimes for these measurements are listed in Table A1. To get reliable measurements for these two blue lines, we merely measured EWCaK for stars with g-band signal-to-noise ratio SNRg >1.0, and measured EWH β for those with SNRg >3.0. Such SNRg trim criteria cut off about 18 per cent and 3 per cent of total sample stars, for Ca ii K and H β, respectively. Fig. A1 shows the measurements of EWCaK and EWH β as a function of effective temperature. Note that the corresponding basal values are also plotted, which were derived based on a large sample of inactive reference dwarf stars with solar metallicity, following the similar procedure that used to derive mean basal values for H α (see Paper I for more details). Similar to EW$$^{\prime }_{\text{H}\,\alpha}$$, we obtained the excess equivalent widths of the Ca ii K and H β lines, EW$$^{\prime }_{\text{Ca K}}$$ and EW$$^{\prime }_{\text{H}\,\beta}$$, respectively. Further, we derived the excess fractional luminosities, $$R^{\prime }_{\text{Ca K}}$$($$\equiv L^{\prime }_{\text{Ca K}}/L_{\text{bol}}$$) and $$R^{\prime }_{\text{H}\,\beta}$$ ($$\equiv L^{\prime }_{\text{H}\,\beta }/L_{\text{bol}}$$), by using corresponding χ ratios (see Appendix C). These measurements are shown in bottom panels of Fig. A1. Figure A1. View largeDownload slide Emissions of H β and Ca ii K lines for all sample stars. Figure A1. View largeDownload slide Emissions of H β and Ca ii K lines for all sample stars. Table A1. Equivalent width measurements of H α, H β, and Ca ii K lines. Line  Line bandpass (Å)  Pseudo-continua (Å)  H α  6557–6569  6547–6557, 6570–6580  H β  4855–4867  4842–4852, 4873–4883  Ca ii K  3930–3937  3910–3915, 3950–3955  Line  Line bandpass (Å)  Pseudo-continua (Å)  H α  6557–6569  6547–6557, 6570–6580  H β  4855–4867  4842–4852, 4873–4883  Ca ii K  3930–3937  3910–3915, 3950–3955  View Large Note that the derived mean value of inactive reference stars cut off at 6500 K. However, a very small fraction of sample stars are hotter than 6500 K (see Fig. 1). We simply linearly extrapolated the reference value when necessary. In order to get reliable measurements for sample stars with Teff > 6500 K, we merely derived the excess equivalent widths for those stars with Teff less than 6900 K in this work. APPENDIX B: MEASUREMENT UNCERTAINTIES The uncertainties in the equivalent width measurements mainly result from the errors in the RV-corrected LAMOST spectra. The uncertainty in the LAMOST spectra at each wavelength (Poisson noise, and any noises during data reduction) is available and was directly propagated to the equivalent width measurements. Fig. B1 shows the errors (σEWH α) in EWH α for stars in each open cluster, and the increment in errors with decreasing temperature is seen as expected. The typical error of EWH α is less than 0.1 Å for GK-type sample stars, e.g. ∼0.05 Å for Pleiades members with temperatures around 5000 K; however, it become large for M-type stars, e.g. typically greater than 0.1 Å and up to 0.5 Å. Also, the lower SNR in the blue band of spectra increases the uncertainties in the measurements of Ca ii K and H β lines, especially for those M-type stars, e.g. the typical errors of EWH β are 0.05, 0.15, 1.0 Å, for Pleiades members with temperatures around 5500, 4500, 3500 K, respectively. Figure B1. View largeDownload slide Errors in measurements of EWH α. Note Hyades members have smaller errors because of the proximity. Figure B1. View largeDownload slide Errors in measurements of EWH α. Note Hyades members have smaller errors because of the proximity. The spectral resolution of LAMOST spectra (R ∼ 1800) limits the accuracy of RV measurements. To check how RV uncertainties affect EWH α measurements, we have taken PHOENIX model spectra and degraded the resolution to 1800 by convolving a Gaussian profile, found that the differences of EWH α due to RV shifts of ±20 km s−1 are less than 0.02 Å (a relative difference of less than 3 per cent), as shown in Fig. B2. In addition, many cool Pleiades members are fast rotators (Rebull et al. 2016), we thus evaluated the effects of rotational broadening on measurements of EWH α, and found the differences of EWH α due to a vsin i of ∼100 km s−1 are typically less than 0.04 Å (see Fig. B2), which is slightly larger than the errors due to RV uncertainties, but typically less than average errors of EWH α measurements propagated from the uncertainties of the observed flux. So the uncertainties from RV-correction and rotational broadening are marginal compared to those due to uncertainty of observed spectra. Figure B2. View largeDownload slide Top: Estimated differential EWH α values due to wavelength shifts, RV = ±20 km s−1, and rotational broadening, vsin i = 50 and 100 km s−1. Bottom: Relative differences compared to EWH α scales, note that the large values around 4300 K result from the fact that EWH α is near to zero for late K-type stars. Figure B2. View largeDownload slide Top: Estimated differential EWH α values due to wavelength shifts, RV = ±20 km s−1, and rotational broadening, vsin i = 50 and 100 km s−1. Bottom: Relative differences compared to EWH α scales, note that the large values around 4300 K result from the fact that EWH α is near to zero for late K-type stars. The uncertainties in EW$$^{\prime }_{\text{H}\,\alpha}$$ are mainly propagates from the errors of EWH α measurement and corresponding basal value used at a given temperature. The basal value of equivalent width for a given line was derived based on a large sample of reference dwarfs with solar metallicity (e.g. for H α line, see Paper I). The scatter of EWH α of reference stars in each 50 K bin is typically smaller than the errors of measured EWH α, as demonstrated in the upper panels of Fig. 2, where three times of the scatter (±3σ) are shown in grey lines. To estimate the uncertainties on the derived EW$$^{\prime }_{\text{H}\,\alpha}$$ resulting from the flux errors of spectra and uncertainties of temperature, we performed a Monte Carlo simulation 1000 times for six representative EWH α values with Teff from 6000 to 3500 K, using the typical measurement error (see Fig. B1) as the standard deviation of EWH α, and setting 100 K as the standard deviation of temperature, as shown in Fig. B3 and listed in Table B1. It is clear that the EW$$^{\prime }_{\text{H}\,\alpha}$$ uncertainties are dominated by the uncertainties of temperatures for GK-type stars. It is not the same case with M-type stars since the reference EWH α value is not so sensitive to temperature compared to hotter stars. Figure B3. View largeDownload slide Top: Examples of typical values of EWH α with measured error bars σ(EWH α). Middle and bottom: Resulted excess emissions with estimated error bars from different uncertainty sources by using Monte Carlo simulations. The data points at the same temperature are separated by 50 K for display purpose: uncertainties from EWH α (left), from Teff (centre), and from both of them (right). Figure B3. View largeDownload slide Top: Examples of typical values of EWH α with measured error bars σ(EWH α). Middle and bottom: Resulted excess emissions with estimated error bars from different uncertainty sources by using Monte Carlo simulations. The data points at the same temperature are separated by 50 K for display purpose: uncertainties from EWH α (left), from Teff (centre), and from both of them (right). Table B1. Estimated uncertainties from Monte Carlo simulations. Teff (σ)  EWH α(σ)  σ1  σ2  σ3  σ4  σ5  (K)  (Å)  (Å)  (Å)  (dex)  (dex)  (dex)  6000 (100)  −2.0 (0.03)  0.09  0.10  0.08  0.19  0.20  5500 (100)  −1.0 (0.04)  0.12  0.12  0.03  0.08  0.08  5000 (100)  0.0 (0.05)  0.09  0.10  0.02  0.04  0.04  4000 (100)  2.0 (0.10)  0.09  0.13  0.02  0.04  0.05  3500 (100)  4.0 (0.20)  0.06  0.22  0.02  0.06  0.07  4500 (100)  0.0 (0.07)  0.08  0.11  0.04  0.06  0.08  4500 (100)  1.0 (0.07)  0.08  0.11  0.02  0.03  0.04  4500 (100)  2.0 (0.07)  0.08  0.11  0.01  0.03  0.03  Teff (σ)  EWH α(σ)  σ1  σ2  σ3  σ4  σ5  (K)  (Å)  (Å)  (Å)  (dex)  (dex)  (dex)  6000 (100)  −2.0 (0.03)  0.09  0.10  0.08  0.19  0.20  5500 (100)  −1.0 (0.04)  0.12  0.12  0.03  0.08  0.08  5000 (100)  0.0 (0.05)  0.09  0.10  0.02  0.04  0.04  4000 (100)  2.0 (0.10)  0.09  0.13  0.02  0.04  0.05  3500 (100)  4.0 (0.20)  0.06  0.22  0.02  0.06  0.07  4500 (100)  0.0 (0.07)  0.08  0.11  0.04  0.06  0.08  4500 (100)  1.0 (0.07)  0.08  0.11  0.02  0.03  0.04  4500 (100)  2.0 (0.07)  0.08  0.11  0.01  0.03  0.03  σ1: $$\sigma (\text{EW}^{\prime }_{\text{H}\,\alpha}$$) due to σ(Teff). σ2: $$\sigma (\text{EW}^{\prime }_{\text{H}\,\alpha}$$) due to both σ(EWH α) and σ(Teff). σ3: $$\sigma (\log {{R}}^{\prime }_{\rm{H}\,\alpha })$$ due to σ(EWH α). σ4: $$\sigma (\log {{R}}^{\prime }_{\rm{H}\,\alpha })$$ due to σ(Teff). σ5: $$\sigma (\log {{R}}^{\prime }_{\rm{H}\,\alpha })$$ due to both σ(EWH α) and σ(Teff). View Large The excess fractional luminosity $$R^{\prime }_{\text{H}\,\alpha}$$ depends on both derived EW$$^{\prime }_{\text{H}\,\alpha}$$ and estimated value of χ at a given temperature. To check the temperature selection effect on derived $$R^{\prime }_{\text{H}\,\alpha}$$, we performed a simple Monte Carlo simulations (see Table B1), which show that temperature uncertainty dominates the derived $$R^{\prime }_{\text{H}\,\alpha}$$ uncertainty, as shown in Fig. B3. APPENDIX C: THE RATIO OF THE CONTINUUM FLUX TO THE BOLOMETRIC FLUX To quantify the strength of CA as the ratio between H α emission luminosity and bolometric luminosity, Walkowicz, Hawley & West (2004) introduced a distance-independent value χ, the ratio between the apparent continuum flux near H α, f(λ6560), and the apparent bolometric flux, fbol, namely, χ = f(λ6560)/fbol. Then the measured equivalent width of H α emission line could be converted to the fractional luminosity by using the χ value, i.e. LH α/Lbol = χ × EWH α. χ-method is widely used where flux-calibrated spectrum and distance information is not available for quantifying CA level (e.g. Douglas et al. 2014; West et al. 2015). χ values could be estimated based on observed spectra (e.g. West & Hawley 2008) or model spectra (e.g. Reiners et al. 2012). In this work, we estimated the χ values as the ratio of the surface continuum flux, F(λc), to the stellar surface bolometric flux, Fbol, using the following formula:   $$\chi = \frac{F(\lambda _{c})}{F_{\text{bol}}} = \frac{F(\lambda _{c})}{\sigma T_{\text{eff}}^4},$$ (C1)where σ is Stefan–Boltzmann constant, σ ≈ 5.6704 × 10−5 erg s−1 cm−2 K−4. We evaluated F(λc) for the lines of interest using PHOENIX ACES model spectra with solar metallicity (Husser et al. 2013). We then obtained χ values for H α, Ca ii K, and H β lines, as shown in Fig. C1. Figure C1. View largeDownload slide χ values of H α, H β, and Ca ii K lines, χH α, χH β, and χCaK, respectively, derived based on SDSS spectra, and PHOENIX model spectra with log g = 4.5 (black solid lines) and log g = 5.0 (black dashed lines). The red long dashed line denotes the relation derived from empirically determined H α line surface flux provided by Soderblom et al. (1993). Figure C1. View largeDownload slide χ values of H α, H β, and Ca ii K lines, χH α, χH β, and χCaK, respectively, derived based on SDSS spectra, and PHOENIX model spectra with log g = 4.5 (black solid lines) and log g = 5.0 (black dashed lines). The red long dashed line denotes the relation derived from empirically determined H α line surface flux provided by Soderblom et al. (1993). To check the agreement of these χ values derived from models with those from the observations, we collected the flux calibrated SDSS spectra of stars in NGC 2420, M 67, and Praesepe. 140 members of NGC 2420 with log g > 4.0 and 52 dwarfs in M 67 were selected from the catalogue of Lee et al. (2008). We estimated their radii using the temperatures from Lee et al. (2008), adopting the PARSEC model (Chen et al. 2014) with a metallicity of [Fe/H] ∼ −0.2 (Jacobson, Pilachowski & Friel 2011) and an age of ∼2 Gyr (Salaris, Weiss & Percival 2004), and the PARSEC model with a solar metallicity (Jacobson et al. 2011) and an age of ∼4 Gyr (Önehag et al. 2011; Geller, Latham & Mathieu 2015), for NGC 2420 and M 67, respectively. We collected SDSS spectra with good quality (SNR > 10) for 36 Praesepe stars of Douglas et al. (2014) based on the SDSS DR12 SSPP output catalogue (Lee et al. 2008; Smolinski et al. 2011). We estimated their effective temperatures using the colours (V − Ks and/or r − Ks), and derived corresponding radii adopting the PARSEC model with a metallicity of [Fe/H] ∼ +0.16 (Carrera & Pancino 2011) and an age of about 650 Myr. The distances adopted from WEBDA (http://www.univie.ac.at/webda/) for NGC 2420, M 67, and Praesepe are 3085, 908, and 187 pc, respectively. Based on the observed continuum flux, f(λc), and their corresponding radii (R) and distances (D), the χ ratios could be derived using the following relation:   $$\chi = \frac{F(\lambda _{c})}{F_{\text{bol}}} = \frac{(D/R)^{2}f(\lambda _{c})}{\sigma T_{\text{eff}}^4},$$ (C2)as shown in Fig. C1, wherein the error bar is equivalent to its uncertainty due to 10 per cent error of distance. There exists a slight offset between members in NGC 2420 and M 67 at the same effective temperatures, which may be due to the uncertainties in distance, and/or the difference in basic physical parameters between these two clusters. Also shown in Fig. C1 by the red dashed line is a relation for H α line, which was derived from empirically determined surface flux of H α line provided by Soderblom et al. (1993). As shown by the figure, the χ values from model spectra are overall in good agreement with those from observations, there however exist some small discrepancies between them, e.g. ∼0.1 dex in H α line around 4500 K. APPENDIX D: ESTIMATE OF QUIESCENT PHOTOSPHERE TEMPERATURE AND MASS The cool spots on the stellar surface contribute light that would affect spectral features and broad-band colours, for which the derived effective temperature becomes the average surface temperature (flux-weighted mean temperature) in spotted stars. Therefore, no direct way to derive a quiescent temperature for active stars, particularly those Pleiades members with large spot coverages (see Paper I for more discussion on this topic). For Pleiades members, we followed the procedures in Paper I to estimate the quiescent temperature (see Paper I for more details). For M34 members, we collected V-band magnitudes from literature (Jones & Prosser 1996; Irwin et al. 2006; Meibom et al. 2011), and got their Ks-band magnitudes from 2MASS archive (Skrutskie et al. 2006) via Vizier x-match service. We estimated their quiescent temperatures based on V − Ks colours using PARSEC isochrones (Chen et al. 2014) with an age of 220 Myr and solar metallicity, wherein we adopted the reddening of E(B − V) = 0.07 mag (Canterna, Crawford & Perry 1979). For Praesepe and Hyades candidates, first we derived their effective temperatures based on V − Ks colours, where V magnitudes taken from the fourth U.S. Naval Observatory CCD Astrograph Catalogue (Zacharias et al. 2013). A small fraction of sample stars have no V magnitudes, we then derived their temperatures based on r − Ks, where r magnitudes from the Carlsberg Meridian Catalogue 15 (Copenhagen University et al. 2011). The PARSEC models with an age of 650 Myr and a metallicity of [Fe/H] ∼ 0.1 were used to convert colours to temperatures, where the reddening $$E(B-V) = 0{^{\rm m}_{.}}027$$ and $$E(B-V)=0{^{\rm m}_{.}}001$$ (Taylor 2006) were adopted for Praesepe and Hyades, respectively. We then compared the colour-based temperatures with the spectra-based temperatures (Tsp, for FGK-type stars, we adopt the values provided by LASP, for M-type stars, we derived it via the molecular band CaH using the method in Paper I), and found that there exist complex differences between colour-based temperatures and spectra-based temperatures, as shown in Fig. D1. We finally adopted the effective temperatures from spectral features as the quiescent temperatures for stars in these two clusters. Figure D1. View largeDownload slide The differences between V − Ks-based temperature (TVK), r − Ks-based temperature (TrK), and the temperature from spectral features (Tsp) among stars in Praesepe and Hyades. Figure D1. View largeDownload slide The differences between V − Ks-based temperature (TVK), r − Ks-based temperature (TrK), and the temperature from spectral features (Tsp) among stars in Praesepe and Hyades. For the M-type slow rotators in Kepler field, first we estimated temperature for each target based on the molecular bands of CaH and TiO following the method in Paper I, and then took the average of our measured value and the value provided by McQuillan et al. (2014), as finally adopted quiescent temperature. To identify potential binaries and non-members for the member candidates in open clusters, we followed the procedures used in Paper I by using available CMDs, e.g. we identified these stars as photometrically binaries/non-members that have larger deviation compared with empirical single-star loci in V versus V − Ks diagram, namely, brighter stars with $$-1{^{\rm m}_{.}}0 < \Delta V < -0{^{\rm m}_{.}}5$$ were classified as probable binary members, while the fainter stars with $$\Delta V >0{^{\rm m}_{.}}3$$ or very much brighter stars with $$\Delta V < -1{^{\rm m}_{.}}0$$ were identified as probable non-member stars. Once the effective temperatures are obtained for each star, we estimated the mass from temperature based on mass–Teff relation from PARSEC isochrones (Chen et al. 2014). Note that, for these M-stars in Kepler field, we used the isochrones with a solar metallicity and an age of 2 Gyr. APPENDIX E: ROTATION PERIODS COLLECTION For Pleiades, Hartman et al. (2010) detected rotation periods for nearly 400 members using the HATNet transit survey data. Recent studies based on dedicated photometric surveys, e.g. the Palomar Transient Factory (PTF) (Covey et al. 2016), and the K2 (Rebull et al. 2016), increased the measurements of rotation periods for more than 500 new members of this open cluster, vastly expanding the number of Pleiades members with rotation periods. For these members with multirotation measurements, we adopted values from the K2 (Rebull et al. 2016), and for others, we adopted their rotation periods provided by Covey et al. (2016) and Hartman et al. (2010). Totally, 231 Pleiades members in our sample have available rotation periods. For M34, we collected rotation periods from literature, Irwin et al. (2006), James et al. (2010), and Meibom et al. (2011). The latest measurements were adopted in this paper. Totally, we found 21 members have rotation periods. For stars in Praesepe, their rotation periods were collected from Scholz & Eislöffel (2007), Delorme et al. (2011), Scholz et al. (2011), Agüeros et al. (2011), and Kovács et al. (2014). Totally, we got their rotation periods for 76 stars in Praesepe sample. For Hyades, we collected the rotation periods from Radick et al. (1987), Delorme et al. (2011), Hartman et al. (2011), and Douglas et al. (2016). Totally, 24 stars in Hyades sample have rotation periods. APPENDIX F: ON-LINE TABLES We compiled the measurements for stars in Pleiades, M34, Praesepe, and Hyades in on-line Table F1–F4, respectively, including effective (quiescent) temperatures, equivalent widths of H α, H β, and Ca ii K lines, their excess equivalent widths, excess fractional luminosities, and spot filling factors. © 2018 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society

### Journal

Monthly Notices of the Royal Astronomical SocietyOxford University Press

Published: May 1, 2018

## You’re reading a free preview. Subscribe to read the entire article.

### DeepDyve is your personal research library

It’s your single place to instantly
that matters to you.

over 18 million articles from more than
15,000 peer-reviewed journals.

All for just $49/month ### Explore the DeepDyve Library ### Search Query the DeepDyve database, plus search all of PubMed and Google Scholar seamlessly ### Organize Save any article or search result from DeepDyve, PubMed, and Google Scholar... all in one place. ### Access Get unlimited, online access to over 18 million full-text articles from more than 15,000 scientific journals. ### Your journals are on DeepDyve Read from thousands of the leading scholarly journals from SpringerNature, Elsevier, Wiley-Blackwell, Oxford University Press and more. All the latest content is available, no embargo periods. DeepDyve ### Freelancer DeepDyve ### Pro Price FREE$49/month
\$360/year

Save searches from
PubMed

Create lists to

Export lists, citations