FliPer: A global measure of power density to estimate surface gravities of mainsequence solarlike stars and red giants
^{1}
IRFU, CEA, Université ParisSaclay, 91191 GifsurYvette, France
email: lisa.bugnet@cea.fr
^{2}
AIM, CEA, CNRS, Université ParisSaclay, Université Paris Diderot, Sorbonne Paris Cité, 91191 GifsurYvette, France
^{3}
School of Physics and Astronomy, University of Birmingham, , Edgbaston, Birmingham B15 2TT, UK
^{4}
Stellar Astrophysics Centre, Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, 8000 Aarhus C, Denmark
^{5}
Instituto de Astrofísica de Canarias, 38200 La Laguna, Tenerife, Spain
^{6}
Universidad de La Laguna, Dpto. de Astrofísica, 38205 La Laguna, Tenerife, Spain
^{7}
Space Science Institute, 4750 Walnut Street Suite 205, Boulder, CO 80301, USA
^{8}
INAF – Osservatorio Astrofisico di Catania, Via S. Sofia 78, 95123 Catania, Italy
Received:
26
March
2018
Accepted:
13
September
2018
Asteroseismology provides global stellar parameters such as masses, radii, or surface gravities using mean global seismic parameters and effective temperature for thousands of lowmass stars (0.8 M_{⊙} < M < 3 M_{⊙}). This methodology has been successfully applied to stars in which acoustic modes excited by turbulent convection are measured. Other methods such as the Flicker technique can also be used to determine stellar surface gravities, but only works for log g above 2.5 dex. In this work, we present a new metric called FliPer (Flicker in spectral power density, in opposition to the standard Flicker measurement which is computed in the time domain); it is able to extend the range for which reliable surface gravities can be obtained (0.1 < log g < 4.6 dex) without performing any seismic analysis for stars brighter than Kp < 14. FliPer takes into account the average variability of a star measured in the power density spectrum in a given range of frequencies. However, FliPer values calculated on several ranges of frequency are required to better characterize a star. Using a large set of asteroseismic targets it is possible to calibrate the behavior of surface gravity with FliPer through machine learning. This calibration made with a random forest regressor covers a wide range of surface gravities from mainsequence stars to subgiants and red giants, with very small uncertainties from 0.04 to 0.1 dex. FliPer values can be inserted in automatic global seismic pipelines to either give an estimation of the stellar surface gravity or to assess the quality of the seismic results by detecting any outliers in the obtained ν_{max} values. FliPer also constrains the surface gravities of mainsequence dwarfs using only longcadence data for which the Nyquist frequency is too low to measure the acousticmode properties.
Key words: asteroseismology / methods: data analysis / stars: oscillations
© ESO 2018
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
The precise knowledge of stellar parameters is crucial for a very broad range of fields in astrophysics. It helps us understand stellar evolution; this knowledge also provides important information needed for planetary searches and for studying the chemical and dynamical evolution of our Galaxy. In the last decade, the NASA mission Kepler (Borucki et al. 2013) continuously collected very highquality photometric data for almost 200 000 stars (Mathur et al. 2017) over ∼4 years. These observations not only revolutionized the search for exoplanets, but also opened a window into stellar physics. Asteroseismology proved to be a very powerful tool able to better characterize the stars in terms of mass, radii, and age (Metcalfe et al. 2010; Mathur et al. 2012; Silva Aguirre et al. 2017; Serenelli et al. 2017), and also in terms of their rotation and magnetic activity (McQuillan et al. 2014; García et al. 2014a; Davies et al. 2015; Ceillier et al. 2017; Kiefer et al. 2017). However, stellar oscillations have not been detected in all red giant (RG) stars (∼16 000 reported out of the ∼24 000 in the latest Kepler starproperties catalog Mathur et al. 2017) or in all the mainsequence (MS) solarlike stars. Around 135 000 MS dwarfs have only been observed in long cadence (LC) with a sampling time of 29.4 min, by Kepler preventing any direct asteroseismic analyses because their acousticmode frequencies are well above the Nyquist frequency and can only be seismically studied with shortcadence (SC) data with a sampling time of 58.85 s (Chaplin et al. 2011b).
To circumvent this, new techniques are being developed to extract precise surface gravities (log g) directly from the photometric data. This is the case of the Flicker method, i.e., the measurement of the brightness variations on timescales shorter than 8 h (Bastien et al. 2013, 2016), the variance in the flux (Hekker et al. 2012), and the granulation (Mathur et al. 2011; Kallinger et al. 2014), and the analysis of the timescales of convectivedriven brightness variations (Kallinger et al. 2016). However all these techniques have limitations. Flicker is restricted by construction to stars with 4500 < T_{eff} < 7150 K and 2.5 < log g < 4.6dex, preventing the study of highluminosity redgiant branch (RGB) and asymptotic redgiant branch (AGB) stars. To obtain the granulation properties it is necessary to fit a complicated model including different scales of convection with many free parameters (for more details see the discussions in Mathur et al. 2011; Kallinger et al. 2014; Corsaro et al. 2017). The final method requires that the oscillation signal be temporally resolved, which prevents extending the analysis to MS dwarfs only observed in LC data. It has also been shown that instead of using classic seismic methods, it is possible to apply machine learning algorithms directly on the data. For instance, Hon et al. (2018a) apply a convolutional neural network on spectra to classify stars. This method gives good results for about 99% of their sample of RGs, including some stars that were not already characterized with seismic pipelines. A random forest regression model (Breiman 2001) applied directly on the photometric light curves of variable stellar sources can also estimate their surface gravity with a 0.42 dex uncertainty (Miller et al. 2014).
We present here a new metric called FliPer (Flicker in power, in opposition to the standard Flicker measurement which is computed in the time domain), which links the variability of a star to its surface gravity in a wider range than Flicker, starting at a log g ∼ 0.1 and similar effective temperatures (4500 < T_{eff} < 7150 K) covering solarlike pulsating stars. We are limited in the 0.1 < log g < 4.6 dex range of surface gravity because of the lack of information we have on extreme surface gravity solarlike stars. There is no intrinsic limits of applicability to the FliPer calculation. We decide to combine powerful methods: we include FliPer values from different lower frequency boundaries into a supervised machine learning random forest algorithm in order to get even more accurate results on the surface gravity estimation. This way, we obtain information about the impact of the lower frequency boundaries and the effective temperature on the estimation of surface gravity.
2. Observations, data selection, and preparation
In this work, LC data (Gilliland et al. 2010) obtained by NASA’s Kepler main mission are used. The light curves have been corrected and the different quarters concatenated following García et al. (2011). Two highpass filters have been used with cutoff frequencies corresponding to 20 and 80 days. To minimize the effects of the gaps in the observations (García et al. 2014b) the missing observations have been interpolated using inpainting techniques (Pires et al. 2015). The power spectrum density is then computed for each star (calibrated as a singlesided spectrum). Data are corrected for apodization following Chaplin et al. (2011a).
We selected ∼15 000 RG stars among the ones in Mathur et al. (2017) showing stellar pulsations and characterized using the A2Z asteroseismic pipeline (Mathur et al. 2010). These stars have 0.1 < log g < 3.4 dex and 3285 K < T_{eff} < 7411 K. In addition, 254 MS stars with 4951 K < T_{eff} < 6881 K are used to extend the study towards a higher surface gravity range, reaching 4.5 dex. These stars have Kepler magnitudes brighter than 14 (Kp < 14).
It is important to note that the values of ν_{max} computed by A2Z do not show any systematic biases at a level of ∼1% when compared to other seismic pipelines as shown by Pinsonneault et al. (2018).
3. The new metric: FliPer
The complete power spectrum contains contributions from the stellar variability at all timescales, such as oscillation modes, surface granulation, and rotation. We define FliPer as
where represents the averaged value of the power spectrum density from a given frequency (see Sect. 3.1) to the Nyquist frequency and P_{n} is the photon noise. This noise can be calculated by taking the average value of the PSD over a range of frequencies close to the Nyquist frequency, but this method leads to biased estimation of FliPer for stars that oscillate with a frequency close to the Nyquist frequency, as explained in detail by Bugnet et al. (2017). Then the photon noise is computed following the empirical expression obtained by Jenkins et al. (2010).
The value of FliPer is dominated by a combination of the granulation and the oscillation modes, which both depend on the evolutionary stage of the star. The more evolved the star, the larger their oscillation and granulation amplitudes (e.g., Mosser et al. 2012; García & Stello 2015), while the frequency of maximum power ν_{max} decreases (e.g., Bedding 2014).
It is important to note that the signature of strong rotation (and its harmonics) would bias FliPer. This does not have a large impact in the case of RGs because a very small fraction of them show signatures of the rotation in the PSD, as shown by Ceillier et al. (2017), but needs to be studied in detail for MS solarlike stars (see Sect. 3.3).
3.1. Computing FliPer from data
The observational frequency range used to compute is limited at high frequency by the Nyquist frequency. For most stars (those observed in LC) we cannot get information above ∼283 μHz. Therefore, we selected a first set of calibrator stars including a RG pulsating at a frequency lower than 300 μHz and for which asteroseismic parameters were available. A second set of known seismic MS dwarfs was used to study FliPer with LC data only.
The lowfrequency limit of is given by the cutoff frequency used in the calibration of the data. For most of the stars, a 20day highpass filter light curve is used. The associated cutoff frequency of the signal is 0.58 μHz. We thus establish a lowfrequency limit for the analysis at 0.7 μHz. As MS stars can rotate with a period shorter than 20 days, FliPer is computed with a lowfrequency limit at 7 μHz (i.e., ∼1.6 days) avoiding most of the pollution induced by rotation signals (see Sect. 3.3). For stars showing rotation harmonics at higher frequencies, the lowfrequency boundary should be taken even higher (e.g., 20 μHz) to avoid any additional impact on FliPer from the peaks associated with rotation. Finally, a small number of RGs in our sample are either highluminosity RGB or AGB stars (log g < 1.2 dex) pulsating at frequencies smaller than the 20day cutoff frequency of the calibrated data. For these stars, an 80day filter is used in the calibration process. It allows us to properly measure the stellar signal down to 0.2 μHz (which is the limit frequency utilized in this analysis) and to include oscillationmode power into the FliPer value.
3.2. First surface gravity estimator
For stars with solarlike oscillations, seismic surface gravities are directly obtained from the frequency of maximum oscillation power ν_{max} computed with the A2Z pipeline (Mathur et al. 2010) and effective temperatures from the kepler DR25 catalog (Mathur et al. 2017). Knowing seismic surface gravities with their uncertainties allows us to study the behavior of FliPer with the evolutionary state of the stars using only LC light curves even for MS stars. It is important to note that for MS stars the seismic log g has been seismically inferred using SC data although FliPer has been computed using LC data.
In Fig. 1 the seismic log g is represented as a function of FliPer (black dots). Three different areas have been identified depending on the evolutionary state of the star: MS and subgiant stars (gray shaded region), RG (red), and highluminosity RG stars from the branch and from the asymptotic branch (blue). For each of these star categories FliPer was computed with a different lowfrequency limit of 7 μHz (avoiding in most cases the region of possible pollution by rotation signatures present on data filtered with a highpass filter at 20 μHz), 0.7 μHz (20day filter), and 0.2 μHz (80day filter), respectively. The color scheme is the same in the captions.
Fig. 1. Seismic log g VS FliPer for ∼15 000 stars observed with Kepler longcadence observational mode (black dots). The gray, red, and blue shaded regions correspond to MS and subgiant stars, RG, and highluminosity RG stars or AGB, respectively. White diamonds represent the weighted mean value of log g (see Eq. (2)) computed with 30 MS stars and 300 RGs each. The white area delimits the location where 68% of the stars in our sample are around the mean value. The white boundary represents the equivalent of a 1σ uncertainty (standard deviation) on the surface gravity obtained from FliPer. 

