Inflection point in the power spectrum of stellar brightness variations III: Faculae vs. Spot dominance on stars with known rotation periods

Stellar rotation periods can be determined by observing brightness variations caused by active magnetic regions transiting visible stellar disk as the star rotates. The successful stellar photometric surveys stemming from the Kepler and TESS observations led to the determination of rotation periods in tens of thousands of young and active stars. However, there is still a lack of information about rotation periods of older and less active stars, like the Sun. The irregular temporal profiles of light curves caused by the decay times of active regions, which are comparable to or even shorter than stellar rotation periods, combine with the random emergence of active regions to make period determination for such stars very difficult. We tested the performance of the new method for the determination of stellar rotation periods against stars with previously determined rotation periods. The method is based on calculating the gradient of the power spectrum (GPS) and identifying the position of the inflection point (i.e. point with the highest gradient). The GPS method is specifically aimed at determining rotation periods of low activity stars like the Sun. We applied the GPS method to 1047 Sun-like stars observed by the Kepler telescope. We show that the GPS method returns precise values of stellar rotation periods. Furthermore, it allows us to constrain the ratio between facular and spot areas of active regions at the moment of their emergence. We show that relative facular area decreases with stellar rotation rate. Our results suggest that the GPS method can be successfully applied to retrieve periods of stars with both regular and non-regular light curves.


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).
Full Table 2 is only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc. u-strasbg.fr/viz-bin/cat/J/A+A/642/A225 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 lowerthan-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 A&A 642, A225 (2020)  Huber et al. (2014). (2) Variability range (Var) and rotation periods (P rot ) are taken from McQuillan et al. (2014). We take the solar synodic (Sy) and sidereal (Sid) Carrington rotation period values as a reference.
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, P HFIP , 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, P rot , via the simple scaling relation: 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, S fac /S spot . 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-tospot area ratio leads to a certain degree of uncertainty (up to 25%) when determining stellar rotation periods since the value of S fac /S spot 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 spotdominated 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.

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 T eff 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 (Q 0 -Q 17 ) due to the telescope re-orienting itself every 90 days. The Kepler observing quarters resulted in Q 0 of 10 days, and Q 1 of 33 days for the commissioning phase and, segmented 90-day LCs for Q 2 to Q 16 (see public data release 25, Thompson et al. 2016;Van Cleve & Caldwell 2016). The second month of Q 17 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 T eff ∈ (5000-6000) K. These selection criteria yielded a sample consisting of 686 stars. Among this sample, 282 stars also have   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 . 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 15 and 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 Q 1 -Q 14 processed with the pre-search conditioning and Bayesian maximum a posteriori approach (PDC-MAP, see Smith et al. 2012). For quarters Q 15 −Q 17 , 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 fundamental parameters 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 (see e.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 A225, page 3 of 10 A&A 642, A225 (2020) (1) ( 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 photometric variability are typical for stars having near-solar fundamental parameters but unknown rotation periods.

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 Q 1 -Q 17 . We determined the period corresponding to the high-frequency inflection point, P HFIP (Qn ) , per quarter and calculate the mean value for P HFIP over all 17 quarters for each star. This allows us to obtain a unique representative value of P HFIP per star. The uncertainty is calculated using 2-σ of the distribution of the obtained P HFIP values. Finally, we used the P HFIP to calculate the stellar rotation period, P rot (see Table 2 and on-line reference for a compilation of GPS outputs and comparison with GLS and ACF reference values). In Fig. 3 The scatter around the linear fits has multiple sources. First, the calibration coefficient between rotation period and inflection point, α = P HFIP /P rot , depends on the relative roles that bright faculae and dark spots play in generating stellar brightness variations. According to the model presented in Paper I, these roles are regulated by the ratio between facular and spot areas of active regions at the time of emergence, S fac /S spot (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 S fac /S spot 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)  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 S fac /S spot 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 nearsolar 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 independent of 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 S fac /S spot ratio in terms of Paper I) with rotation period. The increase of the α factor with rotation rate implies that the S fac /S spot ratio (and, consequently, the contribution of faculae to the rotational variability) is lower in faster A225, page 5 of 10 A&A 642, A225 (2020) Fig. 4. Top panel: α factor versus T eff 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 T eff = 5500 K as a threshold. Bottom panel: α factor versus P rot 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 P rot 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), respectively 1 .  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 S fac /S spot ratio decreases with the size of active regions and their spot components. An 1 A logarithmic visualization of the relation α factor versus P rot is available in the appendix, see Fig. A.3. 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 S fac /S spot 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 A225, page 6 of 10 E. M. Amazo-Gómez et al.: GPS III. Faculae vs. Spot dominance on stars with known rotation periods Extreme-case limits for spot-and faculae-dominated stars are shown as horizontal dashed lines in 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. 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 S fac /S spot 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 S fac /S spot ratios are smaller and, consequently, their α factors are larger.

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 brightnessvariability 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 nearsolar 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).