Open Access
Issue
A&A
Volume 642, October 2020
Article Number A225
Number of page(s) 10
Section Stellar atmospheres
DOI https://doi.org/10.1051/0004-6361/202038926
Published online 23 October 2020

© E. M. Amazo-Gómez et al. 2020

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Open Access funding provided by Max Planck Society.

1 Introduction

Rotation periods in cool main-sequence stars can be traced by observing the brightness modulation caused by the presence of active regions on stellar surfaces. Those active regions are generated by the emergence of strong localized magnetic fields approximately described by flux tubes (see e.g. Solanki 1993). Large flux tubes form dark spots, while ensembles of smaller flux tubes form bright faculae (see, e.g. Solanki et al. 2006, for a detailed review of the solar case). The active regions usually consist of a sunspot group surrounded by faculae. The transits of such active regions over the visible disk as the star rotates would cause brightness variability. Consequently, the stellar light curves (LCs) contain information about both the rotation periods and the properties of active regions. However, retrieving this information from the LCs often proves a daunting task (see e.g. Basri 2018).

The Kepler mission (Borucki et al. 2010) has provided the community with records of photometric observations with unprecedented precision and cadence. The Kepler LCs have been widely used to determine stellar rotation periods (e.g. Walkowicz & Basri 2013; Reinhold & Gizon 2015; Nielsen et al. 2013; García et al. 2014; McQuillan et al. 2014; Buzasi et al. 2016; Angus et al. 2018; Santos et al. 2019). Despite considerable success in determining the rotation periods of many fast-rotating and active stars (see, e.g. McQuillan et al. 2014, who published rotation periods of about 34030 stars identified as being located on the main sequence) there is a lack information on periods of slowly rotating stars, that is, stars with near-solar and longer rotation periods. For example, the rotational period of the Sun may not be detectable during intermediate and high levels of solar activity (see Lanza & Shkolnik 2014; Aigrain et al. 2015).

The difficulties in detecting periods of slowly rotating stars might be an important contribution to the explanation of lower-than-expected numbers of G-type stars with near-solar rotation periods (van Saders et al. 2019). The difficulty in obtaining a reliable measurement of the rotation periods of stars with variability patterns similar to that of the Sun can also affect solar-stellar comparison studies (see e.g. Witzke et al. 2020; Reinhold et al. 2020).

Within this context, we have developed a method aimed at determining rotation periods of low-activity stars like the Sun. In Shapiro et al. (2020) (hereinafter, Paper I), we found that the power spectra of brightness variations of such stars are strongly affected by the evolution of active regions. In particular, the rotation peak can be significantly weakened or it may even disappear from the power spectrum if the lifetimes of starspots are too short. Furthermore, the delicate balance between spot and facular contributions to the variability might lead to the appearance of spurious peaks, which do not correspond to the rotation period but could be easily mistaken for one (see also Shapiro et al. 2017).

In Paper I, we showed that the high-frequency tail of the power spectrum is much less affected by the evolution of magnetic features than frequencies near the rotation period. Consequently, we proposed to use information in the high-frequency tail for the determination of stellar rotation periods. In particular, we suggested that the period, PHFIP, corresponding to the maximum of the gradient of the power spectrum (GPS) (i.e. to the inflection point) in the high-frequency tail could be used to identify the stellar rotation period, Prot, via the simple scaling relation:  Prot=PHFIP/α.\begin{equation*}~P_{\textrm{rot}}=P_{\textrm{HFIP}} / \alpha. \end{equation*}(1)

Here, α is a calibration factor which is independent of the evolution of active regions. It shows only a very weak dependence on the stellar inclination. For example, the inclination dependence can be neglected for inclinations of 45° and greater;see Fig. 9 from Paper I. Statistically, this corresponds to roughly 70% of stars.

The model developed in Paper I indicated that the value of α shows a moderate dependence on the ratio between facular and spot areas of the individual active regions at the moment of emergence, SfacSspot. This ratio was assumed to be the same for all active regions (see a detailed discussion in Paper I). The dependence of the inflection point position on the facular-to-spot area ratio leads to a certain degree of uncertainty (up to 25%) when determining stellar rotation periods since the value of SfacSspot for a given star is a priori unknown. At the same time, it allows us to retrieve valuable information about facular versus spot-dominated regimes with regard to the variability of stars where rotation periods can be determined using other methods (see Amazo-Gómez et al. 2020 hereafter, Paper II).

A first test of the gradient of the power spectrum method (hereafter, GPS) was performed in Paper II, where we applied it to solar brightness variations. We showed that in contrast to other methods, GPS allows for an accurate determination of the solar rotation period at all levels of solar activity. Additionally, we analysed time intervals where solar variability was spot-dominated and when it was faculae-dominated. We showed that these regimes can be distinguished in the GPS profile due to the substantially different centre-to-limb variations of faculae and spots.

In this study, we apply the GPS method to stars with determined rotation periods from Kepler photometry. The goal here is twofold: firstly, we test the GPS method further before applying it to stars with unknown rotation periods; secondly, we investigate whether the α factor and, consequently, the facular or spot composition of stellar active regions, is dependent on the rotation period. In Sect. 2, we describe the stellar sample we used. In Sect. 3, we present the main results. Our conclusions are summarized in Sect. 4.