Open with DEXTER 
In order to characterize the relationship between FliPer and log g represented in Fig. 1, we calculate an average value of log g for each bin of n stars (n = 300 for RG and n = 30 for MS and subgiant stars) as
where log g_{i} represents the value of surface gravity for each star and δlog g_{i} the corresponding uncertainty.
These values are represented by the white diamonds, and are located at the averaged value of FliPer over each bin. To define the 1σ uncertainties, we compute the area containing 68% of the stars of the sample (marked by a white contour region in Fig. 1). Mean values and their corresponding ±1σ uncertainties are given in Table. A.1. By using these mean values it is possible to estimate the stellar surface gravity directly from the FliPer estimator. The uncertainties obtained on log g extend from 0.05 to 0.2 dex, depending on the evolutionary state of the star. Because of the calculation of the mean values, the boundaries in log g are reduced to the range 4.35–0.38 dex, as shown in Table A.1.
3.3. Disentangling mainsequence stars from red giants
As defined, FliPer is mostly dominated by a combination of the power coming from granulation and oscillation modes (when the latter are below the Nyquist frequency). The limitation in the use of the calibrated values from Table A.1 to directly estimate surface gravity of stars appears when the spectrum shows a specific behavior that strongly modifies the mean value of the power density. For instance, in stars showing large excess of power (e.g., due to spikes at thrusters frequency in K2 data or to pollution from a background binary), the value of FliPer is biased towards high power density (Bugnet et al. 2017). On the contrary, in stars with a low signaltonoise ratio the value of FliPer is biased towards lower values because most of the spectrum is dominated by the instrumental noise. As a consequence, FliPer is higher than expected for fast rotating MS stars due to the rotation peaks and their harmonics, which can be particularly high for young MS dwarfs. For these stars, the log g inferred from Table A.1 could be such that it corresponds to a RG star and not to a MS star, even if we calculate FliPer with the 7 μHz frequency limit.
To avoid this problem and to disentangle any MS stars from RGs, we need an additional parameter that takes into account the power due to rotation. The simplest solution is to combine different FliPer values, including some at higher frequencies than the 7 μHz limit. For each star in our sample, we then calculate FliPer with several lowfrequency limits (e.g., F_{p0.2} from 0.2, F_{p0.7} from 0.7, F_{p7} from 7, F_{p20} from 20, and F_{p50} from 50 μHz). For MS stars with small rotation signatures the value of FliPer is almost the same for all the lowfrequency boundaries (see panel a in Fig. 2). However, when rotation peaks are present, there is a large difference between the FliPer parameters, depending on the frequency of the rotational peaks (see Fig. 2). This is the case for the stars KIC 8298090 and KIC 5357446 represented in panels b and c. In panel b all the rotational components are below the 0.7 μHz boundary, meaning that parameters F_{p20} and F_{p50} were not necessary to classify this star as a MS star. However, in panel c the rotation peaks reach higher frequencies; in order to estimate the surface gravity of this star the two new highfrequency parameters are needed. Panel d shows a RG star for comparison. In the regime of RG stars, all the FliPer values are very similar, except the lowest ones coming from the calculation with the highpass filter that does not include the range of frequency of oscillation modes. By comparing the values of FliPer computed with different lowfrequency limits, it is then possible to disentangle MS stars with a high rotation signature from RG stars. This can be done in a starbystar analysis (see Fig. 2). However, it is possible to automatize this procedure, as is explained in the following section.
Fig. 2. Impact of the lower frequency limit in the FliPer calculation on the estimation of surface gravity for different types of star. Left panels: power density spectra of four Kepler targets. Colored areas (blue, red, black, orange, and green) represent the different ranges of frequency used for FliPer calculation (respectively from 0.2, 0.7, 7, 20, and 50 μHz to the Nyquist frequency). The color scheme is the same in the captions. Right panels: all studied Kepler stars from Fig. 1 are represented in gray in the log g V.S FliPer diagram. Colored stars (blue, red, black, orange, and green) show the positions in the diagram of the four stars from the left panels with their color corresponding to the lowfrequency boundaries used to compute the FliPer value. Panel a represents a MS star without any visible rotation component, panel b a MS star showing rotation, panel c a highfrequency rotating MS star, and panel d a RG star. 

