A study of pulsation & rotation in a sample of A-K type stars in the Kepler field

We present the results of time-series photometric analysis of 15106 A-K type stars observed by the Kepler space mission. We identified 513 new rotational variables and measured their starspot rotation periods as a function of spectral type and discuss the distribution of their amplitudes. We examined the well-established period-color relationship that applies to stars of spectral types F5-K for all of these rotational variables and, interestingly, found that a similar period-color relationship appears to extend to stars of spectral types A7 to early-F too. This result is not consistent with the very foundation of the period-color relationship. We have characterized 350 new non-radial pulsating variables such as A- and F-type candidate \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\delta $\end{document}δ Scuti, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\gamma $\end{document}γ Doradus and hybrid stars, which increases the known candidate non-radial pulsators in the Kepler field significantly, by ∼20%. The relationship between two recently constructed observables, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\mathit{Energy}$\end{document}Energy and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\mathit{Efficiency}$\end{document}Efficiency, was also studied for the large sample of non-radial pulsators, which shows that the distribution in the logarithm of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\mathit{Energy}$\end{document}Energy (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\log (\mathit{En})$\end{document}log(En)) can be used as a potential tool to distinguish between the non-radial pulsators, to some extent. Through visual inspection of the light curves and their corresponding frequency spectra, we found 23 new candidate red giant solar-like oscillators not previously reported in the literature. The basic physical parameters such as masses, radii and luminosities of these solar-like oscillators were also derived using asteroseismic relations. Electronic Supplementary Material The online version of this article (10.1007/s10509-018-3480-1) contains supplementary material, which is available to authorized users.