Table 1

Stellar parameters for samples A and B.

2 Stellar sample selection

In this study, we consider stars in the field of view (FOV) of the Kepler telescope for which McQuillan et al. (2014) managed to determine rotation periods using the auto-correlation function (ACF). To ensure that the main source of the variability for the selected stars is magnetic activity, we only selected stars on the main-sequence, using Teff and log g values from the Huber et al. (2014) catalogue to exclude giants (see Table 1). We note that Huber et al. (2014) calibrated effective temperatures to the infrared flux temperature scale. This resulted in an approximately 200 K offset from the original Kepler Input Catalogue (KIC) (Pinsonneault et al. 2012). We also precluded stars flagged in the KIC as giant (GS), eclipsing binary (EB), or host stars with planetary transits confirmed (PTC), with planet candidates (PC), and false-positive planets (FP).

We selected two sets of stars with near-solar parameters. The selection criteria for both samples (A and B) are illustrated in Fig. 1 and given in Table 1. Figures 1 and 2 show the variability ranges of the set of selected stars. Kepler observatory provided 4 yr of photometric information, from 2009 to 2013, segmented in 18 quarters (Q0Q17) due to the telescope re-orienting itself every 90 days. The Kepler observing quarters resulted in Q0 of 10 days, and Q1 of 33 days for the commissioning phase and, segmented 90-day LCs for Q2 to Q16 (see public data release 25, Thompson et al. 2016; Van Cleve & Caldwell 2016). The second month of Q17 was terminated after less than five days of observation, following the failure of reaction wheel 4.

Sample A (see Fig. 1, panel II) was selected to test the performance of the GPS method for stars with near-solar rotation periods. Hence, in this sample, we considered stars with a narrow range of rotational periods between 24.0 and 27.4 days (i.e. encompassing the sidereal Carrington rotation period of the Sun at 25.4 days) and a broad range of effective temperatures Teff ∈ (5000–6000) K. These selection criteria yielded a sample consisting of 686 stars. Among this sample, 282 stars also have rotation periods obtained by Reinhold & Gizon (2015) using the generalized Lomb–Scargle periodograms (hereafter, GLS).

Sample B (see Fig. 1, panel III) was selected to study the dependence of the inflection point position on the rotation period. Therefore, in contrast to sample A, we considered stars with a broad range of rotation periods (between 15and 40 days) but a narrow range of effective temperatures (5700–5900 K, encompassing the solar value of 5778 K). These criteria led to the selection of 361 stars for sample B after removing the overlapping targets the with initial sample A. Rotation periods of 172 stars in this sample were also reported in Reinhold & Gizon (2015).

Between our two samples, we thus considered 1047 Kepler stars in all. The LCs were acquired in the long-cadence mode (i.e. with a cadence of 29.42 min). Following McQuillan et al. (2014) and Reinhold & Gizon (2015), we utilized LCs from Q1Q14 processed with the pre-search conditioning and Bayesian maximum a posteriori approach (PDC-MAP, see Smith et al. 2012). For quarters Q15Q17, only processing with multiscale MAP (PDC-msMAP Stumpe et al. 2014) is available.

In Fig. 2, we plot the distribution of variability ranges in our samples A and B. These variability values are defined by computing the difference between the 95th and 5th percentiles of the sorted flux values for each of the Kepler observing quarters (see Basri et al. 2011) and then taking the median value among the quarters. This defined variability range was chosen versus the approaches based on the standard deviation analysis by Mathur et al. (2014); He et al. (2015) or the smoothed amplitude (10th to 90th) method presented in Douglas et al. (2017), given the higher range of amplitude used in Basri et al. (2011). The selection of the methods mentioned are not expected to compromise the analysed outcome. Additionally, we show solar variability ranges computed using total solar irradiance data (TSI, i.e. total radiative flux from the Sun at 1 A.U.) for 1996–2017 obtained by the Variability of solar IRradiance and Gravity Oscillations (VIRGO; Fröhlich et al. 1997) experiment on the SOlar and Heliospheric Observatory SoHO mission. For these VIRGO data, the entire 1996–2017 observation period was split into 6 Kepler-like time ranges (five 1530-day periods and one 712-day period for a total of 7787-days starting 28 January 1996; see Fig. 3 from Paper II). The solar variability value for each of the time ranges is represented in Fig. 2 as vertical green dashed lines. This gives a range of solar variability of Var∈ (400–1300) ppm.

Figure 2 shows that most of the stars in our samples are much more variable than the Sun. This agrees with García et al. (2014); Buzasi et al. (2016); Reinhold et al. (2020), who showed that solar-type stars (i.e. stars with near-solar fundamentalparameters and rotation periods) are, on average, significantly more variable than the Sun. Furthermore, our samples A and B also contain stars which are cooler and rotate faster than the Sun. These stars are also expected to be more variable than the Sun (seee.g. McQuillan et al. 2014, for a discussion of the dependence of the variability on the rotation period and temperature).