Open with DEXTER 
4. Seismic independent surface gravity prediction from 0.1 to 4.5 dex
The direct estimation of surface gravity from Table A.1 gives good results only when the evolutionary state of the star is already known, and when the spectrum does not show a specific behavior that strongly modifies the mean value of the power density (e.g., when the PSD is polluted by spikes of a background binary or a classical pulsator). The reason is that we only use one value of FliPer computed from one lower frequency limit. Estimating surface gravities of unclassified or complex stars requires a different use of the FliPer method.
4.1. Using machine learning
As explained above, combining different FliPer values is a powerful way to detect MS stars showing high rotation signal among RGs. It also means that by using different highpass filters in the calculation of FliPer we are sensitive to different physical signatures in the PSD. Combining them in the study thus improves the characterization of the star, and we intend to use this wisely to predict surface gravities. To do so, we train a random forest regressor algorithm (Breiman 2001) on a random subsample representing 80% of our set of stars. The random forest method is based on the aggregation of a large number of decision trees (see Appendix B for a description of the method) that has already been proved to be useful in asteroseimology (e.g., Miller et al. 2014). The trees are constructed from a training data set and internally validated to give a prediction based on the predictor for future observations. The random forest method not only allows the use of a large number of parameters but also estimates their individual impact on the regression. The parameters used to estimate surface gravities are F_{p0.2}, F_{p0.7}, F_{p7}, F_{p20}, F_{p50}, T_{eff}, and Kp. They represent the values of FliPer calculated from a lowfrequency limit (0.2, 0.7, 7, 20, 50 μHz), the effective temperature, and the Kepler magnitude of the star.
4.2. Building training and test sets of parameters for the random forest algorithm
We intend to take into account the uncertainties on the parameters during the testing of the algorithm to estimate the intrinsic bias and/or uncertainties of our methodology. The uncertainties on effective temperature δT_{eff} are taken directly from the Mathur et al. (2017) catalog. The error on the surface gravity δlog g comes from the uncertainty on ν_{max} from the seismic analysis of the stars. We can estimate the uncertainty due to the photon noise in the spectra (following a chisquared distribution with two degrees of freedom) impacting the determination of FliPer by considering negligible the uncertainty made on the photon noise. Hence, the uncertainty on FliPer can be explicitly written as
where δP_{i} stands for the error made on the power contained in each bin and N_{bin} is the total number of bins in the power density spectrum. The individual δP_{i} cannot be extracted directly because the noise distribution does not have Gaussian errors. We thus use the central limit theorem and we rebin the spectrum by a factor of n = 50. The total amount of power in the spectrum is then expressed as
where P_{n,j} follows a quasinormal distribution with 2n degrees of freedom. This assumes that the signal does not change dramatically over this range of 50 bins, which is consistent with the shape of the spectra in solarlike stars. The uncertainty on the mean from each n bins is then expressed as
leading to a global uncertainty on FliPer values of
Then we include the effect of these errors on the different parameters on the testing of the algorithm. To do so, we perform a Monte Carlo simulation by generating for each star in our test sample (representing 20% of the total number of stars in our study) 100 artificial sets of parameters from their corresponding normal distributions. With being 100 random values following the standard distribution, we calculate for each X parameter (F_{p0.2}, F_{p0.7}, F_{p7}, F_{p20}, F_{p50},T_{eff}, and log g) 100 new values X_{0} ≤ i ≤ 100 following Eq. (7) below. However, the Kepler magnitude of the star remains constant as it has no uncertainties, and it completes each of the 100 new sets of parameters:
Each new group of 100 sets of parameters is considered in the following study as a representation of a hundred stars to test the algorithm.
4.3. Impacts of parameters on the training
We used 80% of our stars to train the algorithm to estimate the surface gravity. The remaining 20% was used to test the performance of the algorithm by taking into account uncertainties on the different parameters as explained in Sect. 4.2. The impacts of the different parameters on the training process are represented in Fig. 3 A.
Fig. 3. Left panel: importance of the different parameters F_{p0.2}, F_{p0.7}, F_{p7}, F_{p20}, F_{p50},T_{eff}, and Kp in the training process. The color scheme is the same in the captions. Middle panel: histogram of the residuals of the estimated log g_{RF} values from the references log g_{A2Z}. Right panel: histogram of the mean absolute deviation from the expected value MAD(log g_{A2Z} – log g_{RF}). 

