This is our reply to the comment of T. Metcalfe and J. van Saders on the Science report "The Sun is less active than other solar-like stars" by T. Reinhold, A. I. Shapiro, S. K. Solanki, B. T. Montet, N. A. Krivova, R. H. Cameron, E. M. Amazo-Gomez. We hope that both the comment and our reply lead to fruitful discussions which of the two presented scenarios is more likely.

Obtaining measurements of chromospheric and photometric activity of stars with near-solar fundamental parameters and rotation periods is important for a better understanding of solar-stellar connection. We select a sample of 2603 stars with near-solar fundamental parameters from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST)-Kepler field and use LAMOST spectra to measure their chromospheric activity and Kepler light curves to measure their photospheric activity (i.e., the amplitude of the photometric variability). While the rotation periods of 1556 of these stars could not be measured due to the low amplitude of the photometric variability and highly irregular temporal profile of light curves, 254 stars were further identified as having near-solar rotation periods. We show that stars with near-solar rotation periods have chromospheric activities that are systematically higher than stars with undetected rotation periods. Furthermore, while the solar level of photospheric and chromospheric activity appears to be typical for stars with undetected rotation periods, the Sun appears to be less active than most stars with near-solar rotation periods (both in terms of photospheric and chromospheric activity).

Magnetic activity of the Sun and other stars causes their brightness to vary. We investigate how typical the Sun's variability is compared to other solar-like stars, i.e. those with near-solar effective temperatures and rotation periods. By combining four years of photometric observations from the Kepler space telescope with astrometric data from the Gaia spacecraft, we measure photometric variabilities of 369 solar-like stars. Most of the solar-like stars with well-determined rotation periods show higher variability than the Sun and are therefore considerably more active. These stars appear nearly identical to the Sun, except for their higher variability. Their existence raises the question of whether the Sun can also experience epochs of such high variability.

Young and active stars generally have regular, almost sinusoidal, patterns of variability attributed to their rotation, while the majority of older and less active stars, including the Sun, have more complex and non-regular light-curves which do not have clear rotational-modulation signals. Consequently, the rotation periods have been successfully determined only for a small fraction of the Sun-like stars observed by transit-based planet-hunting missions, such as CoRoT, Kepler, and TESS. This suggests that only a small fraction of such systems have been properly identified as solar-like analogs. We apply a new method for determining rotation periods of low-activity stars, like the Sun. The method is based on calculating the gradient of the power spectrum (GPS) of stellar brightness variations and identifying a tell-tale inflection point in the spectrum. The rotation frequency is then proportional to the frequency of that inflection point. In this paper test this GPS method against available photometric records of the Sun. We apply GPS, autocorrelation functions, Lomb-Scargle periodograms, and wavelet analyses to the total solar irradiance (TSI) time series obtained from the Total Irradiance Monitor (TIM) on the Solar Radiation and Climate Experiment (SORCE) and the Variability of solar IRradiance and Gravity Oscillations (VIRGO) experiment on the SOlar and Heliospheric Observatory (SoHO) missions. We analyse the performance of all methods at various levels of solar activity. We show that the GPS method returns accurate values of solar rotation independently of the level of solar activity. In particular, it performs well during periods of high solar activity, when TSI variability displays an irregular pattern and other methods fail. Our results suggest that the GPS method can successfully determine the rotational periods of stars with both regular and non-regular light-curves.

Rotation period measurements of stars observed with the Kepler mission have revealed a lack of stars at intermediate rotation periods, accompanied by a decrease of photometric variability. Whether this so-called dearth region is a peculiarity of stars in the Kepler field, or reflects a general manifestation of stellar magnetic activity, is still under debate. Our goal is to measure stellar rotation periods and photometric variabilities for tens of thousands of K2 stars, located in different fields along the ecliptic plane, to shed light on the relation between stellar rotation and photometric variability. We use Lomb-Scargle periodograms, auto-correlation and wavelet functions to determine consistent rotation periods. Stellar brightness variability is assessed by computing the variability range from the light curve. We further apply Gaussian mixture models to search for bimodality in the rotation period distribution. Combining measurements from all K2 campaigns, we detect rotation periods in 29,860 stars. For effective temperatures below 6000K, the variability range shows a local minimum at different periods, consistent with an isochrone age of 750 Myr. Additionally, the K2 rotation period distribution shows evidence for bimodality, although the dearth region is less pronounced compared to the Kepler field. The period at the dip of the bimodal distribution shows good agreement with the period at the local variability minimum. We conclude that the period bimodality is present in different fields of the sky, and is hence a general manifestation of stellar magnetic activity. The reduced variability in the dearth region is interpreted as a cancelation between dark spots and bright faculae. Our results strongly advocate that the role of faculae has been underestimated so far, suggesting a more complex dependence of the brightness variability on the rotation period.