We note that an anomalous low variability of solar-type stars found by Reinhold et al. (2020) does not necessarily imply that the Sun is an outlier. An alternative explanation is that by comparing solar variability to the sample of stars with known rotation periods, we focus only on a small sub-sample of stars for which the ACF method could return rotation periods (and the Sun most probably would not belong to such a sample). Along this line of reasoning, Reinhold et al. (2020) found that solar levels of photometricvariability are typical for stars having near-solar fundamental parameters but unknown rotation periods.

thumbnail Fig. 1

Panel I: temperature–rotation diagram for a sample of 34030 stars (coloured circles indicating the variability range) with rotation periods determined by McQuillan et al. (2014), along with 55501 stars where the authors found a period but evaluated it as not significant (beige dots; see panel III for a better visualization). Panel II: sample of 34030 stars with determined rotation periods in grey and 55501 stars with non-significant rotation period determinations in beige. For panels II and III, only the stars from sample A (panel II) and sample B (panel III) are shown in colour; see Table 1 for the properties of samples A and B. Panel II illustrates the stellar sample A, selected by near-solar rotation period and temperatures from 5000 to 6000 K. Panel III illustrates stellar sample B, which contains stars with near-solar effective temperature and a broad range of rotation periods. The Sun is represented by the solar symbol ⊙.

thumbnail Fig. 2

Histograms for variability ranges from samples A (violet outline) and B (blue rectangles). Vertical green dashed lines represent solar variability calculated for six activity seasons using 21 yr of VIRGO TSI data from Paper II (see text for more details).

Table 2

GPS outcome values.

3 Results and discussion

In this section, we calculate the position of the inflection point for each star in the samples A and B defined in Sect. 2. Following the methodology described in Papers I and II, we first calculate the power spectra of the stellar brightness variations using a Paul wavelet on the order of six (see Torrence & Compo 1998) for the Kepler observing quarters Q1Q17. We determined the period corresponding to the high-frequency inflection point, PHFIP_(Q_n), per quarter and calculate the mean value for PHFIP over all 17 quarters for each star. This allows us to obtain a unique representative value of PHFIP per star. The uncertainty is calculatedusing 2-σ of the distribution of the obtained PHFIP values. Finally, we used the PHFIP to calculate thestellar rotation period, Prot (see Table 2 and on-line reference for a compilation of GPS outputs and comparison with GLS and ACF reference values).

In Fig. 3, we plot the mean values of the PHFIP positions for each of the stars against the rotation periods from Reinhold & Gizon (2015) (left panel, GLS) and McQuillan et al. (2014) (right panel, ACF). The rotation periods and positions of the inflection points are well-correlated. A linear fit constrained to go through the origin of the coordinate system gives PHFIP = 0.19 × Prot with Pearson coefficients for periods of 0.81 from Reinhold & Gizon (2015) and 0.80 from McQuillan et al. (2014).

The scatter around the linear fits has multiple sources. First, the calibration coefficient between rotation period and inflection point, α = PHFIPProt, depends on the relative roles that bright faculae and dark spots play in generating stellar brightness variations. According to the model presented inPaper I, these roles are regulated by the ratio between facular and spot areas of active regions at the time of emergence, SfacSspot (i.e. zero ratio would lead to a purely spot-dominated star, while very large ratios would correspond to a faculae-dominated star). Secondly, there is an intrinsic statistical uncertainty of the GPS method. For example, in Paper I, we found that even for a star with a fixed SfacSspot ratio, the factor α showed 5–10% variations from one realization of active regions emergence to another. Finally, there is also an uncertainly in the determination of rotation periods by Reinhold & Gizon (2015) and McQuillan et al. (2014) (see e.g. Fig. A.2, where we compare the periods from these two sources for the 172 stars of sample B that are common to both).

In Fig. 4, we show calibration factors, α, for samples A (top panel) and B (bottom panel). The rotation periods of stars in both samples are taken from McQuillan et al. (2014). In Paper I, we demonstrated that the profile of the high-frequency tail of the power spectrum and, consequently, the values of α depend on the centre-to-limb variations (CLVs) of the brightness contrasts of magnetic features. Since spots and faculae have different CLVs, the value of α depends on their relative contributions to the stellar brightness variations. For the extreme cases, we found that α is about 0.14 for simulated stellar light curves with variability solely determined by faculae and about 0.21 for simulated stars with variability dominated by spots. These values are designated, respectively, by the red and green horizontal dashed lines in Fig. 4. It is reassuring to see that most of the α values for samples A and B appear between these two extreme-cases. Stars with values of α outside of this range (in particular, with α > 0.21) are likely due to: inclination angles below 45°, which can lead to a shift of the inflection point to lower frequencies (see Fig. 9 from Paper I); statistical noise of the GPS method; and possible uncertainties in rotation periods from McQuillan et al. (2014).