Open with DEXTER 
A predictable result is that the F_{p0.7} parameter largely dominates the training. It comes from the fact that this is the most suitable parameter to study RG, representing more than 90% of the total number of stars. Other relevant values of the filtering appear to be 7 and 0.2 μHz. Indeed, F_{p7} plays an important role in the training because of its ability to distinguish MS stars from RGs, and the F_{p0.2} parameter helps in the prediction of surface gravity for highluminosity stars. The other parameters F_{p20} and F_{p50} have lower impacts on the training, but still help the learning for high rotating MS stars. Impacts of the effective temperature and Kp do not exceed a small fraction. We confirm from Fig. 3 that combining different lower frequency boundaries in the FliPer calculation makes a great difference for the estimation of robust surface gravities.
4.4. Results
To evaluate the performance of the algorithm, the estimate of surface gravity from the test sample is compared to the corresponding A2Z estimation of surface gravity. The mean absolute deviation (MAD) of the random forest surface gravities from reference values is reported in Table 1. This estimator of the deviation is chosen to be robust against outliers to avoid any issue coming from an eventual remaining error in the A2Z estimation of surface gravity.
The estimation of surface gravity resulting from the machine learning on the test sample has an averaged deviation of ∼0.046 dex from our reference values (see Table 1). We also obtain surface gravity deviation and errors from the reference values for different ranges of surface gravity in Table 1. We conclude that for all our stars, the new method gives a very good precision on surface gravity; the Flicker method in the range of log g of 2.5–4.6 dex has typical errors between 0.1 and 0.2 dex. Here, errors are in the range from ∼0.04 dex for RGs to ∼0.09 dex for highluminosity stars (see Fig. 3 right panel). Our estimates are on average centered around the log g reference values (see Table 1). There is a small bias for HL and MS stars (see Fig. 3 middle panel) because the estimation of extreme surface gravities is the hardest for the algorithm, which sometimes becomes slightly biased by the presence of many RGs in the sample. This bias, which depends on the evolutionary state of the star, should be taken into account, but it remains smaller than the uncertainties on the original surface gravity values.
Our algorithms are available on GitHub^{1} where the functions to calculate FliPer and the random forest algorithm are provided. We also provide the already trained algorithms for the estimation of surface gravities. They can be directly applied to any solartype star to estimate its surface gravity from 0.1 to 4.6 dex.
5. Discussion and conclusion
In this work we present a new method used to estimate surface gravity of solarlike stars that extracts information from global power in their spectra. The sample of ∼15 000 stars is constituted of MS and subgiant stars, stars on the redgiant branch (RGB) and clump stars, and also highluminosity stars on the asymptotic giant branch (AGB). This way, we study stars with 0.1 < log g < 4.5 dex in which mode oscillations are expected to arise from surface convection. Power spectra should then present patterns of granulation power, rotation components, and oscillationmode power.
FliPer values are calculated by taking the average power density normalized by the photon noise of the star from different lower frequency limits to the Nyquist frequency. Our first method consists of calibrating surface gravity of stars from their FliPer value with a 1σ uncertainty (see Table A.1). We explained how these values can be used directly to give a first estimate of surface gravity; however, it works well only on stars that are already characterized. The evolutionary state has to be known or the star must have a weak rotational signature in order to distinguish MS stars from RGs. To give estimations of surface gravities for any star, we introduce a second method. A random forest regressor algorithm is trained to estimate surface gravity on a sample of our stars. We use FliPer values computed with different frequency ranges, spectroscopic effective temperatures, and seismic surface gravities. This way, stars are better characterized during the process, and no additional information is needed to provide accurate estimations of surface gravity, even for highly rotating MS stars. By testing the algorithm on the rest of our sample, we obtain estimates of surface gravity with a mean absolute deviation of 0.046 dex from seismic log g. The training relies on seismic observations of solarlike stars representing 80% of our sample. However, there is no need for additional seismic measurements to obtain precise estimations of surface gravities on the test set of stars. The uncertainty on our results largely improves upon previous nonseismic estimations of surface gravity. Spectroscopic estimations are known to have 0.1–0.3 dex error bars (Frasca et al. 2016; García Pérez et al. 2016). Recent methods such as Flicker (Bastien et al. 2016) give estimates with errors higher than 0.1 dex, while the study of the granulation timescale (Kallinger et al. 2016) is limited to stars showing a visible oscillation pattern but with better uncertainties, around 0.018 dex. In addition, FliPer is extended to a wider range of surface gravities, reaching log g as small as 0.1 dex with a mean absolute deviation of 0.046 dex comparable to the other RGs.
For MS stars that oscillate at high frequency (above the Kepler LC Nyquist frequency), FliPer computed from LC data does not contain mode power, but only granulationrelated power (Corsaro et al. 2017) and rotation signals. However, Fig. 1 clearly shows that FliPer values for MS stars are still correlated with surface gravity. This is new evidence of the link between granulation and asteroseismic properties (Mathur et al. 2011; Kallinger et al. 2014), allowing us to estimate ν_{max} or rather surface gravities on LC data for which highfrequency modes are not measured. Thus, proper surface gravities can be precisely inferred for any Kepler LC solarlike target, from MS to highluminosity stars, without using direct seismic analysis.
Lots of studies concern the estimation of seismic parameters of stars with new techniques directly from the properties of the time series or the power spectrum density (Bugnet et al. 2017; Hon et al. 2018b; Pande et al. 2018; Bell et al. 2018). We thus adapt our methodology to estimate ν_{max} instead of the surface gravity based on the same sample of stars (see Appendix C). The results are of course very similar to those on the surface gravity, with uncertainties on ν_{max} about 0.044 dex and a mean distance to references ν_{max} values (∼1.3 × 10^{−3} dex) negligible. Moreover, FliPer has already been included by Bell et al. (2018) as a validation procedure to their seismic results and it is also being implemented as one of the parameters to be used in the classification algorithm that is being developed for the NASA TESS mission (Ricker et al. 2014) using a random forest classifier (Tkachenko et al., in prep.).
Acknowledgments
We thank the anonymous referee for the very useful comments. We would also like to thank M. Benbakoura for his help in analyzing the light curves of several binary systems included in our set of stars. L.B. and R.A.G. acknowledge the support from PLATO and GOLF CNES grants. S.M. acknowledges support from the National Aeronautics and Space Administration under Grant NNX15AF13G, the National Science Foundation grant AST1411685, and the Ramon y Cajal fellowship no. RYC201517697. E.C. is funded by the European Union’s Horizon 2020 research and innovation program under the Marie SklodowskaCurie grant agreement no. 664931. O.J.H and B.M.R. acknowledge the support of the UK Science and Technology Facilities Council (STFC). Funding for the Stellar Astrophysics Centre is provided by the Danish National Research Foundation (Grant DNRF106). This research has made use of NASA’s Astrophysics Data System. Data presented in this paper were obtained from the Mikulski Archive for Space Telescopes (MAST). STScI is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS526555.
References
 Bastien, F. A., Stassun, K. G., Basri, G., & Pepper, J. 2013, Nature, 500, 427 [NASA ADS] [CrossRef] [Google Scholar]
 Bastien, F. A., Stassun, K. G., Basri, G., & Pepper, J. 2016, ApJ, 818, 43 [NASA ADS] [CrossRef] [Google Scholar]
 Bedding, T. R. 2014, Solarlike oscillations: An observational perspective, eds. P. L. Pallé, & C. Esteban, 60 [Google Scholar]
 Bell, K. J., Hekker, S., & Kuszlewicz, J. S. 2018, MNRAS, accepted [arXiv:1809.09135] [Google Scholar]
 Borucki, W. J., Agol, E., Fressin, F., et al. 2013, Science, 340, 587 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Breiman, L. 2001, Mach. Learn., 45, 5 [CrossRef] [Google Scholar]
 Bugnet, L., Garcia, R. A., Davies, G. R., Mathur, S., & Corsaro, E. 2017, Proceedings SF2A 2017 [Google Scholar]
 Ceillier, T., Tayar, J., Mathur, S., et al. 2017, A&A, 605, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chaplin, W. J., Kjeldsen, H., Bedding, T. R., et al. 2011a, ApJ, 732, 54 [NASA ADS] [CrossRef] [Google Scholar]
 Chaplin, W. J., Kjeldsen, H., ChristensenDalsgaard, J., et al. 2011b, Science, 332, 213 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Corsaro, E., Mathur, S., García, R. A., et al. 2017, A&A, 605, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Davies, G. R., Chaplin, W. J., Farr, W. M., et al. 2015, MNRAS, 446, 2959 [NASA ADS] [CrossRef] [Google Scholar]
 Frasca, A., MolendaŻakowicz, J., De Cat, P., et al. 2016, A&A, 594, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 García, R. A., & Stello, D. 2015, Asteroseismology of Red Giant Stars in Extraterrestrial Seismology (Cambridge: Cambridge University Press) [Google Scholar]
 García, R. A., Hekker, S., Stello, D., et al. 2011, MNRAS, 414, L6 [NASA ADS] [CrossRef] [Google Scholar]
 García, R. A., Ceillier, T., Salabert, D., et al. 2014a, A&A, 572, A34 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 García, R. A., Mathur, S., Pires, S., et al. 2014b, A&A, 568, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 García Pérez, A. E., Allende Prieto, C., Holtzman, J. A., et al. 2016, AJ, 151, 144 [NASA ADS] [CrossRef] [Google Scholar]
 Gilliland, R. L., Jenkins, J. M., Borucki, W. J., et al. 2010, ApJ, 713, L160 [NASA ADS] [CrossRef] [Google Scholar]
 Hekker, S., Elsworth, Y., Mosser, B., et al. 2012, A&A, 544, A90 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hon, M., Stello, D., & Yu, J. 2018a, MNRAS, 475, 3233 [NASA ADS] [CrossRef] [Google Scholar]
 Hon, M., Stello, D., & Zinn, J. C. 2018b, ApJ, 859, 64 [NASA ADS] [CrossRef] [Google Scholar]
 Jenkins, J. M., Caldwell, D. A., Chandrasekaran, H., et al. 2010, ApJ, 713, L120 [NASA ADS] [CrossRef] [Google Scholar]
 Kallinger, T., De Ridder, J., Hekker, S., et al. 2014, A&A, 570, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kallinger, T., Hekker, S., Garcia, R. A., Huber, D., & Matthews, J. M. 2016, Sci. Adv., 2, 1500654 [NASA ADS] [CrossRef] [Google Scholar]
 Kiefer, R., Schad, A., Davies, G., & Roth, M. 2017, A&A, 598, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kotsiantis, S. B. 2007, Proceedings of the 2007 Conference on Emerging Artificial Intelligence Applications in Computer Engineering: Real Word AI Systems with Applications in eHealth, HCI, Information Retrieval and Pervasive Technologies (Amsterdam, The Netherlands: IOS Press), 3 [Google Scholar]
 Mathur, S., García, R. A., Régulo, C., et al. 2010, A&A, 511, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mathur, S., Hekker, S., Trampedach, R., et al. 2011, ApJ, 741, 119 [NASA ADS] [CrossRef] [Google Scholar]
 Mathur, S., Metcalfe, T. S., Woitaszek, M., et al. 2012, ApJ, 749, 152 [NASA ADS] [CrossRef] [Google Scholar]
 Mathur, S., Huber, D., Batalha, N. M., et al. 2017, ApJS, 229, 30 [NASA ADS] [CrossRef] [Google Scholar]
 McQuillan, A., Mazeh, T., & Aigrain, S. 2014, ApJS, 211, 24 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Metcalfe, T. S., Monteiro, M. J. P. F. G., Thompson, M. J., et al. 2010, ApJ, 723, 1583 [NASA ADS] [CrossRef] [Google Scholar]
 Miller, A., Richards, J., Bloom, J. S., & a larger Team 2014, Amer. Astron. Soc. Meet. Abstr., 223, 125.01 [NASA ADS] [Google Scholar]
 Mosser, B., Elsworth, Y., Hekker, S., et al. 2012, A&A, 537, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pande, D., Bedding, T. R., Huber, D., & Kjeldsen, H. 2018, MNRAS, 480, 467 [NASA ADS] [CrossRef] [Google Scholar]
 Pedregosa, F., Varoquaux, G., Gramfort, A., et al. 2011, J. Mach. Learn. Res., 12, 2825 [Google Scholar]
 Pinsonneault, M. H., Elsworth, Y. P., Tayar, J., et al. 2018, ApJS, submitted [arXiv:1804.09983] [Google Scholar]
 Pires, S., Mathur, S., García, R. A., et al. 2015, A&A, 574, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ricker, G. R., Winn, J. N., & Vanderspek, R. 2014, SPIE Conf. Ser., 9143, 20 [Google Scholar]
 Serenelli, A., Johnson, J., Huber, D., et al. 2017, AJSS, 233, 23 [NASA ADS] [CrossRef] [Google Scholar]
 Silva Aguirre, V., Lund, M. N., Antia, H. M., et al. 2017, ApJ, 835, 173 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Surface gravity as a function of FliPer