Accurate determination of stellar rotation periods is important for estimating stellar ages as well as for understanding stellar activity and evolution. While for about thirty thousand stars in the Kepler field rotation periods can be determined, there are over hundred thousand stars, especially with low photometric variability and irregular pattern of variations, for which rotational periods are unknown. Here, we investigate the effect of metallicity on the detectability of rotation periods. This is done by synthesising light curves of hypothetical stars, which are identical to our Sun, with the exception of the metallicity. These light curves are then used as an input to the period determination algorithms. We find that the success rate for recovering the rotation signal has a minimum close to the solar metallicity value. This can be explained by the compensation effect of facular and spot contributions. In addition, selecting solar-like stars with near-solar effective temperature, near solar photometric variability, and with metallicity between M/H = -0.35 and M/H = 0.35 from the Kepler sample, we analyse the fraction of stars for which rotational periods have been detected as a function of metallicity. In agreement with our theoretical estimate we found a local minimum for the detection fraction close to the solar metallicity. We further report rotation periods of 87~solar-like Kepler stars for the first time.

The study of stellar activity cycles is crucial to understand the underlying dynamo and how it causes activity signatures such as dark spots and bright faculae. We study the appearance of activity signatures in contemporaneous photometric and chromospheric time series. Lomb-Scargle periodograms are used to search for cycle periods present in both time series. To emphasize the signature of the activity cycle we account for rotation-induced scatter in both data sets by fitting a quasi-periodic Gaussian process model to each observing season. After subtracting the rotational variability, cycle amplitudes and the phase difference between the two time series are obtained by fitting both time series simultaneously using the same cycle period. We find cycle periods in 27 of the 30 stars in our sample. The phase difference between the two time series reveals that the variability in fast rotating active stars is usually in anti-phase, while the variability of slowly rotating inactive stars is in phase. The photometric cycle amplitudes are on average six times larger for the active stars. The phase and amplitude information demonstrates that active stars are dominated by dark spots, whereas less active stars are dominated by bright faculae. We find the transition from spot to faculae domination at the Vaughan-Preston gap, and around a Rossby number equal to one. We conclude that faculae are the dominant ingredient of stellar activity cycles at ages >2.55 Gyr. The data further suggest that the Vaughan-Preston gap can not explain the previously detected dearth of Kepler rotation periods between 15-25 days. Nevertheless, our results led us to propose an explanation for the rotation period dearth to be due to the non-detection of periodicity caused by the cancellation of dark spots and bright faculae at 800 Myr.

In recent years it has been claimed that the length of stellar activity cycles is determined by the stellar rotation rate. It is observed that the cycle period increases with rotation period along the so-called active and inactive sequences. In this picture the Sun occupies a solitary position in between the two sequences. Our goal is to measure cyclic variations of the stellar light curve amplitude and the rotation period using four years of Kepler data. Periodic changes of the light curve amplitude or the stellar rotation period are associated with an underlying activity cycle. Using the McQuillan et al. 2014 sample we compute the rotation period and the variability amplitude for each Kepler quarter and search for periodic variations of both time series. To test for periodicity in each stellar time series we consider Lomb-Scargle periodograms and use a selection based on a False Alarm Probability (FAP). We detect amplitude periodicities in 3203 stars between 0.5-6 years covering rotation periods between 1-40 days. Given our sample size of 23,601 stars and our selection criteria that the FAP is less than 5%, this number is almost three times higher than that expected from pure noise. We do not detect periodicities in the rotation period beyond those expected from noise. Our measurements reveal that the cycle period shows a weak dependence on rotation rate, slightly increasing for longer rotation period. We further show that the shape of the variability deviates from a pure sine curve, consistent with observations of the solar cycle. Our measurements do not support the existence of distinct sequences in the P_rot-P_cyc plane, although there is some evidence for the inactive sequence for rotation periods between 5-25 days. Unfortunately, the total observing time is too short to draw sound conclusions on activity cycles with similar length as the solar cycle.