Introduction
In the past, many photometric and spectroscopic surveys were conducted in both the Northern and Southern hemisphere to search for low amplitude variability in A-to F-type stars (Martinez et al. 1991;Nelson and Kreidl 1993;Paunzen et al. 2012;Kochukhov 2006;Kurtz et al. 2006). At the Aryabhatta Research Institute of Observational Sciences (ARIES), Nainital, India, we initiated a similar research project to detect low-amplitude photometric variability and study it in more detail using spectroscopy (Ashoka et al. 2000;Martinez et al. 2001;Joshi et al. 2003Joshi et al. , 2006Joshi et al. , 2009Joshi et al. , 2010Joshi et al. , 2012Joshi et al. , 2016Joshi et al. , 2017.
The best precision ever obtained from ground-based photometric observations is 14 µmag (Kurtz et al. 2005) which is not even sufficient to detect the amplitude variation generally found in solar-like stars. After the launch of various space missions, much lower detection levels of the order of µmag have been achieved: WIRE (Buzasi et al. 2000), MOST (Matthews 2004), CoRoT (Baglin et al. 2009

), and
Kepler Koch et al. 2010). This motivated us to extend our survey program from ground-based to space-based observations. This resulted in the formation of a new collaboration with South African astronomers to detect and study photometric variability using the Kepler archival data. The first two papers of this collaboration are published by Balona et al. (2013aBalona et al. ( , 2016. In this third paper, we study the time-series photometric data of 15106 A-K type stars observed by the Kepler space mission.
Space-based observations of variable stars have revolutionized the field of variability studies. Dedicated missions such as CoRoT and Kepler have duty cycles which are unachievable from the ground. It effectively solves many problems prevalent in ground-based observation campaigns. In particular, the Kepler mission in its original operational mode collected data of unprecedented photometric precision for almost 200000 stars in a large field of 105 square degrees in the direction of Cygnus and Lyra. In addition to its primary objective of detecting planetary transits, Kepler has provided high quality data for many areas of investigation, particularly asteroseismology and rotation-related activity. We have taken advantage of the photometric quality of data from the Kepler mission to study variable A-K type stars in the original Kepler field.
The long cadence mode (LC) of Kepler (one observation every 29.4 min) is perfectly adequate for detecting variability in red giants. Solar-like oscillations have been discovered in thousands of them Huber et al. 2010;Kallinger et al. 2010a), including some in open clusters (Hekker et al. 2011;Stello et al. 2010;Basu et al. 2011). It is possible to independently estimate the mass M, radius R and surface gravity log g of each star, in the case of pulsating red giants. In these stars, the mean density, mass and radius can be obtained from the large separation ν (i.e. frequency spacing between modes of consecutive overtones n but of same spherical harmonic degree ) and frequency of maximum amplitude ν max . Stello et al. (2010) was the first to apply this method to red giants in the cluster NGC 6819, using Kepler data.
In addition to its primary objective of detecting planetary transits, Kepler's high quality data provide a valuable opportunity for specifically investigating relationships between several areas of research in stellar structure and evolution. The main sequence δ Scuti (δ Sct) and γ Doradus (γ Dor) stars with masses between 1.2 and 2.5M are particularly useful for pulsation studies, to infer the interior structure of stars and to test theoretical models describing them. The δ Sct stars pulsate in pressure modes (p-modes), giving us the opportunity to study the stellar envelope while the γ Dor stars pulsate in gravity modes (g-modes), that help to study the regions near the stellar core. Hybrid objects are pulsating in both types of modes and are hence of great asteroseismic interest because they can give constraints on the full stellar structure. Debosscher et al. (2011) used automated algorithms to classify the variable stars in the Kepler field. On combining the first two quarters of data for hundreds of variable stars, Grigahcène et al. (2010) showed that the frequency spectra are so rich that there are practically no pure δ Sct or γ Dor pulsators. Essentially all the stars in the corresponding area of the Hertzsprung-Russell (H-R) diagram show pulsational frequencies in both the δ Sct and γ Dor frequency ranges. To categorize these non-radial pulsators, the authors proposed a new observational classification scheme, taking into account the amplitude of the dominant mode as well as the range of observed frequencies.
This paper reports on the outcome of a study for which the key goals are threefold: (i) to identify additional solar-like oscillators. For the solar-like oscillators found in our sample, we derive their basic parameters (e.g. mass, radius and luminosity) using the asteroseismic scaling relations. (ii) to explore the nature of a rotation period-color relation (if any) that might extend blueward of the canonical limit at (B − V ) = 0.47 (see Mamajek and Hillenbrand 2008). This is motivated by Balona et al. (2013a,b), who found a continuation of the canonical period-color relation to stars with (B − V ) > 0.2, which demands a re-evaluation of the physical cause of this relationship (or a re-evaluation of the presence of convection in Atype stars). We therefore searched for spotted stars in our sample to investigate this relationship. (iii) to search for other non-radial pulsators in our sample, using some automated constraints and independently confirming them by visual inspection of all the light curves and corresponding Discrete Fourier transforms (DFTs). Finding new non-radial pulsators is imperative to advance our understanding of such variable stars, through asteroseismology. Two recently constructed observables, Energy and Efficiency  have also been studied for our large sample of non-radial pulsators to explore the relationship between δ Sct, γ Dor and candidate hybrid variables in the parameter space of these new observables.
Apart from the DFT technique, the analysis of variance (AoV) method (Schwarzenberg-Czerny 1989, 1996 was also adopted to estimate the period of variability and verify the classification in all cases. We study the Kepler field stars using the Sloan g, r, i, z photometry retrieved from the Kepler Input Catalogue (KIC; Brown et al. 2011) and timeseries photometry from the Kepler archive at the Mikulski Archive for Space Telescopes (MAST). The work presented here is organized as follows: In Sect. 2 we give the sample selection and description of the Kepler data archive. In Sects. 3, 4 and 5 the different classes of pulsating and rotating stars are discussed, followed by the conclusion drawn from our study in Sect. 6.

Data description
The data presented in this work include all available LC data from quarter 0 (Q0) to quarter 17 (Q17) for 15106 Kepler light curves. We selected stars of spectral types ranging from A0 to K7 as this range included the cooler red giant stars where solar-like oscillations are likely to be present. The hotter stars in the spectral region (B − V ) < 0.47 were also included for the investigation of the period-color relationship in rotational variables. The complete list of stars that were studied in this paper is provided online.
For the vast majority of stars, Kepler photometry is available only in the LC mode, which consists of practically uninterrupted exposures with a cadence of 29.4 minutes. For a few thousand stars, the short cadence (SC) mode of ∼1-minute exposures is available, but these usually cover only one or two months of observations. Kepler observed in white light, with a passband with a FWHM covering the 430-890 nm wavelength range. The observed stars have Kepler magnitudes ranging from 9 to 16 mag. The Kepler data are available in either 'uncalibrated' or 'calibrated' form. The calibrated data suffer from artefacts caused by the processing and are not used here. Details of the technique used to correct the uncalibrated data and how the corrections affect the derived frequencies can be found in . As stated in , the trend-removing procedure for Kepler data ". . . tends to dampen or remove very low frequencies. . . ", such that, in practice, the lowest frequency that can be detected in these data is about 0.1 d −1 .  go on to say "but frequencies above about 0.1 d −1 are not affected". We took this as a vote of confidence in the determination of rotational periods shorter than 10 days and therefore confined our study to rotational periods below this upper limit. The Nyquist frequency of the Kepler LC data is around 24 d −1 . All the data used in this paper are publicly available and can be retrieved using MAST.

Solar-like oscillations
Solar-like oscillations can be expected in low-mass mainsequence stars, sub-giants, stars on the red-giant branch (RGB), horizontal branch and asymptotic-giant branch stars (Christensen-Dalsgaard and Frandsen 1983;Houdek et al. 1999;Dziembowski et al. 2001). Following the success of helioseismology, the large volumes of long time-series data of exquisite quality, in particular from the CoRoT and Kepler space missions, have been particularly important for the study of solar-like oscillations in stars, owing to the small amplitudes of these oscillations. In particular, after the launch of Kepler in 2009, there has been tremendous progress in asteroseismology of G-and K-type red giants, with the detection of solar-like oscillations in thousands of stars. These oscillations are easily identified in the periodogram because of the localized comb-like structure where amplitudes decrease sharply from a central maximum.
Following the approach by Kallinger et al. (2010b) we use the values of ν, ν max and the effective temperatures from the Kepler Input Catalog in combination with the scaling relations as described by Brown et al. (1991), Kjeldsen and Bedding (1995), and Chaplin et al. (2008) to compute the masses and radii of the stars directly: where ν max = 3120 µHz and ν = 134.88 µHz are solar values as determined by Kallinger et al. (2010a). From the radii and T eff we computed the luminosity as L ∝ R 2 T 4 eff . ν max values were derived by fitting a Gaussian to the highest peaks in the periodogram at an appropriate frequency range of the comb-like structure, after removing a linear fit to the background noise level. For deriving ν, we used the autocorrelation function of the periodogram over the range of the oscillations (Huber et al. 2009). The frequency spacing between consecutive l = 0 and 1 modes also gives rise to peaks in the autocorrelation function at 1 2 ν. We selected those peaks closest to the value given by with a = 0.266 ± 0.004 and b = 0.761 ± 0.004.
This empirical relationship is based on Kepler data of 662 field red giants (Hekker et al. 2011). As an example of our analysis, the periodogram of KIC 002568575, together with the fitted Gaussian envelope and the autocorrelation function, is shown in Fig. 1. The values of ν max and ν derived from our fits were used to calculate mass, radii and, consequently, luminosities, for the 23 newly-identified red giants for which we detected solar-like oscillations. Their  Table 1. We have also shown the radius and luminosity from the GAIA Data Release 2 (DR2) (Andrae et al. 2018) in the same table. We compared our ν max and ν values with Yu et al. (2018) for the common objects, and calculated the percentage difference for every common object by dividing the difference with our estimated values. We then took a standard deviation of these percentages differences and used it to find the propagation errors in radius and luminosity. The radius and luminosity, along with their roughly estimated errors are then plotted against GAIA DR2 values, as shown in Fig. 2. There is no significant difference between our estimated radius and luminosity from the solar-like oscillations and GAIA DR2 radius and luminosity. The average difference between our calculated and GAIA DR2 values, is 5% in radius and 3% in luminosity. Figure 3 shows the relationship between ν max and ν for our stars showing solar-like oscillations and how well they fit the known empirical relation given in Eq. (3). The increased spread around the correlation at the high-end of the ν max relation is consistent with the expected spread, which increases with ν max (see e.g. Hekker et al. 2011;Stello et al. 2009a;Huber et al. 2010). Since the large separation, frequency of maximum amplitude, luminosity and the effective temperature are related to each other, one can use stellar evo-lution and pulsation models to constrain the solution. The stellar parameters are usually estimated by searching among a grid of stellar models to get a best fit for observed values of ν, ν max , T eff and metallicity. This is usually referred to as 'grid' asteroseismology (Stello et al. 2009b;Basu et al. 2010;Quirion et al. 2010;Gai et al. 2011).
In our study, we have used the direct method to estimate the stellar parameters. This so-called direct method provides estimates that are independent of stellar evolutionary theory. The luminosities from the KIC and asteroseismic luminosity log(L/L ) ν derived using scaling relations are not in good agreement because the estimated luminosity from the scaling relations is highly dependent on the accuracy of T eff . A typical error of 150 K in the determination of T eff leads to an error of about 0.2 mag in the bolometric absolute magnitude. Therefore, spectroscopic observations are essential for the accurate determination of T eff and log g (Catanzaro et al. 2011;Niemczura et al. 2015). The error in mass, radius and luminosity also depends on the error in estimating ν max and ν. We manually selected the fits for each star that describes the overall shape of the envelope very well. It is difficult to estimate the uncertainties in ν max and ν. Hekker et al. (2012) have investigated the uncertainties in these two quantities by means of simulations using different methods and found that for a time series of duration about Table 1 A list of 23 new stars with solar-like oscillations discovered by visual inspection. Kp is the Kepler Magnitude. The first five columns are KIC parameters, except the second column, which is revised effective temperature from Pinsonneault et al. (2012). The calculated parameters i.e., frequency of maximum amplitude, the large separation, mass, radius and luminosity calculated from the solar-like oscillations are shown from sixth to tenth column, respectively. The last two columns are stellar parameters from GAIA DR2. The complete table is available online  Fig. 3 The relationship between ν max and ν for stars given in Table 1 overplotted with the function ν ≈ 0.266ν 0.761 max (solid black line), with the 1 sigma errors included in blue and 3 sigma errors included in magenta 500 d, the uncertainty in ν max is about 2% and in ν it is about 1%.
Comparative studies between the direct method and independently determined properties from binaries, shows that the estimated parameters are found to be consistent at the level of precision of the uncertainties, i.e. up to 10% or better Miglio 2012). Silva Aguirre et al.
(2012) used the brightest solar-like oscillators from the Kepler data, and found excellent agreement between scaling relation inferred radius and those inferred from using Hipparcos parallaxes (at the level of a few percent). Huber et al. (2012) used combined Hipparcos parallaxes and interferometric observations of some of the brightest Kepler and CoRoT targets and also found excellent agreement with the stellar radii inferred from scaling relations, at the 5% level. Gaulme et al. (2016) compared the masses and radii of 10 red giants obtained by combining the radial velocities and eclipse photometry with the estimates from the asteroseismic scaling relations, and found that the asteroseismic scaling relations overestimate red-giant radii by about 5% on average and masses by about 15% for stars at various stages of red-giant evolution.

δ Scuti, γ Doradus and hybrid stars
A detailed analysis of stellar pulsation frequencies can be used to infer interior stellar structure and test theoretical models. The main-sequence γ Dor and δ Sct stars are particularly useful for asteroseismic studies. The δ Sct stars pulsate in both radial and non-radial p-modes with periods ranging from 18 min to 7 hrs (see Breger 2000 and Bowman and Kurtz 2018, for a review). The pulsations in the δ Sct variables are driven by the κ-mechanism operating in the HeII partial ionization zone with some additional contribution from the HI ionization zone. The γ Dor variables are generally cooler than δ Sct stars and are mostly located near the cool edge of the δ Sct instability strip. The pulsations in γ Dor variables are characterized by high-order, low-degree and multiple non-radial g-modes with periods of 0.3 to 3 days (Kaye et al. 1999;. The pulsations in these stars are excited by the convective blocking mechanism at the base of their envelope convection zone (Guzik et al. 2000). The distinction between the two classes is clearer if we consider the value of the pulsation constant, Q, where all γ Dor stars have Q > 0.09 d and all δ Sct stars have Q < 0.09 (see Fig. 9 of Handler and Shobbrook 2002, for details), accommodating possible uncertainties in the determination of Q (Breger 1990). Note that it is not always obvious to distinguish γ Dor variables from spotted stars as the expected pulsation and rotation periods are of the same order.
Due to the overlapping of T eff and log g ranges for these two types of variables, there is a possibility for stars showing both types of pulsation simultaneously, i.e. δ Sct-type combined with γ Dor-type pulsations, and are thus known as the A/F-type hybrid stars. Before the launch of the Kepler and CoRoT space missions, only a few stars pulsating in both p and g-modes were discovered and interpreted using theoretical models (Bouabid et al. 2009). Ground-based observations discovered only 4 hybrid stars (Henry and Fekel 2005;Uytterhoeven et al. 2008;Handler 2009). The first published analysis on 234 targets by the Kepler Asteroseismic Science Consortium (KASC) confirmed the pulsations of either type (δ Sct or γ Dor) and revealed the hybrid behavior in essentially all of them (Grigahcène et al. 2010). In a study of 750 A-F stars observed for four quarters, 475 stars showed either δ Sct or γ Dor variability, and 36% of these were hybrids ). However, Balona (2014) found that essentially all δ Sct stars in the Kepler field contain low-frequency modes typically associated with g-modes and concluded that all δ Sct stars should therefore be considered as 'hybrids'. He also concluded that the presence of δ Sct pulsations is strongly connected with the low frequencies in δ Sct stars. In accordance with this finding, it is clear that so-called δ Sct/γ Dor hybrids (according to Grigahcène et al. 2010) actually occur throughout the δ Sct instability strip. This supports Balona's assertion that the δ Sct/γ Dor hybrid classification is redundant. In this paper, we classify δ Sct/γ Dor stars as δ Sct variables and we still refer to stars with their dominant frequencies in the low (g-mode) range but also displaying high-frequency modes as γ Dor/δ Sct hybrids. Hybrid stars are among the most interesting and vital targets for asteroseismology for providing additional constraints on stellar structure, because the γ Dor stars pulsate in g-modes which have high amplitudes deep in the star and allow us to probe the stellar core, while the p-modes, efficient in δ Sct stars, have high amplitudes in the outer regions of the star and probe the stellar envelope. Using the Kepler archive time-series data, Kurtz et al. (2014) determined the rotation rate in the deep core and surface of the main sequence hybrid star KIC 11145123 and found that the surface rotates slightly faster than the core. However, our understanding of the structure of the outer layers, the convective core and the outer convective zone is still far from complete, hence a better understanding of the driving mechanisms in these stars is very important.
Some major previous efforts to classify the non-radial pulsators in the Kepler and K2 fields include those mentioned in Table 2. Combining all these studies results in around 1893 candidate non-radial pulsators in the nominal Kepler field (counting the overlapping objects once). This is an approximation, since it includes the 403 candidate δ Sct/β Cephei (β Cep) stars and 441 candidate γ Dor/slowly pulsating B (SPB) stars given by Debosscher et al. (2011), who could not reliably distinguish between δ Sct and β Cep stars, nor between γ Dor and SPB stars, because their classifier only used information from the Kepler light curves. Moreover, since there are no other previous major studies dedicated to the classification of non-radial pulsators in the Kepler field, we are certain that the error in our approximation is small. Very recently, Ibanoglu et al. (2018) compiled a comprehensive catalog of γ Dor pulsators, listing 109 genuine variables, 291 candidate stars in the CoRoT fields, 318 candidate variables in Kepler field including 11 objects that are in a binary system. They have included one star which was classified as binary+hybrid by Uytterhoeven et al. (2011), in their list of 11 candidate binary+γ Dor variables. In addition, they also compiled a list of 233 hybrids of which 205 are in the Kepler field.
We identified the candidate δ Sct and γ Dor variables from a large sample of Kepler field stars, using automated constraints (frequency and amplitude), followed by a visual inspection of all the light curves and the frequency spectra (DFT). In addition to frequency, the pulsation amplitudes were also used to potentially filter out the high degree modes that generally would be expected to have a lower photometric amplitude. We also used the AoV algorithm in multi-harmonic mode (Schwarzenberg-Czerny 1989, 1996 to estimate and confirm the periods of these stars, but list the periods derived from the DFT in Table 3. We characterized the stars as δ Sct if most of the detected frequencies in the DFT are between 5 d −1 and 24 d −1 (being the Nyquist frequency for LC Kepler data). Stars where the lower frequencies (less than 5 d −1 ) are also present, are also classified as δ Sct variables. Stars for which the detected frequencies were between 0.30 and 5 d −1 are characterized as γ Dor variables. By carefully examining the periodograms of δ Sct stars and γ Dor stars we were also able to identify the γ Dor/δ Sct hybrids.
This search led to the discovery of 168, 110, and 72 new candidate δ Sct stars, γ Dor stars, and γ Dor/δ Sct variables, respectively. They are listed in Table 3 and their positions in the H-R diagram are shown in Fig. 4. Since we included the δ Sct/γ Dor variables in the δ Sct variable definition, we do not compare our number of candidate δ Sct variables to those mentioned in the literature. Comparing the number of candidate non-radial pulsators detected in our study with the    Bouabid et al. (2013) showed that some stars classified as a potential hybrid star from the frequency distribution could in fact be rapidly rotating γ Dor stars whose g-modes have been shifted to higher frequencies. To confirm the classification of the hybrid stars and to determine the nature of their pulsations as p-or g-modes, multi-color photometry and/or a spectroscopic study is required for the determination of the atmospheric parameters (Catanzaro et al. 2010). The candidate δ Sct, γ Dor and γ Dor/δ Sct hybrids found in our study could be the potential targets for detailed ground-based photometry and mid-to high-resolution spectroscopic observations because the fundamental parameters (T eff , log g, v sin i, and log(L/L )) are essential for computing asteroseismic models and interpreting Kepler data.
To understand the relationship between δ Sct, γ Dor and the hybrid pulsators, Uytterhoeven et al. (2011) first introduced two new empirical observables, Energy (related to the energy contained in the oscillation) and Efficiency (related to the convective efficiency of the outer convective zone), using the available observational data: where A max and ζ max refer to the highest amplitude mode of the star (in ppm) and its frequency (in d −1 ), respectively. Uytterhoeven et al. (2011) found that the distribution in the logarithm of Efficiency and Energy peaks at different values for δ Sct and γ Dor stars. We further examined the peak values of the distribution in the logarithm of Efficiency (log(Eff )), Energy (log(En)) and Amplitude (log(Amp)), for our large sample of non-radial pulsators, where Amplitude here refers to the highest amplitude mode of the star in the DFT of the light curve. The results are shown in Table 4. The values in the table are derived by fitting a gaussian to the histograms. The distributions of log(En), log(Amp), and log(Eff ) along with the fitted gaussians are shown in Fig. 5. The majority of δ Sct stars have log(En) > 7.5 while, for the γ Dor variables, the majority of stars in the distributions almost all have log(En) < 7.5. The log(En) versus log(Eff ) diagram is plotted in Fig. 6. Although not much can be inferred about these observables and its physical basis at this moment, it is certain that these observables peak at different values in their respective distributions (see Fig. 5 and Table 4). As a result, it can be used to some extent, as a method to distinguish the non-radial pulsators, especially using the distribution in the logarithm of Energy (log(En)). Our study confirmed and improved the previous results obtained by Uytterhoeven et al. (2011), who used fewer non-radial pulsators in deriving the distribution peaks, based on their new observables.

Rotational variables
It was found a long time ago that the hot B stars are the most rapidly rotating main sequence stars and that rotation slows down for cooler stars. For stars cooler than F5, the rate of rotation drops below 10 km s −1 which is difficult to measure. The reason why rotation is so slow for cool stars like the Sun is thought to be due to convection. The stars cooler than F5 have outer convective envelopes, while hotter stars have radiative envelopes. Convection is thought to be necessary for the generation of a magnetic field by the dynamo mechanism (Reiners 2012). If a star has a magnetic field with surface convective envelopes, that means the ionized gases (which are conducting) are trapped by the magnetic lines of flux. Ionized gas in the stellar wind is forced to move along the magnetic field lines, which carry away the angular momentum and hence there is a braking action on rotation. Another product of a magnetic field is that it can create starspots and flares. Therefore one only expects stars cooler than F5 to have starspots and to show flares.
It was shown by Skumanich (1972) that the rotation periods of solar-type stars decrease over time, t, such that the rotational velocity, v rot ∝ 1/ √ t . In a similar way, the chromospheric activity, which is a proxy for magnetic field strength,

Fig. 5 Distribution of Efficiency (log(Eff )), Energy (log(En)) and
Amplitude (log(Amp)). The first and third row shows the distribution for δ Sct and γ Dor variables, respectively. Distribution for γ Dor/δ Sct hybrids are shown in the second row. The Y-axis represents the number of stars belonging to each bin(N). The dashed black line is the gaussian fit to each of the histogram, for the peak estimation also decreases over time. The cause of the Skumanich law is believed to be angular momentum loss due to stellar winds, but the exact dependence of v rot on age is still not understood properly. However, this has not hindered the development and use of 'gyrochronology', a technique to determine the ages of field stars which is based upon the period-age-mass (PtM) relationship (Barnes 2007;Mamajek and Hillenbrand 2008;Kawaler 1989;Barnes and Kim 2010).
Using observed rotation periods, it was shown by Barnes (2001), that the age dependence of rotation for these stars is of Skumanich type (P ∝ √ t) and that the correlation of stellar mass (or, more specifically, its proxy (B −V ) 0 ), with surface rotation rate that is seen in the Hyades cluster, also extends to stars of types F5 and later in general. Barnes (2003) demonstrated the consistency of this correlation in a detailed study including many open clusters, finding that the correlation becomes tighter with increasing cluster age. The PtM relationship is not unique for stars of young (≤100 Myr) clusters, and two (fast and slow) branches emerge from the PtM surface (Barnes 2003). The fast branch vanishes with increasing age and the stars come together in their rotational evolution onto the slow branch, leading to an unique relationship between P and M at an age of about 600 Myr.
The Kepler data are ideal for measuring the rotational period of spotted stars. Rotation can be measured at the surface of individual stars using either spectroscopy or periodic variations in photometric light curves due to the presence of these starspots (Mosser et al. 2009). Previously, stellar rotation was measured using spots, rotational broadening of absorption lines (e.g., Kaler 1989) and chromospheric activity, specifically Ca II emission (Noyes et al. 1984). With the launch of the CoRoT and Kepler telescopes, the study of stellar rotation has seen significant advances and has reached a new level of understanding. The high precision data from these space missions allow to estimate the rotational period for a well-defined large sample of the slowly rotating stars with low-amplitude modulation. Some previous studies to derive the rotational periods focusing on a broader Kepler sample are those of: Reinhold et al. (2013) (RRB13 in what follows), with an emphasis on differential rotation, who derive the rotation period of ∼24000 stars using data from Q3; Nielsen et al. (2013), who measured the rotation period of ∼12000 stars with Q2-Q9 data, comparing their results with previous spectroscopic studies; McQuillan et al. (2014), who derived rotation periods for ∼34000 main-sequence stars of temperature less than 6500 K. Meibom et al. (2011), who exploited the Kepler data to determine rotation periods in the cluster NGC 6811, confirmed the existence of a unique functional relationship between rotation period, color and age. However the exact form of this relationship is yet to be investigated. More recently, Reinhold and Gizon (2015) derived rotation and differential rotation periods for more than 18500 and 12300 stars respectively, using different approaches, thereby confirming that the relative shear α (a measure of differential rotation), increases with rotation period for stars with T eff < 6700 K, while hotter stars show the opposite behavior. From the mean rotation periods and their uncertainties, they further infer stellar ages between 100 Myr and 10 Gyr for more than 17000 stars using different gyrochronology relations. Gyrochronology is calibrated using derived rotation period and age for cool stars of different masses, for stars in clusters with known ages. Until recently, it was not possible to measure rotation periods and test the gyrochronology relation for stars older than about one billion years, so model predictions were used to infer gyrochronology ages. Now, the rotation periods are known for stars in an open cluster of intermediate age (NGC 6819, 2.5 Gyr old, Meibom et al. 2015), and for old field stars with asteroseismologically determined ages. This confirmed the expected relationship between rotation period and stellar mass at the cluster age with a precision of order 10%, but failed to describe the asteroseismic sample (Angus et al. 2015). van Saders et al. (2016) essentially confirmed the presence of unexpectedly rapid rotation in stars that are more evolved than the Sun and showed that after incorporating dramatically weakened magnetic braking in the stellar evolutionary models, for old stars, both the cluster and asteroseismic data can be explained. This weakened braking has a significant impact on the gyrochronology relation for stars, including our Sun, that are more than midway through their main-sequence lifetimes.
We examined 15106 Kepler field stars showing light variations that could be attributed to starspots. This was done by examining the periodogram (DFT) and light curve of each star and assigning a variability type. The DFT is crossmatched with the AoV periodograms to make sure peaks are not spurious. Often the peaks are not single, but multiple. If the frequency ratios are equal to integer numbers, the extra peaks represent harmonics. Spots on a rotating stellar surface will exhibit harmonic frequencies because they act like low-level eclipses which consist of a flat light curve with dips. Since this is not a sinusoidal signal, harmonics will be present, depending on the spot location (its latitude), and the inclination angle. Sometimes, peaks are multiple without harmonics, which probably means there are multiple spots migrating with respect to each other and/or changing in size. We took the presence of harmonics as a good indicator of starspots. In doubtful cases, we also checked the folded light curves of the stars, based on the derived rotation frequencies.
In this way we can distinguish between pulsational variables of the γ Dor type and spotted stars.
We measured the rotation periods of 513 stars where we assume that the variation is due to starspots and rotation. The periodograms of some of the rotational variables are shown in Fig. 7. There is a possibility of contamination of light curves by background stars. However, out of a sample of 513 stars, it is not statistically possible that more than a minute fraction of these stars could have their values of (B − V ) 0 artificially raised by hotter (albeit more distant) background stars, for the following reasons: A background star will only rarely coincide closely enough with a target star to contaminate its light curve; background stars will be more distant and will have a small effect (if any) on a target star's light curve and the number density of stars falls off rapidly with temperature-so that the likelihood of contamination by a hotter star is small. It is also unlikely that SPB stars have been misidentified as rotational variables, since SPB stars are rare and do not tend to show the clear harmonic patterns displayed in Fig. 7. The possibility remains that cool background stars or close companions could be responsible for some of the detections of rotational modulation. We do not expect the incidence of such contamination to be a common occurrence in our sample, as cooler background stars would already provide a weaker total light signal than our candidate stars due to cooler surface temperature and greater distance, so that the very small effects (on top of the total light signal) that are due to rotational modulation (see Fig. 8 in this paper) would not always show up as a measurable effect in our candidate stars' light curves. Still, it is possible that the light curves of some of our candidate stars with rotational modulation have been contaminated by distant stars in the same line of sight, or by close companions. We checked contamination flags in the MAST archive and found that most of our 513 candidates had contamination values less than 0.1 and only 29 objects had values between 0.1 and 0.55.
The newly discovered rotational variables are listed in Table 5 (complete table is provided online). Figure 8 shows the distribution of amplitudes of the peak at the rotational frequency. It appears that the number of stars with starspots increases with decreasing amplitude. The number decreases below 0.03 ppt-that could be due to the difficulty in detecting peaks below this amplitude. There is also a very long tail towards the higher amplitudes. In fact, about 4.4% of stars have amplitudes larger than 1 ppt and there are 2 stars where amplitudes are higher than 10 ppt. Such large amplitudes are seen even in B type stars with spots, for example, HD131120 (Briquet et al. 2001).

The period-color relationship
In Fig. 9 we show the photometric period as a function of (B − V ) 0 for all the stars that we classified as spotted stars. To transform (g − r) 0 to (B − V ) 0 we used the equation (Jester et al. 2005): The figure clearly shows a correlation between rotational period and color index. There is a clear trend for hotter stars (especially for (B − V ) 0 between 0.2 and 0.4) to have shorter rotation periods. This is in agreement with many past studies (see, e.g., Bouvier (2013) for a review). According to Barnes (2007), the period P is related to the color (B − V ) 0 Fig. 8 Distribution of amplitudes of stars with a low frequency peak and harmonics. Only amplitudes smaller than 1.5 ppt are shown; the distribution has a long tail where a, b, c, and n are called the gyrochronology parameters. For n, a value of 0.5 is usually assumed, following the findings of Skumanich (1972). The most recent determinations for the other gyrochronology parameters are given by Meibom et al. (2009), who found a = 0.770 ± 0.014, b = 0.553 ± 0.052, and c = 0.472 ± 0.027. We will take the expression  Barnes and Kim 2010). The problem with this function is that it is not defined for (B − V ) 0 < 0.47 and therefore does not adequately describe the observations, as there is an unambiguous extension of a strong correlation of rotation period with color index at least down to (B − V ) 0 = 0.17. Our plot in Fig. 9 is modelled on Fig. 6 in Reinhold and Gizon (2015) (RG15 in what follows). In that paper, the authors work with a value of c (see above) equal to 0.4, quoting Barnes (2007). The asymptote at c = 0.4 is also clearly indicated in Fig. 6 of RG15. Barnes (2007) actually uses c = 0.4 as well as c = 0.5 in his paper, according to which, the sample of stars are being fitted. For reasons stated in Sect. 2 of this paper, we have limited our study to stars with rotation periods below 10 days. That difference notwithstanding, our Fig. 9 shows both an important agreement and an important difference with Fig. 6 in RG15: We see the same sparse population of stars with shorter periods at high values of (B − V ) 0 as they do; however, the clear and sharp break to the left of the (B − V ) 0 = 0.47 line in RG15 is not seen at all in our data. This is an intriguing finding and needs further exploration. In order to examine the nature of this observed extension of the canonical rotation period-color relation, we adapted Eq.  Fig. 9. Rather, the stars seem to possess rotation speeds evenly scattered between 0.3 d and 10 d (our upper limit). Admittedly, only few dozen stars appear blueward of (B − V ) 0 = 0.17 in our figure, too small a number to make a definitive statement.
We wish to emphasize the conservative approach we have taken in selecting candidate stars showing rotational modulation-as discussed earlier in this section. This gives us confidence in the reality of the break around (B − V ) 0 = 0.17. The most intriguing, and potentially most important feature of the distribution of the 513 stars plotted in Fig. 9, is that the relation between rotation period and intrinsic color is neither linear nor random, but appears to follow the same Skumanich-type law that has been associated with the stars redward of (B − V ) 0 = 0.47 in the past (as discussed in detail in Reinhold and Gizon (2015)). This relation is clearly visible even in the absence of the exploratory curves plotted onto Fig. 9. The most reasonable conclusion we can draw from this, is that our results, in agreement with the results independently presented in Balona et al. (2013a), indicate an extension of the canonical relation between rotation period and intrinsic color, blueward of the canonical cut-off. It is interesting to note that, even in RG15, as well as in Fig. 8 of RRB13, a clear relation between rotation period and color is also seen for a population of stars blueward of (B − V ) 0 = 0.4: In RRB13 this population extends to (B − V ) 0 = 0.3 and in RG15 it extends to (B − V ) 0 = 0.2. The large scatter of these populations for any fixed value of (B − V ) 0 appears more akin to the classical Skumanich age-dependent scatter than to a simple decline of surface rotation speed with decreasing mass for stars blueward of F5 (the canonical view). We have now detected a similar mode of behavior blueward of (B − V ) 0 = 0.4. Note that the stars plotted in Fig. 9 are limited to the new rotational variables we discovered in the Kepler sample, which is why we do not reproduce the large populations of points redward of (B − V ) 0 seen in RRB13 and RG15.
The very clear presence of a correlation between rotation period and color index down to (B − V ) 0 = 0.17 is a much stronger confirmation of the trend already noted in Balona et al. (2013a,b) for the clusters NGC 6866 and NGC 6819. It was assumed that the low-frequency variation seen in the Kepler observations of these stars is the rotational frequency, as found in cooler stars. This fact is supported by the agreement between the distribution of equatorial velocities resulting from the assumption and the distribution of equatorial velocities derived from line broadening (Balona 2011). The trend displayed in Balona et al. (2013a,b), where the period-color correlation extends to (B − V ) 0 = 0.17, is substantially strengthened by the results we show in Fig. 9. This clear correlation needs explanation. Since there is less brightness difference between spot and photosphere for hotter stars, leading to lower photometric amplitudes, it is very difficult to observe rotational modulation in stars earlier than mid-F type, from the ground. Hence, for normal A type stars with (B − V ) 0 < 0.3, rotational periods derived from rotational modulation are not available in the literature.
In fact, the empirical function used by Barnes (2007) is undefined for (B − V ) 0 < 0.40 and the rotational evolution model developed by Barnes and Kim (2010) is not applicable for these stars. Our study as well as some earlier studies suggests that there might be no definitive break in the relationship at the boundary between convective and radiative photospheres as previously suggested by Barnes (2003). We conclude that some additional mechanisms must also play a role in the PtM relationship. Hence, we reiterate the observation by Balona et al. (2013a), that a re-evaluation of the physical cause of the PtM relationship (or a re-evaluation of the presence of convection in A-type stars) is required.

Conclusion
The time-series photometry from the Kepler archive at MAST was used to investigate solar-like oscillations, pulsations and rotational behavior of A-K stars in the Kepler field. This study builds on prior work done on solar-like oscillations in red giants in open clusters and in the Kepler field in general (Stello et al. 2009bBasu et al. 2010;Huber et al. 2014). We independently identified solar-like oscillations in 23 new red giants. The basic astrophysical parameters such as mass, radius and luminosity for these stars were determined using their pulsational properties. On comparison with the luminosities and radius from GAIA DR2, we found that the GAIA-based results agree well with our estimated values from the solar-like oscillations.
Observations from space are changing the outlook of γ Dor and δ Sct variables. In the present study, a global analysis of a sample of 15106 stars has been performed to search for non-radial pulsators and hybrids. A careful seismic study of individual stars is needed to confirm their classification and fully characterize their properties and understand the relationship between them. In this regard, our list of candidate non-radial pulsators will serve as an useful resource, for any further detailed analysis of individual objects. The γ Dor stars are difficult to distinguish from rotational variability caused by starspots, but in some stars the periodogram morphology can help to distinguish between the two scenarios. By carefully looking for the presence or absence of harmonics, we attributed our targets respectively to the class of rotationally variable stars (those stars showing harmonics of the rotation period in their periodograms) or to the class of the γ Dor stars (those without harmonics).
Until recently, the γ Dor stars were not considered very useful as asteroseismic tools, owing to the unknown effects of rotation on the frequencies and the lack of mode identification. A study by  suggested that pulsation and rotation periods might be very closely related to each other. Recently, Ouazzani et al. (2017) investigated the effect of uniform rotation on the g-mode pulsation spectra, and more specifically, on the g-modes period spacings. It led to the introduction of a new observable Σ. This is the slope of the period spacing when plotted as a function of period. It is based on the non-uniformity of the period spacings of γ Dor stars which is related to the internal rotation. They further studied the relationship between this new observable and different stellar parameters, such as metallicity, centrifugal distortion and type of mixing. Many of the newly discovered γ Dor stars in our study can be suitable targets for a similar analysis, but also including differential rotation. Moreover, with missions such as TESS and PLATO, γ Dor variables will provide a useful testing bed for this new observable and asteroseismology in general, because it needs to be further investigated to what extent the differential rotation influences the excitation of the observed modes.
In this study, several δ Sct pulsators were identified and further study on identification of modes by means of groundbased multi-color photometry will be vital to reveal the potential for asteroseismology in these stars. To firmly characterize the pulsation properties of the candidate γ Dor/δ Sct hybrid stars, a more detailed analysis, also involving a spectroscopic study, is required. We studied the two new empirical observables of Energy and Efficiency, first introduced by Uytterhoeven et al. (2011), for a larger sample of non-radial pulsators. The peak values of the distributions in Energy, Efficiency and Amplitude for δ Sct, γ Dor/δ Sct hybrid and γ Dor variables are presented according to the gaussian fits to the histograms. The peaks vary in the log(En) distribution by more than two orders of magnitude, where the δ Sct stars have the higher values. Due to this, we suggest that the distribution in the logarithm of Energy (log(En)) can be used as a potential tool to distinguish the non-radial pulsators, to some extent, along with some other classification scheme. Further investigations are required to understand the physical basis of these empirical observables and its impact on the relationship between δ Sct, γ Dor and hybrid variables.
A potentially very significant result of the present work is the detection of a correlation of rotation period with color index of Kepler field stars down to (B − V ) 0 = 0.2 or even lower. This correlation is based on the measurement of rotation periods for 513 stars which are likely to be rotational variables owing to the presence of harmonics of the dominant frequency in the periodogram (or peak periods too long to be attributed to γ Dor pulsations). There is no known explanation for this observation. We did detect the expected correlation between rotational period and color index, which is explained as the result of magnetic braking of rotation via the interaction of convective envelopes with stellar rotation in cooler stars and which is generally applied to stars with (B − V ) 0 > 0.47 and used in the technique known as 'gyrochronology' (Barnes 2001(Barnes , 2003(Barnes , 2007Barnes and Kim 2010). We were surprised, however, to find an additional correlation for hotter stars, mimicking the established relation to some extent.
It is believed that for stars hotter than mid-F, rotational braking is negligible [hence, the limit of applicability of the PtM relation to (B − V ) 0 > 0.47]. If the cause of the PtM relationship is rotational braking due to mass-loss arising from convection, one is led to conclude that mass-loss and convection are important in A and early F-type stars as well. However, the convective zones in the envelopes of these hot stars are believed to be very thin. Hence, we conclude that either the role of convection in A-type stars is not fully understood or that something else is responsible for, at least, a period-color relationship very similar to the standard PtM relationship, but for hotter stars.