Weighted mean value of log g (diamonds in Fig. 1) with their 1σ uncertainties for each bin of 30 (for MS and HL stars) or 300 (for RGs) stars.
Appendix B: Random forest regressor
B.1. Supervised machine learning
A random forest algorithm is a supervised machine learning (ML) method (Kotsiantis 2007). It learns how to predict an output variable (Y_{predicted}) from some training data (X) for which the corresponding result (Y_{known}) is already known. It thus learns a mapping function f from the input(s) to the output:
The algorithm iteratively makes predictions (Y_{predicted}) on the training data (X). These predicted values are corrected to achieve a maximum level of performance by comparing them with the Y_{known} values. We use a supervised ML algorithm for our study because we have input variables X (which are F_{p0.2}, F_{p0.7}, F_{p7}, F_{p20}, F_{p50}, T_{eff}, and Kp) and an output variable Y_{known} (representing our surface gravity log g).
B.2. Regression trees
The regression tree method is part of the Classification and Regression Trees (CART) introduced by Breiman (2001). A decision tree algorithm constructs a binary tree during the training, with each node representing a split point on a single input variable (X) (a numerical value for regression algorithms, or a class name for classification algorithms). The leaf nodes of the tree contain the output possible predictions (Y_{predicted}), as shown in Fig B.1.
Fig. B.1. Schematic representation of a regression decision tree. At each node (green points) one variable is split at a value such that the cost function (see Eq (B.2)) is minimized. 