The high-precision photometry from the CoRoT and Kepler satellites has led to measurements of surface rotation periods for tens of thousands of stars. Our main goal is to derive ages of thousands of field stars using consistent rotation period measurements in different gyrochronology relations. Multiple rotation periods are interpreted as surface differential rotation (DR). We re-analyze the sample of 24,124 Kepler stars from Reinhold et al. (2013) using different approaches based on the Lomb-Scargle periodogram. Each quarter (Q1-Q14) is treated individually using a prewhitening approach. Additionally, the full time series, and different segments thereof are analyzed. For more than 18,500 stars our results are consistent with the rotation periods from McQuillan et al. (2014). Thereof, more than 12,300 stars show multiple significant peaks, which we interpret as DR. Gyrochronology ages between 100 Myr and 10 Gyr were derived for more than 17,000 stars using different gyrochronology relations. We find a bimodal age distribution for Teff between 3200-4700 K. The derived ages reveal an empirical activity-age relation using photometric variability as stellar activity proxy. Additionally, we found 1079 stars with extremely stable (mostly short) periods. Half of these periods may be associated with rotation stabilized by non-eclipsing companions, the other half might be due to pulsations. The derived gyrochronology ages are well constrained since more than 93.0 % of the stars seem to be younger than the Sun where calibration is most reliable. Explaining the bimodality in the age distribution is challenging, and limits accurate stellar age predictions. The existence of cool stars with almost constant rotation period over more than three years of observation might be explained by synchronization with stellar companions, or a dynamo mechanism keeping the spot configurations extremely stable.