For sample A, the ratios are shown as a function of stellar effective temperature from Huber et al. (2014), while for sample B, they are plotted as a function of stellar rotation period from McQuillan et al. (2014). The upper panel of Fig. 4 shows that for near-solar rotation periods (the rotation periods in sample A were constrained between 24 and 27.4 days; see Table 1), the position of the inflection point shows no significant dependence on the effective temperature (e.g. the fitting of a slope gives a value of 7.36 × 10−7, which is well below the 1 σ uncertainty of 1.8 × 10−6). We also note that the mean value of α =  0.19 is equal to the slope of the regression shown in Fig. 3. This implies that neither the SfacSspot value nor CLVs of facular and spot contrast change significantly within the 5000–6000 K domain of sample A. We note, however, that we cannot conclusively exclude the improbable scenario that the effect from the change of the facular and spot contributions to brightness variability on α is compensated by a change of facular and spot CLVs, such that the net effect on the inflection point is very small.

The bottom panel of Fig. 4 shows that for stars with near-solar effective temperatures there is a rather weak but statistically significant dependence of the α factor on the rotation period. For example, fitting a linear dependence returns a slope value of 9.3 × 10−4 which is 3.8 times larger than its 1σ uncertainty of 2.5 × 10−4. However, the value of the slope is strongly affected by a couple of slowly rotating stars and, thus, might not represent a trend in the full sample. To better characterize such a trend, we calculated the mean value of the calibration factor in several bins of rotation period values. We compiled the mean α values per several bins of rotation periods; see Table 3 for details. To further illustrate the trend of α values with rotation period, the histogram to the right side of the panel shows the distributions of α values for two rotation periods – one for stars with rotation periods below 21 days and another with rotation periods above 21 days. We can see that the two distributions are clearly shifted relative to each other and the α-values of faster rotating stars are larger than those slow-rotating stars.

We note that the “n” number of stars and the amplitude of photometric variability in our samples decreases with rotation period. Consequently, slow rotators might be more affected by photometric noise. We investigated the possible effect of Kepler white noise on the deduced positions of inflection points for the stars in our samples. In Fig. 5, we plot the dependence of the α factor values on the expected Kepler noise levels for each of the stars, calculating the amplitude of the Kepler noise as a function of the Kepler magnitude (following Lammer 2013). The derived precision (called the noise in the context of Kepler) accounts for noise introduced by the instrument and gives it as a function of the Kepler magnitude of the source and the variability of sources (see Fig. A.4).

As much as 99.9% of the stars in our sample present a Kepler magnitude of 16 mag or fainter. We find that values of the α factor are independentof the Kepler noise, with fits of a linear dependence to samples A and B giving slope values well below their 1σ uncertainties (7.4 × 10−7 and 4.9 × 10−8, respectively). Consequently, we do not expect the Kepler noise to affect the positions of inflection points determined for stars in our samples.Furthermore, we note that photometric noise would shift the position of the inflection point to lower frequencies (see Paper II), meaning that it would lead to a trend that is opposite to what we see in the bottom panel of Fig. 4.

A possible explanation of the observed tendency in Fig. 4 is a change in the relative contribution of faculae and spots to stellar rotation variability (or SfacSspot ratio in terms of Paper I) with rotation period. The increase of the α factor with rotation rate implies that the SfacSspot ratio (and, consequently, the contribution of faculae to the rotational variability) is lower in faster rotating and, therefore, more active stars. Such a trend is consistent with an extrapolation to higher activities of observed solar behaviour. Indeed, the mean size of spots on the Sun increases during periods of high solar activity (Hathaway 2015; Mandal et al. 2020). At the same time the SfacSspot ratio decreases with the size of active regions and their spot components. An extrapolation of these trends to activity levels higher than seen in the Sun results in an increase of the α factor with activity, and, consequently, with rotation rate, as indicated by the bottom panel of Fig. 4.

We note that the ratio SfacSspot between facular and spot areas of the individual magnetic features at the moment of their emergence discussed until now is different from the ratio between “instantaneous” stellar disk coverage by faculae and spots. The former is a property of a magnetic feature during its emergence onto the surface of the star, while the latter is strongly affected by the evolution of the magnetic flux after emergence. For example, in the hypothetical case of facular portions of active regions evolving exactly as spot portions, these two ratio remain the same. In reality, the “instantaneous” ratio is generally significantly larger than that “at the time of emergence” since faculae live longer than spots.

Solar observations show that the ratio between such instantaneous solar-disk coverage by faculae and spots decreases as solar activity increases (Chapman et al. 1997; Foukal 1998). The observed patterns of stellar-brightness variability indicate that this trend also extends to activity values that are significantly higher than those observed on the Sun (Shapiro et al. 2014). Our result indicates that not only the ratio between instantaneous facular and spot disc coverage shows this trend. Also facular to spot area ratio corresponding to individual active regions “at the time of emergence” continues to decrease with increasing level of activity, also beyond the level of solar activity observed until now. We note that this result is not a simple consequence of the drop in the instantaneous ratio. Simulations with a surface flux transport model by Cameron et al. (2010) show that the origin of the decrease in the instantaneous ratio with increasing activity is rather complex. It is, to a great extent, caused by a stronger cancellation of small-scale magnetic field associated with faculae. Consequently, it does not necessarily demand any changes in the structure of the emerging magnetic flux which defines the ratio corresponding to the individual active region at the time of emergence; see the discussion in Paper I.