Open with DEXTER 
The tree is built in such way that the cost function is minimized. Equation (B.2) is the cost function used for the regressor, with N_{training} being the number of stars in our training sample:
Once the tree is built on the training sample, it is used to evaluate Y_{predicted} for new X_{new} data.
B.3. Ensemble method random forest regressor
An ensemble method combines the prediction from multiple ML algorithms together. It aims to make even more accurate predictions than any individual model. The random forest regressor is an ensemble method that combines regression trees. It consists in the following:
Creating many subsamples of the training sample;
Training a regression tree on each subsample, keeping a low number of variables that can be looked at for each split point. It aims to decrease the correlation between the different trees. For a regression algorithm, the typical number of features that can be searched is with p the number of input (X) variables;
Calculating the average prediction from each model for the new test sample: this averaged value is taken as the estimate for the output variable (Y_{predicted}).
In our work we use the “RandomForestRegressor” function from the “sklearn.ensemble” Python library (Pedregosa et al. 2011) to compute the training on surface gravity.
Appendix C: Automatic estimation of ν_{max}
As a complementary study we also trained our algorithm to estimate the frequency of maximum power ν_{max}. The training is made following Sect. 4 by using ν_{max} instead of log g as the predicted parameter Y_{predicted}. The training input values are computed as in Sect. 4.2 by combining the ν_{max} values estimated by the A2Z global seismic pipeline for our sample of ∼15 000 stars along with their uncertainties.
Results are very similar to the estimation of surface gravity, and are given in Table C.1. The estimation of ν_{max} can be made for any star with solarlike oscillations with 0.1 < log g < 3.4 dex, 3285 < T_{eff} < 7411 K, Kp < 14, and provide a very good prior for any more complex seismic analysis of the star. The complete algorithm for the ν_{max} estimation can be found on the Git repository^{1}.
Summary of statistical results on ν_{max} on the test set from Fig. 3.
All Tables
Weighted mean value of log g (diamonds in Fig. 1) with their 1σ uncertainties for each bin of 30 (for MS and HL stars) or 300 (for RGs) stars.
All Figures
Fig. 1. Seismic log g VS FliPer for ∼15 000 stars observed with Kepler longcadence observational mode (black dots). The gray, red, and blue shaded regions correspond to MS and subgiant stars, RG, and highluminosity RG stars or AGB, respectively. White diamonds represent the weighted mean value of log g (see Eq. (2)) computed with 30 MS stars and 300 RGs each. The white area delimits the location where 68% of the stars in our sample are around the mean value. The white boundary represents the equivalent of a 1σ uncertainty (standard deviation) on the surface gravity obtained from FliPer. 