We present the results of a blind exercise to test the recoverability of stellar rotation and differential rotation in Kepler light curves. The simulated light curves lasted 1000 days and included activity cycles, Sun-like butterfly patterns, differential rotation and spot evolution. The range of rotation periods, activity levels and spot lifetime were chosen to be representative of the Kepler data of solar like stars. Of the 1000 simulated light curves, 770 were injected into actual quiescent Kepler light curves to simulate Kepler noise. The test also included five 1000-day segments of the Sun's total irradiance variations at different points in the Sun's activity cycle. Five teams took part in the blind exercise, plus two teams who participated after the content of the light curves had been released. The methods used included Lomb-Scargle periodograms and variants thereof, auto-correlation function, and wavelet-based analyses, plus spot modelling to search for differential rotation. The results show that the `overall' period is well recovered for stars exhibiting low and moderate activity levels. Most teams reported values within 10% of the true value in 70% of the cases. There was, however, little correlation between the reported and simulated values of the differential rotation shear, suggesting that differential rotation studies based on full-disk light curves alone need to be treated with caution, at least for solar-type stars. The simulated light curves and associated parameters are available online for the community to test their own methods.

Surface differential rotation (DR) is one major ingredient of the magnetic field generation process in the Sun and likely in other stars. The term solar-like differential rotation describes the observation that solar equatorial regions rotate faster than polar ones. The opposite effect of polar regions rotating faster than equatorial ones (termed as antisolar DR) has only been observed in a few stars, although there is evidence from theoretical dynamo models. We present a new method to detect the sign of DR (i.e. solar-like or antisolar DR) by analyzing long-term high-precision light curves with the Lomb-Scargle periodogram.We compute the Lomb-Scargle periodogram and identify a set of significant periods $P_k$, which we associate with active regions located at different latitudes on the the stellar surface. If detectable, the first harmonics ($P_k'$) of these periods were identified to compute their peak-height-ratios $r_k:=h(P_k')/h(P_k)$. Spots rotating at lower latitudes generate less sine-shaped light curves, which requires additional power in the harmonics, and results in larger ratios $r_k$. Comparing different ratios $r_k$ and the associated periods $P_k$ yields information about the spot latitudes, and reveals the sign of DR. We tested our method on different sets of synthetic light curves all exhibiting solar-like DR. The number of cases where our method detects antisolar DR is the false-positive rate of our method. Depending on the set of light curves, the noise level, the required minimum peak separation, and the presence or absence of spot evolution, our method fails to detect the correct sign in at most 20%. We applied our method to 50 Kepler G stars and found 21-34 stars with solar-like DR and 5-10 stars with antisolar DR, depending on the minimum peak separation.

We present rotation periods for thousands of active stars in the Kepler field derived from Q3 data. In most cases a second period close to the rotation period was detected, which we interpreted as surface differential rotation (DR). Active stars were selected from the whole sample using the range of the variability amplitude. To detect different periods in the light curves we used the Lomb-Scargle periodogram in a pre-whitening approach to achieve parameters for a global sine fit. The most dominant periods from the fit were ascribed to different surface rotation periods, but spot evolution could also play a role. Due to the large number of stars the period errors were estimated in a statistical way. We thus cannot exclude the existence of false positives among our periods. In our sample of 40.661 active stars we found 24.124 rotation periods $P_1$ between 0.5-45 days. The distribution of stars with 0.5 < B-V < 1.0 and ages derived from angular momentum evolution that are younger than 300 Myr is consistent with a constant star-formation rate. A second period $P_2$ within $\pm30$% of the rotation period $P_1$ was found in 18.619 stars (77.2%). Attributing these two periods to DR we found that the relative shear $\alpha=\Delta\Omega/\Omega$ increases with rotation period, and slightly decreases with effective temperature. The absolute shear $\Delta\Omega$ slightly increases between $T_{eff}=3500-6000$ K. Above 6000 K $\Delta\Omega$ shows much larger scatter. We found weak dependence of $\Delta\Omega$ on rotation period. Latitudinal differential rotation measured for the first time in more than 18.000 stars provides a comprehensive picture of stellar surface shear, consistent with major predictions from mean-field theory. To what extent our observations are prone to false positives and selection bias is not fully explored, and needs to be addressed using more Kepler data.

Co-rotating spots at different latitudes on the stellar surface generate periodic photometric variability and can be useful proxies to detect Differential Rotation (DR). DR is a major ingredient of the solar dynamo but observations of stellar DR are rather sparse. In view of the Kepler space telescope we are interested in the detection of DR using photometric information of the star, and to develop a fast method to determine stellar DR from photometric data. We ran a large Monte-Carlo simulation of differentially rotating spotted stars with very different properties to investigate the detectability of DR. For different noise levels the resulting light curves are prewhitened using Lomb-Scargle periodograms to derive parameters for a global sine fit to detect periodicities. We show under what conditions DR can successfully be detected from photometric data, and in which cases the light curve provides insufficient or even misleading information on the stellar rotation law. In our simulations, the most significant period P1_{out} could be detected in 96.2% of all light curves. Detection of a second period close to P1_{out} is the signature of DR in our model. For the noise-free case, in 64.2% of all stars such a period was found. Calculating the measured latitudinal shear of two distinct spots \alpha_{out}, and comparing it to the known original spot rotation rates shows that the real value is on average 3.2% lower. Comparing the total equator-to-pole shear $\alpha$ to $\alpha_{out}$ we find that $\alpha$ is underestimated by 8.8%, esp. the detection of DR for stars with $\alpha$ < 6% is challenging. Finally, we apply our method to four differentially rotating Kepler stars and find close agreement with results from detailed modeling. Our method is capable of measuring stellar rotation periods and detecting DR with relatively high accuracy and is suitable for large data sets.

The Kepler space telescope monitors over 156.000 stars with an unprecedented photometric precision. We are interested in stellar rotational periods which we find using Lomb-Scargle periodograms. This work focuses on the 306 exoplanet candidate host stars released on June 15, 2010. We present statistics on how many of them show periodic photometric variability, providing preliminary periods and estimates of stellar activity. In the future, our work will focus on spot evolution and differential rotation.