The bottom panel of Fig. 4 shows that the dependence of α on the rotation period is quite noisy, that is, there is quite a large spread of values for a fixed rotation period. This spread basically covers the entire range of values between faculae- and spot-dominated variability. In particular, it is significantly larger than statistical noise in the inflection point position that we found in Paper I. We speculate that such a large spread implies that the SfacSspot ratio is not uniquely defined by the stellar effective temperature and rotation period.

In Papers I and II, we found that the solar value of the calibration factor (αSun = 0.158) is closer to the faculae-dominated case (α = 0.14) than to the spot-dominated case (α = 0.21). Interestingly, Fig. 4 shows that the solar (αSun) value appears to be rather low relative to that of stars in both of our samples (see also Fig. 3, where the Sun is clearly below the regression line). This is, however, not surprising since most of the stars in our samples are significantly more variable than the Sun even though we selected the stellar sample by the reported detected rotation period and not by the variability (see Fig. 2, cf. Reinhold et al. 2020). This implies that these stars are also more active than the Sun (see also Zhang et al. 2020, who showed that stars with known near-solar rotation periods have systematically higher values of S-index than the Sun). Therefore, we can expect that their SfacSspot ratios are smaller and, consequently, their α factors are larger.

thumbnail Fig. 3

Position of the high-frequency inflection point (PHFIP) plotted against rotation period. Only sample B is shown here. Rotation periods are taken from Reinhold & Gizon (2015) (left panel) and McQuillan et al. (2014) (right panel). Colours represent the stellar effective temperature, Teff. The Sun is represented by the solar symbol ⊙. Dashed lines in both panels indicate a linear fit constrained to go through the origin of the coordinate system. A logarithmic visualization is available in the appendix, see Fig. A.1.

thumbnail Fig. 4

Top panel: α factor versus Teff for sample A exhibits consistency across a broad temperature range. The black line corresponds to the linear fit to values with an uncertainty within 2-σ of the mean of the distribution as shown by coloured squares (grey squares lie outside the 2-σ of the distribution). The histograms to the right side of the panel display the distribution of α values for two effective temperature regimes, using Teff = 5500 K as a threshold. Bottom panel: α factor versus Prot from McQuillan et al. (2014) for sample B shows a slight decrease in α with rotation period. The coloured segments indicate the mean of α for the different Prot ranges as indicated in Table 3. The histograms to the right of the panel indicate the distribution of α values for two rotation period regimens, using 21 d as a threshold. For both panels, the error bars represent 2-σ uncertainties of the α values over all Kepler quarters available for each star. The grey squares lie outside of a 2-σ of the distribution. The dashed red and green horizontal lines represent the α factor values in the extreme cases with all variability being due to spots (α = 0.21) and all due to faculae (α = 0.14), respectively1.

Table 3

Mean α-values in sample B per bin.

thumbnail Fig. 5

α factor versus photometric precision for sample A (left panel) and sample B (right panel) for Kepler observations. Those records for single long observation of isolated stars observed in an uncrowded pixel, i.e. light curves of resolved targets without contamination of additional sources. Extreme-case limits for spot- and faculae-dominated stars are shown as horizontal dashed linesin red and green, respectively. The individual data points are coloured according to the detected variability range for that particular star. This is consistent with the preceding figure, Fig. 4. The error bars represent 2-σ uncertainties of the α values distribution over all Kepler quarters available per star. Gray rhomboids represent data points that lie more than 2-σ from the centre of the distribution.

4 Summary

In this study, we developed the GPS method, which is a novel means for determining stellar rotation periods from photometric time series. Instead of basing this determination on the more traditional means of identifying the strongest peak in a Lomb–Scargle periodogram or a maximum of the auto-correlation function, we identify the steepest point (i.e. the inflection point) in the global wavelet power spectrum of stellar brightness variations. In Paper II, we showed that while the solar brightness contributions from faculae and spots can oppose each other to reduce any peak due to the rotation period from Lomb–Scargle periodograms and auto-correlation functions, it has only a very minor effect on the location of this high-frequency inflection point and the resulting ratio between the period corresponding to the inflection point and the actual rotation period, α. In Paper I, the factor ofα, however, shows a moderate dependence on the relative contribution of faculae to stellar brightness variations. Therefore, identifying the position of the inflection point allowed determination of the rotation period in stars where other methods fail (with an internal uncertainty of about 25%). At the same time, this GPS method allows us to assess the relative role of faculae in stars with known rotation periods. In Paper II, we tested the performance of the GPS method against solar photometric data. We demonstrated that in contrast to other methods, the GPS method allows for an accurate determination of the solar rotation period independently of the solar activity level.

In this study, we applied the GPS method to 1047 F-, G- and K-type stars with rotation periods as reported in McQuillan et al. (2014). We show that the position of the high-frequency inflection point is well-correlated with the rotation periods of stars in the two samples we analysed, providing further validation of the GPS method. We emphasize that the stellar light curves analysed in this study and the solar light curves analysed in Paper II are quite different: the amplitudes of brightness variability in the stellar samples in this study are generally higher than that of the Sun, and the stellar brightness modulation are much more regular over rotational timescales.