Open with DEXTER  
In the text 
Fig. 2. Impact of the lower frequency limit in the FliPer calculation on the estimation of surface gravity for different types of star. Left panels: power density spectra of four Kepler targets. Colored areas (blue, red, black, orange, and green) represent the different ranges of frequency used for FliPer calculation (respectively from 0.2, 0.7, 7, 20, and 50 μHz to the Nyquist frequency). The color scheme is the same in the captions. Right panels: all studied Kepler stars from Fig. 1 are represented in gray in the log g V.S FliPer diagram. Colored stars (blue, red, black, orange, and green) show the positions in the diagram of the four stars from the left panels with their color corresponding to the lowfrequency boundaries used to compute the FliPer value. Panel a represents a MS star without any visible rotation component, panel b a MS star showing rotation, panel c a highfrequency rotating MS star, and panel d a RG star. 

Open with DEXTER  
In the text 
Fig. 3. Left panel: importance of the different parameters F_{p0.2}, F_{p0.7}, F_{p7}, F_{p20}, F_{p50},T_{eff}, and Kp in the training process. The color scheme is the same in the captions. Middle panel: histogram of the residuals of the estimated log g_{RF} values from the references log g_{A2Z}. Right panel: histogram of the mean absolute deviation from the expected value MAD(log g_{A2Z} – log g_{RF}). 

Open with DEXTER  
In the text 
Fig. B.1. Schematic representation of a regression decision tree. At each node (green points) one variable is split at a value such that the cost function (see Eq (B.2)) is minimized. 

Open with DEXTER  
In the text 