Issue 
A&A
Volume 576, April 2015



Article Number  A15  
Number of page(s)  8  
Section  Stellar atmospheres  
DOI  https://doi.org/10.1051/00046361/201425337  
Published online  13 March 2015 
Discriminating solar and antisolar differential rotation in highprecision light curves^{⋆}
^{1} Institut für Astrophysik, Universität Göttingen, FriedrichHundPlatz 1, 37077 Göttingen, Germany
email: reinhold@astro.physik.unigoettingen.de
^{2} LeibnizInstitut für Astrophysik Potsdam (AIP), An der Sternwarte 16, 14482 Potsdam, Germany
Received: 14 November 2014
Accepted: 27 January 2015
Context. Surface differential rotation (DR) is one major ingredient of the magnetic field generation process in the Sun and likely in other stars. The term solarlike 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.
Aims. We present a new method of detecting the sign of DR (i.e., solarlike or antisolar DR) by analyzing longterm highprecision light curves with the LombScargle periodogram.
Methods. We compute the LombScargle periodogram and identify a set of significant periods P_{k}, which we associate with active regions located at different latitudes on the stellar surface. If detectable, the first harmonics () of these periods were identified to compute their peakheightratios r_{k}:= h(P'_{k})/h(P_{k}) . Spots rotating at lower latitudes generate less sineshaped 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.
Results. We tested our method on different sets of synthetic light curves all exhibiting solarlike DR. The number of cases where our method detects antisolar DR is the falsepositive 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 solarlike DR and 5–10 stars with antisolar DR, depending on the minimum peak separation.
Conclusions. The method is able to determine the sign of DR in a statistical way with a low falsepositive rate. Applying our method to real data might suggest that – within the uncertainties – antisolar DR was detected in 5–10 Kepler stars.
Key words: stars: rotation / starspots
Table 3 is available in electronic form at http://www.aanda.org
© ESO, 2015
1. Introduction
The differential rotation (DR) of the solar surface (and its interior) is a crucial part of the generation of the Sun’s magnetic field and the solar dynamo. The differential rotation winds up magnetic field lines, which can be considered as frozen in the plasma, transferring an initial poloidal field into a toroidal one. Speaking of solar or solarlike differential rotation in the following we always think of latitudinaldependent surface rotation rates, with equatorial regions rotating faster than polar ones. This latitudinal shear is usually described by the equation (1)where Ω_{eq} denotes the equatorial angular velocity, θ the latitude, and α a dimensionless parameter describing the amount of surface shear. The Sun has a relative shear of α_{⊙} = 0.2. According to Eq. (1) DR in stars other than the Sun is usually termed solarlike when α> 0, whereas the opposite effect, α< 0, is called antisolar; the case α = 0 describes rigid (surface) rotation. For other purposes DR might also be expressed in terms of the absolute surface shear ΔΩ = Ω_{eq} − Ω_{pole} = αΩ_{eq}.
Measuring DR is challenging because it acts as a subtle perturbation of the mean rotation rate. For stars other than the Sun different techniques have been used. Reiners & Schmitt (2002) applied the Fourier transform method to Doppler broadened spectral lines. Doppler imaging (Donati & Collier Cameron 1997) traces active regions located at different latitudes to infer the absolute surface shear ΔΩ. Furthermore, DR can also be measured from asteroseismology (Gizon & Solanki 2004; Lund et al. 2014). Classical spot models also supply a measure on DR, with recent studies being, e.g., Lanza et al. (2014); Bonanno et al. (2014). Using the shortterm Fouriertransform Vida et al. (2014) provided measurements of activity cycles of fast rotating Kepler stars attributing the observed change in rotation period to DR. Another approach was presented by Reinhold et al. (2013) that measured different periods in highprecision light curves. This method only yields absolute values of the shear, but provides no information of the sign because the latitudes of the active regions are usually unknown.
Figure 1 illustrates that it is difficult to tell from the light curve whether the star exhibits solar or antisolar DR. Both light curves were drawn from the same parameters and show an almost identical beating pattern, with the exception that α is positive (negative) in the upper (lower) panel, respectively. Thus, visual inspection of the light curve does not immediately reveal the underlying rotational behavior.
Fig. 1 Upper panel: simulated light curve of a star with two spots and an equatortopole shear α = 0.085 (solarlike DR) according to Eq. (1). Lower panel: same parameters as above but with a negative value of α = −0.085 instead (antisolar DR). 
Recently, 3D numerical simulations have tabled the question of which physical conditions lead to solarlike or antisolar differential rotation, and which parameters determine the transition region (Guerrero et al. 2013; Gastine et al. 2014; Käpylä et al. 2014; Karak et al. 2015). A theoretical explanation for the existence of antisolar DR is provided by Kitchatinov & Rüdiger (2004) that favors the explanation of fast meridional motion. Nevertheless, observations of antisolar DR are sparse and have only been claimed in some K giants so far. Strassmeier et al. (2003) found strong signatures of antisolar DR on the K2giant HD 31993 supplying the value α = −0.125 ± 0.05. Weber et al. (2005) found marginal evidence for antisolar DR in several K giants. Kővári et al. (2007) detected antisolar DR of the K1 giant σ Geminorum providing a relative shear of α ≈ −0.022 ± 0.006. Later studies of the same object approved antisolar DR that were almost twice as strong as the earlier value (Kővári et al. 2015). Vida et al. (2007) observed antisolar DR in the K0 giant UZ libræ using crosscorrelation of consecutive Doppler images obtained in the years 1998 and 2000, yielding α ≈ −0.03 and α ≈ −0.004, respectively. Dikpati & Cally (2011) analyzed the stability of toroidal magnetic fields and differential rotation and found that kGfields (i.e., activity) are only compatible with small shear if the DR is antisolar. The finding agrees with the small values observed. Furthermore, Ammlervon Eiff & Reiners (2012) applied the Fourier transform method to a sample of more than 110 A–F field stars. For most of the stars in their sample these authors either found no signature of DR or very small values consistent with rigid rotation. For ten stars significant DR (0.10 <α< 0.50) was detected. Five stars in their sample showed signatures that can be explained by antisolar DR. Nevertheless, the authors interpreted the measurements as rigid rotators with cool polar caps.
The goal of this article is to present a new method of measuring the sign of DR from highprecision light curves. We compute their LombScargle periodograms and show how different periods can be associated with higher or lower latitudes using overtone amplitudes. No explicit spot modeling is needed. We explain the method in detail, present results for simulated light curves, and discuss the problems when applying this technique to Kepler data. The paper is organized as follows; in Sect. 2 we motivate the parameter space of the simulated light curves. Section 3 contains a detailed description of the method. The application of the method to simulated and real data is presented in Sect. 4. The last section is dedicated to a detailed discussion of the results.
2. Simulations
In order to test our method we simulated different sets of differentially rotating stars. Each set consists of 100 light curves and was generated by our code MODSTAR (Reinhold & Reiners 2013). The code creates a spherical surface of a rotating spotted star. Each surface element has a fixed intensity value. Spots are modeled as circular areas with a spottophotosphere contrast of 0.67, which roughly corresponds to the solar penumbratophotosphere contrast. For simplicity, umbra and penumbra structures of star spots as well as noncircular spot shapes were neglected.
The crucial parameters to describe latitudinal shear are the equatorial rotation period P_{eq} = 2π/ Ω_{eq}, the equatortopole shear α , and the spot latitudes (see Eq. (1)). Other parameters of interest are the number of spots and the stellar inclination angle i. All light curves were created using the following parameters, which were chosen at random for each light curve. The stellar inclination ranges from 30°<i< 80°, the equatorial rotation period from 1 <P_{eq}< 20 days, and the equatortopole shear covers 0.05 <α< 0.30. We emphasize that all simulated stars exhibit a positive α, and thus solarlike DR!
The first model set consists of 100 stars with only two spots at fixed latitudes θ_{1} = 10° and θ_{2} = 60°. The second set contains 100 model stars with 10 spots located at random latitudes between −30°<θ< 90°. The number of spots and their latitudes is the only difference between the two sets. The spots can cover the full longitude range −180°<ϕ< 180°, with radii between 4°<r< 6°. For the 10spot models these radii correspond to a spot coverage of ~1−3%. Additionally, spots are allowed to overlap without further contrast reduction. Each light curve covers a total length of 1500 days, with a sampling rate of 4 points/day.
The spots in our simple models have infinite lifetimes, i.e., their radii do not change with time. Spot evolution can also be implemented in the models by increasing or decreasing the spot radii. In Sect. 4 we test the robustness of our method under various perturbations, especially probing the effect of spot evolution. Since there are many ways of realizing spot evolution, we chose a simple approach by changing the spot radii when creating the light curves. At the beginning of the simulation the program chooses at random whether all spots should be increased or decreased. Over a course of 30 days, all spot radii are then successively changed by 0.4° in total. Every 30 days, the program chooses again whether the spots should be increased or decreased. Owing to the method it happens that spot radii drop below 0°. In such cases the spots are forced to increase again during the next 30 days period. Both the 30day interval as well as the amount of 0.4° growth or decay in radii was chosen such that the light curves visually exhibit similar shape to what can be found in Kepler data.
The first model set only contains stars with two spots on their surface. These sets might not be very realistic, but are sufficient to probe the signature of higher and lower latitude spots in the light curve. This set is ideal to present the power of our method, and to focus on the primary goal of detecting the correct sign of α. The other parameters were chosen to be comparable with Kepler data. We chose an upper limit of 20 days of the equatorial rotation period because many Kepler stars exhibit very stable rotation periods in this time regime. Longer periods usually become unstable owing to spot evolution, stronger differential rotation, or instrumental effects. From previous simulations we know that DR is hard to detect in stars with small equatortopole shear (Reinhold & Reiners 2013), so a lower limit of α> 0.05 was imposed. The inclination is limited to i> 30° because our method tends to fail for lower values. Furthermore, the total length of the light curves in comparable to Kepler’s observing time span.
3. Method
We compute the generalized LombScargle periodogram (Zechmeister & Kürster 2009) of the full light curves, identify the 30 highest peaks, and sort them by their peak heights^{1} in descending order. The highestpeak period is termed P_{1} and represents the most significant period in the data. To measure the amount of latitudinal DR we search for periods adjacent to P_{1} in the range of P_{1} ± 30% with peak heights greater than 50. This value represents a lower limit of the significance level of the period in the data, and was chosen after byeye inspection of several periodograms. The second highest period satisfying the above conditions is termed P_{2}, the third highest peak P_{3}, and so on.
Ideally these periods belong to certain spot rotation periods, i.e., to spots rotating at certain latitudes. For the 2spot model we found that this assignment is true for almost all cases (see Sect. 4.1). Unfortunately, for the 10spot models (Sect. 4.2) this association is no longer valid because not every spot is resolved by an individual peak. Furthermore, in real data the number of spots on the surface is unknown, and one has to assign the peaks to spots (or spot groups) on the surface to get a measure of the latitudinal DR. Naturally, one would use the highestpeak periods P_{1} and P_{2} to guarantee that these are really present in the data. However, no matter which peaks are used to get a measure on DR, they only provide an absolute value of the surface shear because there is no information about the associated spot latitudes. In the following we show that latitudinal information can be achieved by taking into account the first harmonics (half periods) of the selected periods.
The light curve signature of a single star spot rotating in and out of view depends on its radius, the intensity contrast with its surrounding, the stellar inclination angle, and the spot latitude. In our model we assume circular spot shape and neglect any temperature gradient within the spot. Our method is based on the outcome of the simulations that – for a wide parameter range – lower latitude spots generate less sineshaped light curves compared to higher latitude spots. As a consequence, additional harmonics of the primary period are required for a proper fit. Since the goodness of a sine fit is equivalent to the peak height in the periodogram, one can achieve information about the spot latitude by comparing the heights of selected peaks and their harmonics. We search the periodogram for the first harmonic of all periods P_{k},k = 1,2,... satisfying the above conditions, because in some cases the harmonics of P_{1} and P_{2} cannot be detected. If present, the harmonics are termed and should lie within 1% of the half periods P_{k}/ 2, i.e., they should satisfy . The corresponding peak heights of the periods P_{k} are termed h(P_{k}) and of the harmonics , respectively. For all pairs () satisfying the above relations we compute the peakheightratios (PHR) . Comparing different ratios now reveals information about the spot latitude. The relation r_{k}>r_{k + 1} means that the relative power of the harmonic compared to the power of the corresponding period is higher for k than for k + 1. In this case, we associate the period P_{k} with a lower latitude spot compared to the spot latitude associated with P_{k + 1}. We define (2)and (3)Since we sorted the periods according to their peak heights, we start to compare the PHR with the lowest kvalue, k_{min}, satisfying all conditions, to infer the rotation law. In most cases, k_{min} = 1 is used if is detectable. If P_{low}<P_{high} for k = k_{min} the star exhibits solarlike DR. If the opposite relation is true, antisolar DR is present. Thus, comparing the PHR and the corresponding periods yields the sign of DR!
In some cases contradictory results are returned by comparing the PHR for different kvalues, e.g., when the comparison for k = 1 reveals solarlike DR, and for k = 2 antisolar DR, or vice versa. This happens because the routine matches all peaks in the range of the rotation period to their harmonics, and associates them with spot rotation rates at different latitudes. It is not possible to associate each peak with a certain spot rotation rate, which means that in some cases the periods are associated with the wrong latitudes. To minimize this effect, we only use the relation for k = k_{min}, because k_{min} uses the two highest peaks with detectable harmonics, for comparison.
Fig. 2 Upper panel: light curve of a randomly selected star from the 10spot models without spot evolution. Lower panel: periodogram of the above light curve. The peaks of the pairs () and () are clearly visible and emphasized by black upside down and normal triangles, respectively. 
For the 10spot models we checked under which conditions a negative sign of DR is returned. To minimize the number of wrong detections we imposed some constraints on the program. For stars with too many individually resolved periods P_{k}, we restricted the search to k ≤ 5. Furthermore, the ratios r_{k} and r_{k + 1} must differ by at least 10%. The largest contribution of wrong detections results from choosing peaks very close to each other, so we forced the program to only match peaks that are separated by at least five points on the frequency grid. Relaxing this condition has a large impact on the number of detection in real data (see Sect. 4.3), and is discussed in Sect. 5. To illustrate the method Fig. 2 shows the periodogram of a randomly selected light curve from the 10spot model without spot evolution. In this model the periods and harmonics for k = 1,...,4 satisfy all conditions and are collected in Table 1. To infer the rotation law the peaks of the pairs () and () are selected and emphasized by black upside down and normal triangles, respectively. We detected the periods P_{1} = 15.15 d, P_{2} = 13.76 d, d, d, and the ratios r_{1} = 0.026 and r_{2} = 0.277. Since r_{2}>r_{1} the period P_{2} belongs to a lower latitude spot (see Eq. (3)), and the relation P_{2}<P_{1} reveals solarlike DR! The original spot rotation periods P_{spot} of this model lie in the range 13.71 d <P_{spot}< 16.47 d, showing that the periodogram detects different periods in the right range. Other parameters of the picked model are i = 60.6° and α = 0.168.
4. Results
In this section we present the results of the application of our method to the different artificial and real data sets. To quantify the results we provide the distribution of the observed relative differential rotation α_{obs}, which we define as (4)with the periods P_{high} and P_{low} from Eqs. (2) and (3). This quantity automatically provides the falsepositive rate of our method, being the number of stars with negative α_{obs} values. Negative values of α_{obs} result from associating certain periods with the wrong latitudes as discussed in Sect. 3. Again, we note that all simulated stars exhibit α> 0.
As we will show in the following, the falsepositive rate is rather small for unperturbed data. The above sets of light curves are noisefree, the spots do not evolve with time, and – compared to real data – come without gaps or instrumental effects. To test the robustness of our method we include different perturbations in the light curves. We found that spot evolution has the largest effect, so for each set we varied the spot radii with time (see Sect. 2) using the same parameters as in the unperturbed case, and considered different levels of Poisson noise.
To quantify the impact of each effect we compare the distribution of α_{obs} to the distribution of the measurable shear of the model. Therefore, we define α_{model}: = (P_{max} − P_{min}) /P_{max}, with P_{max} = max(P(θ)) and P_{min} = min(P(θ)), being the maximum and minimum periods of the model, respectively. By definition it holds α_{obs} ≤ α_{model}, because the periods P_{high} and P_{low} do not necessarily belong to the highest and lowest latitude spots, respectively. For each model set we show how well the original spot periods are recovered, directly providing the falsepositive rate of our method.
4.1. 2spot models
The distributions of α_{model} and α_{obs} are shown in Fig. 3 in solid black and dashed red, respectively. The upper panel shows the distributions without spot evolution, whereas in the lower panel spot evolution is included in the models. In both panels, the noise level increases from left to right.
For nonevolving spots we find α_{obs}> 0 in 97 of 100 cases in the noisefree models. For the three remaining stars the conditions were not fulfilled, i.e., the falsepositive rate is zero! Almost the same result was found when adding 500 ppm Poisson noise to the light curves, with 97 solarlike and one antisolar detections. Adding 1000 ppm white noise yields solarlike DR for 85 stars, and no detection for the remaining 15 light curves. Adding even higher noise seems unreasonable because the periodicity in the light curves visible to the naked eye gets corrupted. Although the main rotation periods are still detectable, their harmonics cannot be identified correctly because the large noise is interpreted as short periods in the periodograms. Such noisy light curves have not been considered when searching for DR in Kepler data (Reinhold et al. 2013).
We also considered the same models including spot evolution. For the noisefree, the 500 ppm, and the 1000 ppm models we find 94, 90, and 74 solarlike detections, respectively. Only one case of antisolar rotation was detected in the 1000 ppm model. In all cases we find good agreement between the model periods and the recovered values. Furthermore, the falsepositive rate is zero for almost all cases supporting the power of the method.
Fig. 3 Distributions of α_{model} (solid black) and α_{obs} (dashed red) of the 2spot model for different noise levels. Upper panel: no spot evolution, lower panel: spot evolution included. The dashed blue line indicates α_{obs} = 0. 
4.2. 10spot models
In Fig. 4 we show the distributions of α_{model} and α_{obs} of the 10spot models in solid black and dashed red, respectively. Again, distributions in the upper panel are without spot evolution, whereas spot evolution is included in the lower panel, and the noise level increases from left to right.
The distributions of α_{obs} look very different to the ones of the 2spot models, with α_{obs} being much lower than α_{model}, on average. The shear is underestimated in many cases because the stars exhibit more than two significant periods. The match between α_{obs} and α_{model} can be improved by picking other periods, but our primary goal was detecting the correct sign of α as accurately as possible, rather than recovering the total amount of shear. The values of α_{obs} should instead be considered as lower limits of the, in principle, measurable shear.
The falsepositive rates for the different cases are as follows. For the models without spot evolution we find 2 of 78, 3 of 78, and 5 of 77 antisolar cases for the noisefree, the 500 ppm, and the 1000 ppm models, respectively. Including spot evolution the number of wrong associations increases to 8 of 83, 8 of 82, and 9 of 78 for the noisefree, the 500 ppm, and the 1000 ppm models, respectively. This result clearly shows that spot evolution has a much bigger effect on the falsepositive rate than adding Poisson noise to the light curves.
According to the above results our method exhibits a falsepositive rate of less than ~12%. Relaxing the constraints imposed on the method, especially allowing for smaller peak separations (see Sect. 3), leads to higher falsepositive rates, as discussed in Sects. 4.3 and 5.
Fig. 4 Distributions of α_{model} (solid black) and α_{obs} (dashed red) of the 10spot model for different noise levels. Upper panel: no spot evolution, lower panel: spot evolution included. Falsepositive detections exhibiting α_{obs}< 0 are located left of the dashed blue line (indicating α_{obs} = 0). 
4.3. Kepler data
In this section we apply our method to a sample of 50 selected Kepler stars. As shown previously our method performs well for simulated data. We were curious to see whether it is possible to detect the sign of DR in real data in general, especially answering the question for how many stars antisolar DR is revealed by our method. Therefore, we selected a handpicked sample of solarlike stars from the Kepler archive. Although solarlike DR is expected for these stars, the underlying differential rotation law is still unknown.
The selected stars should hold the following criteria: surface gravities log g> 3.5, color index 0.60 <B − V< 0.92 (roughly corresponding to spectral type G0–G9 using the g − r to B − V color transform from Jester et al. 2005), and Kepler magnitude kepmag< 13, leaving roughly 5700 stars. Since we are only interested in rotation, the above sample was further reduced by discarding eclipsing binaries^{2} and KOIs^{3}. Two stars in the sample were characterized as hybrid pulsators (Uytterhoeven et al. 2011), and one star as γ Doradus pulsator (Tkachenko et al. 2013). These stars were also discarded.
To minimize the effect of possible instrumental flaws we only selected stars with strong intrinsic variability, expressed in terms of the variability range R_{var}> 0.3% (Basri et al. 2011). Additionally, the most significant period in the periodogram should hold 1 <P_{1}< 25 days. These limits should guarantee a stable rotational modulation over many rotation cycles.
We only used Kepler Q1–Q14 data, which were consistently reduced with the PDCMAP pipeline (version 2.1–3.0). Q15−Q17 data have not been considered because these data were reduced with the msMAP pipeline, which applies a 20day highpass filter to certain quarters, thus affecting the signal of slower rotators. The light curves were stitched together in a simple way by dividing each quarter by its median and subtracting unity. The light curves were smoothed with a 30day window to remove remaining longterm trends, likely caused by the instrument.
Furthermore, we visually inspected all light curves and their periodograms. We only selected stars with multiple periodogram peaks, at least one visible harmonic, and clearly visible rotational modulation in the light curve. The above criteria and the visual inspection left only 50 stars in total.
Our method described in Sect. 3 was applied using the same constraints as for the simulated data. We found that the separation of the two peaks P_{1} and P_{2} on the frequency grid has a strong influence on the determination of the correct sign of DR. Therefore, we define (5)To have a dimensionless quantity for the peak separation, we divide this frequency difference by the minimum frequency f_{min} = 1/max(time) of the periodogram. In the following, the number of points on the frequency grid separating two peaks is called the peak separation. Forcing the program only to pick peaks with Δf/f_{min} ≥ n with n ≥ 2 rules out certain peak combinations. This happens, e.g., when more than two distinct rotation peaks are found, but only the two closest peaks exhibit a detectable harmonic. For this reason the total number of detections drops when increasing the minimum peak separation (s. second column in Table 2).
Table 2 summarizes the results and Fig. 5 shows the distribution of α_{obs} using the different constraints. Increasing the minimum peak separation significantly reduces the total number of detections, especially the cases where the highest peaks have been chosen (k_{min} = 1). The number of detections of antisolartype differential rotation varies from 28–37%.
To put the Kepler results into context we applied our method to the 10spot models including spot evolution and 1000 ppm white noise using the same peak separations as above. We found that reducing the minimum peak separation from five to two frequency bins increases the percentage of falsepositives from 11.3% to 20%. The number of cases with k_{min} = 1 fluctuates between 78–91%. This might suggest that the result for Δf/f_{min} ≥ 2 is least reliable, although k_{min} = 1 was used in all cases. Such split peaks with small peak separation are likely caused by spot evolution.
Number of solar and antisolar detections out of the 50 Kepler stars.
Fig. 5 Distribution of α_{obs} of the 50 Kepler stars forcing different minimum peak separations. The dashed blue line indicates rigid rotation (α_{meas} = 0). 
5. Discussion
As shown in the previous section, our method performs well for simulated data, especially the results of the 2spot models are striking. The detection rate for the 10spot models is still high, although the periodograms exhibit more than two significant rotation peaks, as do their harmonics (see Fig. 2).
The falsepositive rate of our method depends on the noise level, the minimum peak separation, and the presence of spot evolution. We found that our method is very robust against adding pure white noise to the light curves, in agreement with other studies (Reinhold & Reiners 2013). Adding higher white noise to the light curves decreases the total number of detections, but does not significantly increase the falsepositive rate. Spot evolution was found to be the major contributor of spurious detections. The largest rate was found for the 10spot models including spot evolution and having small peak separations, increasing up to 20% for the highest noise levels. In certain cases very shortlived spots can produce split rotation peaks, and thus mimic differential rotation.
In real data this effect becomes more important because spot lifetimes are difficult to assess. In some situations spot evolution can barely be distinguished from DR, making the determination of the sign of DR is even harder. In general, periodograms of Kepler stars exhibit many more peaks, with some of them being spurious detections. Furthermore, the periodogram strongly depends on the length of time series. Since the peak width is proportional to the inverse of the total observing time span, a certain length of highquality photometry is required to resolve individual peaks. However, the analysis of two different light curve segments of equal length may also change the result. Spot evolution and differential rotation may be more “visible” in one segment than in the other, as revealed by their periodograms, either showing a single or multiple clear peaks, respectively. This effect immediately connects to the different peak separations (see Fig. 5). Computing the periodogram of the full time series averages out these effects. Therefore, analyzing different light curve segments might increase the sample of suitable candidates to test this method.
Unfortunately, the falsepositive rate of our method cannot be transferred directly to Kepler data because the percentage of real antisolar rotators is unknown. It should rather be interpreted as the number of cases where true solar rotators are misclassified as antisolar rotators, and vice versa. Since the light curves and periodograms of solar and antisolar rotators exhibit similar shapes, one can argue that both cases should be treated equally by the method, i.e., the rate of stars with wrongly determined sign of DR should be the same on either side. The measured number of antisolar detections then consists of the fraction of real antisolar rotators, which are actually detected as such, plus the number of migrated solarlike stars, misclassified as antisolar rotators. Thus, we define S′ and A′ are the measured solar and antisolar detections, respectively, S and A the number of real solar and antisolar rotators, respectively, and f = 0.2 the maximum falsepositive rate. Solving this system of linear equations yields Applying this solution to the numbers given in Table 2 we detect solar DR in 21–34 and antisolar DR in 5–10 stars of our sample of 50 stars. Table 3 lists the KIC numbers of our sample of 50 Kepler stars, providing their differential rotation law found by our method for the different peak separations. Most of the stars are consistently found as solar (s) or antisolar (a) rotators, whereas some switch from one to the other. Additionally, one can see that forcing a larger peak separation leads to a nondetection (–) for some stars.
As discussed above this new method performs well in our simulations, but should be used with caution in real data. Although we tried to successively approach the shape of real light curves by increasing the number of spots, probing the effect of spot evolution, and considering different noise scenarios, some effects have not been considered. The spot evolution as implemented in our model may be overestimating real spot lifetimes. Furthermore, instrumental artifacts in the quarterly delivered Kepler data, and their reduction by the PDCMAP pipeline certainly affect the overall light curve shape, and thus the period determination.
This number of antisolar detections seems quite large, but reflects the outcome of the method. Since the light curves created from the 10spot models look very similar to real Kepler data, we are optimistic that our findings of antisolar DR are real, although independent measurements to either prove or disprove our results are urgently needed.
Online material
Solar (s) and antisolar (a) detections for the selected 50 Kepler stars depending on the peak separation.
We use the normalization from Eq. (22) in Zechmeister & Kürster (2009).
Acknowledgments
T.R. acknowledges the support from Prof. Laurent Gizon.
References
 Ammlervon Eiff, M., & Reiners, A. 2012, A&A, 542, A116 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Basri, G., Walkowicz, L. M., Batalha, N., et al. 2011, AJ, 141, 20 [NASA ADS] [CrossRef] [Google Scholar]
 Bonanno, A., Fröhlich, H.E., Karoff, C., et al. 2014, A&A, 569, A113 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dikpati, M., & Cally, P. S. 2011, ApJ, 739, 4 [NASA ADS] [CrossRef] [Google Scholar]
 Donati, J.F., & Collier Cameron, A. 1997, MNRAS, 291, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Gastine, T., Yadav, R. K., Morin, J., Reiners, A., & Wicht, J. 2014, MNRAS, 438, L76 [NASA ADS] [CrossRef] [Google Scholar]
 Gizon, L., & Solanki, S. K. 2004, Sol. Phys., 220, 169 [NASA ADS] [CrossRef] [Google Scholar]
 Guerrero, G., Smolarkiewicz, P. K., Kosovichev, A. G., & Mansour, N. N. 2013, ApJ, 779, 176 [NASA ADS] [CrossRef] [Google Scholar]
 Jester, S., Schneider, D. P., Richards, G. T., et al. 2005, AJ, 130, 873 [NASA ADS] [CrossRef] [Google Scholar]
 Käpylä, P. J., Käpylä, M. J., & Brandenburg, A. 2014, A&A, 570, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Karak, B. B., Käpylä, P. J., Käpylä, M. J., & Brandenburg, A. 2015, A&A, in press, DOI: 10.1051/00046361/201424521 [Google Scholar]
 Kitchatinov, L. L., & Rüdiger, G. 2004, Astron. Nachr., 325, 496 [NASA ADS] [CrossRef] [Google Scholar]
 Kővári, Z., Bartus, J., Strassmeier, K. G., et al. 2007, A&A, 474, 165 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kővári, Z., Kriskovics, L., Künstler, A., et al. 2015, A&A, 573, A98 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lanza, A. F., Das Chagas, M. L., & De Medeiros, J. R. 2014, A&A, 564, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lund, M. N., Miesch, M. S., & ChristensenDalsgaard, J. 2014, ApJ, 790, 121 [NASA ADS] [CrossRef] [Google Scholar]
 Reiners, A., & Schmitt, J. H. M. M. 2002, A&A, 384, 155 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Reinhold, T., & Reiners, A. 2013, A&A, 557, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Reinhold, T., Reiners, A., & Basri, G. 2013, A&A, 560, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Strassmeier, K. G., Kratzwald, L., & Weber, M. 2003, A&A, 408, 1103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Tkachenko, A., Aerts, C., Yakushechkin, A., et al. 2013, A&A, 556, A52 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Uytterhoeven, K., Moya, A., Grigahcène, A., et al. 2011, A&A, 534, A125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Vida, K., Kővári, Z., Švanda, M., et al. 2007, Astron. Nachr., 328, 1078 [NASA ADS] [CrossRef] [Google Scholar]
 Vida, K., Oláh, K., & Szabó, R. 2014, MNRAS, 441, 2744 [NASA ADS] [CrossRef] [Google Scholar]
 Weber, M., Strassmeier, K. G., & Washuettl, A. 2005, Astron. Nachr., 326, 287 [NASA ADS] [CrossRef] [Google Scholar]
 Zechmeister, M., & Kürster, M. 2009, A&A, 496, 577 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
All Tables
Solar (s) and antisolar (a) detections for the selected 50 Kepler stars depending on the peak separation.
All Figures
Fig. 1 Upper panel: simulated light curve of a star with two spots and an equatortopole shear α = 0.085 (solarlike DR) according to Eq. (1). Lower panel: same parameters as above but with a negative value of α = −0.085 instead (antisolar DR). 

In the text 
Fig. 2 Upper panel: light curve of a randomly selected star from the 10spot models without spot evolution. Lower panel: periodogram of the above light curve. The peaks of the pairs () and () are clearly visible and emphasized by black upside down and normal triangles, respectively. 

In the text 
Fig. 3 Distributions of α_{model} (solid black) and α_{obs} (dashed red) of the 2spot model for different noise levels. Upper panel: no spot evolution, lower panel: spot evolution included. The dashed blue line indicates α_{obs} = 0. 

In the text 
Fig. 4 Distributions of α_{model} (solid black) and α_{obs} (dashed red) of the 10spot model for different noise levels. Upper panel: no spot evolution, lower panel: spot evolution included. Falsepositive detections exhibiting α_{obs}< 0 are located left of the dashed blue line (indicating α_{obs} = 0). 

In the text 
Fig. 5 Distribution of α_{obs} of the 50 Kepler stars forcing different minimum peak separations. The dashed blue line indicates rigid rotation (α_{meas} = 0). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext 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 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.