We find that the α factor increases with rotation rate, indicating that faculae become less important on stars rotating faster than the Sun. We have also found that the facular contribution to solar brightness variability is larger than its contribution to brightness variability in a sample of stars having near-solar rotation periods and temperatures. We attribute this to a selection effect since the rotation periods of stars with brightness-variability patterns similar to that of the Sun are rather difficult to measure via the ACF method for rotation-period determinations and, thus, there is a dearth of such stars in our sample. Consequently, our results indicate that in addition to being more active than the Sun (see also Reinhold et al. 2020; Zhang et al. 2020), the stars with near-solar effective temperatures and near-solar rotation periods determined by McQuillan et al. (2014) have different compositions of active regions (with smaller facular contributions). The GPS method for determining rotation periods could thus prove to be an important contributor to enhance the lower-than-expected number of G-type stars with near-solar rotation periods reported by van Saders et al. (2019). This method can also improve the solar-stellar comparison as in Reinhold et al. (2020). The outcome of GPS might bring a new perspective in understanding stellar activity.

While we focus in this study on applying GPS to stars with known rotation periods, in a forthcoming study, we plan to apply the GPS method to Kepler stars with previously unknown rotation periods as well as to TESS stars. This might help establish a new and more complete sample of stars having near-solar rotation periods, based upon which we can investigate whether solar variability still appears anomalously low in comparison to stars in this broader sample.

This will be of importance to the exoplanet community, since the knowledge of rotation periods will help identify radial velocity jitter from planetary signals (see Oshagh 2018; Faria et al. 2020; Hojjatpanah et al. 2020). The anticipated GPS-determined expanded database of stellar rotation periods could also bring crucial information for ongoing and upcoming surveys such as NIRPS-HARPS and ESPRESSO (see Pepe et al. 2010; Bouchy & Doyon 2018). Additionally, the precise determination of host-star rotation periods is important for recovering accurate exoplanet radii, which will be crucial for future searches of transiting Earths or super-Earths in the light curves of solar twins in the PLATO field (see Rauer et al. 2014).

Acknowledgements

We would like to thank the referee for the constructive comments which helped to improve the quality of this paper. This work was supported by the International Max-Planck Research School (IMPRS) for Solar System Science at the University of Göttingen and European Research Council under the European Union Horizon 2020 research and innovation program (grant agreement by the No. 715947). M.O. acknowledges the support of the Deutsche Forschungsgemeinschft (DFG) priority program SPP 1992, Exploring the Diversity of Extrasolar Planets (RE1664/17-1). E.M.A.G. and M.O. also acknowledge the support of the FCT/DAAD bilateral grant 2019 (DAAD ID: 57453096). Financial support was also provided by the Brain Korea 21 plus program through the National Research Foundation funded by the Ministry of Education of Korea and by the German Federal Ministry of Education and Research under project 01LG1209A. We would like to thank the International Space Science Institute, Bern, for their support of science team 446 and the resulting helpful discussions. This paper includes data collected by the Kepler mission. Funding for the Kepler mission is provided by the NASA Science Mission directorate. The data presented were obtained from the Mikulski Archive for Space Telescopes (MAST). STScI is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555.

Appendix A Additional figures

thumbnail Fig. A.1

Logarithmic visualisation of Fig. 3.

thumbnail Fig. A.2

Comparison for the 172 stars in sample B with reported rotation periods by the ACF, GLS, and GPS methods. Left: ACF versus GLS. Right: HFIP-GPS versus GLS. Scatter colour to visualise the temperature. The comparison is made in a similar range scale for a better visualisation.

thumbnail Fig. A.3

Logarithmic visualisation of bottom panel in Fig. 4.

thumbnail Fig. A.4

Variability range in ppm versus Kepler magnitude Kmag. The colour bar indicates Kepler precision.

References

  1. Aigrain, S., Llama, J., Ceillier, T., et al. 2015, MNRAS, 450, 3211 [NASA ADS] [CrossRef] [Google Scholar]
  2. Amazo-Gómez, E. M., Shapiro, A. I., Solanki, S. K., et al. 2020, A&A, 636, A69 [CrossRef] [EDP Sciences] [Google Scholar]
  3. Angus, R., Morton, T., Aigrain, S., Foreman-Mackey, D., & Rajpaul, V. 2018, MNRAS, 474, 2094 [NASA ADS] [CrossRef] [Google Scholar]
  4. Basri, G. 2018, ApJ, 865, 142 [CrossRef] [Google Scholar]
  5. Basri, G., Walkowicz, L. M., Batalha, N., et al. 2011, AJ, 141, 20 [NASA ADS] [CrossRef] [Google Scholar]
  6. Borucki, W. J., Koch, D., Basri, G., et al. 2010, Science, 327, 977 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  7. Bouchy, F., & Doyon, R. 2018, in European Planetary Science Congress, EPSC2018-1147 [Google Scholar]
  8. Buzasi, D., Lezcano, A., & Preston, H. L. 2016, J. Space Weather Space Climate, 6, A38 [CrossRef] [EDP Sciences] [Google Scholar]
  9. Cameron, R. H., Jiang, J., Schmitt, D., & Schüssler, M. 2010, ApJ, 719, 264 [NASA ADS] [CrossRef] [Google Scholar]
  10. Chapman, G. A., Cookson, A. M., & Dobias, J. J. 1997, ApJ, 482, 541 [NASA ADS] [CrossRef] [Google Scholar]
  11. Douglas, S. T., Agüeros, M. A., Covey, K. R., & Kraus, A. 2017, ApJ, 842, 83 [NASA ADS] [CrossRef] [Google Scholar]
  12. Faria, J. P., Adibekyan, V., Amazo-Gómez, E. M., et al. 2020, A&A, 635, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Foukal, P. 1998, ApJ, 500, 958 [NASA ADS] [CrossRef] [Google Scholar]
  14. Fröhlich, C., Crommelynck, D. A., Wehrli, C., et al. 1997, Sol. Phys., 175, 267 [NASA ADS] [CrossRef] [Google Scholar]
  15. García, R. A., Ceillier, T., Salabert, D., et al. 2014, A&A, 572, A34 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Hathaway, D. H. 2015, Liv. Rev. Sol. Phys., 12, 4 [Google Scholar]
  17. He, H., Wang, H., & Yun, D. 2015, ApJS, 221, 18 [NASA ADS] [CrossRef] [Google Scholar]
  18. Hojjatpanah, S., Oshagh, M., Figueira, P., et al. 2020, A&A, 639, A35 [CrossRef] [EDP Sciences] [Google Scholar]
  19. Huber, D., Silva Aguirre, V., Matthews, J. M., et al. 2014, ApJS, 211, 2 [NASA ADS] [CrossRef] [Google Scholar]
  20. Lammer, H. 2013, Origin and Evolution of Planetary Atmospheres (Springer, Berlin Heidelberg), https://doi.org/10.1007/978-3-642-32087-3 [CrossRef] [Google Scholar]
  21. Lanza, A. F., & Shkolnik, E. L. 2014, MNRAS, 443, 1451 [NASA ADS] [CrossRef] [Google Scholar]
  22. Mandal, S., Krivova, N. A., Solanki, S. K., Sinha, N., & Banerjee, D. 2020, A&A, 640, A78 [CrossRef] [EDP Sciences] [Google Scholar]
  23. Mathur, S., Salabert, D., García, R. A., & Ceillier, T. 2014, J. Space Weather Space Climate, 4, A15 [Google Scholar]
  24. McQuillan, A., Mazeh, T., & Aigrain, S. 2014, ApJS, 211, 24 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  25. Nielsen, M. B., Gizon, L., Schunker, H., & Karoff, C. 2013, A&A, 557, L10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Oshagh, M. 2018, Asteroseismology and Exoplanets: Listening to the Stars and Searching for New Worlds, 49, 239 [NASA ADS] [CrossRef] [Google Scholar]
  27. Pepe, F. A., Cristiani, S., Rebolo Lopez, R., et al. 2010, in Proc. SPIE, Vol. 7735, Ground-based and Airborne Instrumentation for Astronomy III, 77350F [Google Scholar]
  28. Pinsonneault, M. H., An, D., Molenda-Żakowicz, J., et al. 2012, ApJS, 199, 30 [NASA ADS] [CrossRef] [Google Scholar]
  29. Rauer, H., Catala, C., Aerts, C., et al. 2014, Exp. Astron., 38, 249 [NASA ADS] [CrossRef] [Google Scholar]
  30. Reinhold, T., & Gizon, L. 2015, A&A, 583, A65 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Reinhold, T., Shapiro, A. I., Solanki, S. K., et al. 2020, Science, 368, 518 [Google Scholar]
  32. Santos, A. R. G., García, R. A., Mathur, S., et al. 2019, ApJS, 244, 21 [Google Scholar]
  33. Shapiro, A. I., Solanki, S. K., Krivova, N. A., et al. 2014, A&A, 569, A38 [CrossRef] [EDP Sciences] [Google Scholar]
  34. Shapiro, A. I., Solanki, S. K., Krivova, N. A., et al. 2017, Nat. Astron., 1, 612 [NASA ADS] [CrossRef] [Google Scholar]
  35. Shapiro, A. I., Amazo-Gómez, E. M., Krivova, N. A., & Solanki, S. K. 2020, A&A, 633, A32 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Smith, J. C., Stumpe, M. C., Van Cleve, J. E., et al. 2012, PASP, 124, 1000 [Google Scholar]
  37. Solanki, S. K. 1993, Space Sci. Rev., 63, 1 [NASA ADS] [CrossRef] [Google Scholar]
  38. Solanki, S. K., Inhester, B., & Schüssler, M. 2006, Rep. Progr. Phys., 69, 563 [CrossRef] [Google Scholar]
  39. Stumpe, M. C., Smith, J. C., Catanzarite, J. H., et al. 2014, PASP, 126, 100 [Google Scholar]
  40. Thompson, S. E., Caldwell, D. A., Jenkins, J. M., et al. 2016, Kepler Data Release 25 Notes, Kepler Science Document [Google Scholar]
  41. Torrence, C., & Compo, G. P. 1998, Bull. Am. Meteorol. Soc., 79, 61 [Google Scholar]
  42. Van Cleve, J. E., & Caldwell, D. A. 2016, Kepler Instrument Handbook, Tech. rep. [Google Scholar]
  43. van Saders, J. L., Pinsonneault, M. H., & Barbieri, M. 2019, ApJ, 872, 128 [NASA ADS] [CrossRef] [Google Scholar]
  44. Walkowicz, L. M., & Basri, G. S. 2013, MNRAS, 436, 1883 [NASA ADS] [CrossRef] [Google Scholar]
  45. Witzke, V., Reinhold, T., Shapiro, A. I., Krivova, N. A., & Solanki, S. K. 2020, A&A, 634, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Zhang, J., Shapiro, A. I., Bi, S., et al. 2020, ApJ, 894, L11 [CrossRef] [Google Scholar]

1

A logarithmic visualization of the relation α factor versus Prot is available in the appendix, see Fig. A.3.

All Tables

Table 1

Stellar parameters for samples A and B.

Table 2

GPS outcome values.

Table 3

Mean α-values in sample B per bin.

All Figures

thumbnail Fig. 1

Panel I: temperature–rotation diagram for a sample of 34030 stars (coloured circles indicating the variability range) with rotation periods determined by McQuillan et al. (2014), along with 55501 stars where the authors found a period but evaluated it as not significant (beige dots; see panel III for a better visualization). Panel II: sample of 34030 stars with determined rotation periods in grey and 55501 stars with non-significant rotation period determinations in beige. For panels II and III, only the stars from sample A (panel II) and sample B (panel III) are shown in colour; see Table 1 for the properties of samples A and B. Panel II illustrates the stellar sample A, selected by near-solar rotation period and temperatures from 5000 to 6000 K. Panel III illustrates stellar sample B, which contains stars with near-solar effective temperature and a broad range of rotation periods. The Sun is represented by the solar symbol ⊙.

In the text
thumbnail Fig. 2

Histograms for variability ranges from samples A (violet outline) and B (blue rectangles). Vertical green dashed lines represent solar variability calculated for six activity seasons using 21 yr of VIRGO TSI data from Paper II (see text for more details).

In the text
thumbnail Fig. 3

Position of the high-frequency inflection point (PHFIP) plotted against rotation period. Only sample B is shown here. Rotation periods are taken from Reinhold & Gizon (2015) (left panel) and McQuillan et al. (2014) (right panel). Colours represent the stellar effective temperature, Teff. The Sun is represented by the solar symbol ⊙. Dashed lines in both panels indicate a linear fit constrained to go through the origin of the coordinate system. A logarithmic visualization is available in the appendix, see Fig. A.1.

In the text
thumbnail Fig. 4

Top panel: α factor versus Teff for sample A exhibits consistency across a broad temperature range. The black line corresponds to the linear fit to values with an uncertainty within 2-σ of the mean of the distribution as shown by coloured squares (grey squares lie outside the 2-σ of the distribution). The histograms to the right side of the panel display the distribution of α values for two effective temperature regimes, using Teff = 5500 K as a threshold. Bottom panel: α factor versus Prot from McQuillan et al. (2014) for sample B shows a slight decrease in α with rotation period. The coloured segments indicate the mean of α for the different Prot ranges as indicated in Table 3. The histograms to the right of the panel indicate the distribution of α values for two rotation period regimens, using 21 d as a threshold. For both panels, the error bars represent 2-σ uncertainties of the α values over all Kepler quarters available for each star. The grey squares lie outside of a 2-σ of the distribution. The dashed red and green horizontal lines represent the α factor values in the extreme cases with all variability being due to spots (α = 0.21) and all due to faculae (α = 0.14), respectively1.

In the text
thumbnail Fig. 5

α factor versus photometric precision for sample A (left panel) and sample B (right panel) for Kepler observations. Those records for single long observation of isolated stars observed in an uncrowded pixel, i.e. light curves of resolved targets without contamination of additional sources. Extreme-case limits for spot- and faculae-dominated stars are shown as horizontal dashed linesin red and green, respectively. The individual data points are coloured according to the detected variability range for that particular star. This is consistent with the preceding figure, Fig. 4. The error bars represent 2-σ uncertainties of the α values distribution over all Kepler quarters available per star. Gray rhomboids represent data points that lie more than 2-σ from the centre of the distribution.

In the text
thumbnail Fig. A.1

Logarithmic visualisation of Fig. 3.

In the text
thumbnail Fig. A.2

Comparison for the 172 stars in sample B with reported rotation periods by the ACF, GLS, and GPS methods. Left: ACF versus GLS. Right: HFIP-GPS versus GLS. Scatter colour to visualise the temperature. The comparison is made in a similar range scale for a better visualisation.

In the text
thumbnail Fig. A.3

Logarithmic visualisation of bottom panel in Fig. 4.

In the text
thumbnail Fig. A.4

Variability range in ppm versus Kepler magnitude Kmag. The colour bar indicates Kepler precision.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.