Systematic search for stellar pulsators in the eclipsing binaries observed by Kepler

Eclipsing binaries (EBs) are unique targets for measuring precise stellar properties and can be used to constrain stellar evolution models. In particular, it is possible to measure masses and radii of both components of a double-lined spectroscopic EB at the percent level. Since the advent of high-precision photometric space missions (MOST, CoRoT, Kepler, BRITE, TESS), the use of stellar pulsation properties to infer stellar interiors and dynamics constitutes a revolution for studies of low-mass stars. The Kepler mission has led to the discovery of thousands of classical pulsators such as δ Scuti and solar-like oscillators (main sequence and evolved), but also almost 3000 EBs with orbital periods shorter than 1100 days. We report the first systematic search for stellar pulsators in the entire Kepler EB catalog. The focus is mainly aimed at discovering δ Scuti, γ Doradus, red giant, and tidally excited pulsators. We developed a data inspection tool (DIT) that automatically produces a series of plots from the Kepler light curves that allows us to visually identify whether stellar oscillations are present in a given time series. We applied the DIT to the whole Kepler EB database and identified 303 systems whose light curves display oscillations, including 163 new discoveries. A total of 149 stars are flagged as δ Scuti (100 from this paper), 115 as γ Doradus (69 new), 85 as red giants (27 new), and 59 as tidally excited oscillators (29 new). There is some overlap among these groups, as some display several types of oscillations. Despite the likelihood that many of these systems are false positives, for example, when an EB light curve is blended with a pulsator, this catalog gathers a vast sample of systems that are valuable for a better understanding of stellar evolution.


Introduction
It is widely agreed that asteroseismology has become the most reliable way to infer global and internal properties of solar-like stars, from the main sequence (MS) to red-giant (RG) phase, ever since the remarkable success of the space-borne photometers CoRoT, Kepler, and TESS (Baglin et al. 2009;Borucki et al. 2010;Ricker et al. 2015). The simplest and most popular application of asteroseismology consists in comparing the oscillation global properties of a given star to those of the Sun, and retrieving its mass and radius from the asteroseismic scaling relations (Kjeldsen & Bedding 1995). Since masses and radii lead to ages and distances, solar-like oscillators represent key tracers of the properties of stellar populations. This approach has been used to characterize large samples of stars observed by CoRoT and Kepler, which detected thousands of solar-like oscillation spectra for MS and RG stars (Chaplin & Miglio 2013;Chiappini et al. 2015). It will also play a central role in the ESA PLATO mission, for which asteroseismic inference is expected to constrain stellar and planetary properties of tens of thousands of systems.
Besides, classical pulsators such as γ Doradus (γ Dor), δ Scuti (δ Sct), β Cephei, slowly pulsating B-type stars (SPB), pulsating B stars exhibiting emission lines (Be), or rapidly oscillating Ap stars are also becoming of prime importance. These stars are more massive and hotter than the Sun and display pressure modes (e.g., δ Sct), gravity modes (e.g., γ Dor), or both (hybrid). So far, asteroseismology for γ Dor and δ Sct stars has been hindered due to difficulties in mode identifica-tion, rotational splitting, combination frequencies, mode selection, and mismatch between observed mode spectrum and theoretical predictions. However, very significant progress has been made in the last few years in deciphering the pulsation spectra, providing hope that asteroseismology may be feasible with these stars. Bedding et al. (2015) use échelle diagrams to study γ Dor g-mode period spacings, and identify sequences of rotationally split multiplets of degrees 1 and 2. Van Reeth et al. (2015a,b) show how to use nonuniform period spacings in Kepler γ Dor stars as a seismic diagnostic. García  report progress in developing a method to use low-order p modes to determine the mean density of δ Sct stars, and apply it successfully to a CoRoT star. Furthermore, Mirouh et al. (2019) report theoretical studies about period spacing in fast-rotating δ-Sct stars, and Kurtz et al. (2015a) show that complex frequency spectra of γ Dor, SPB, and Be stars can be explained by just a few g-mode frequencies plus their combination frequencies. It appears that the community is on the verge of a breakthrough in interpreting the frequency spectra of MS pulsators and developing asteroseismic techniques to derive interior structure properties.
Given the importance of asteroseismology for both solarlike stars and classical pulsators, it is of fundamental importance to identify a set of benchmarks with which to refine the stellar models. Such benchmarks are stars whose main physical properties are known with high precision, especially mass, radius, metallicity, and temperature. In recent decades, eclipsing binaries (EBs) have become very popular as benchmarks for stellar physics, as they provide accurate ways to measure masses, In the context of the preparation of the new ESA space mission PLATO, where asteroseismology plays a key role, we estimate that it is a good time for an inventory of the stellar pulsators in EBs, a unique class of stellar benchmarks. This motivated us to lead the first systematic search for any stellar pulsators in EBs in the Kepler data. We focus on the original Kepler mission because we can legitimately consider the list of EBs to be complete, whereas the catalog of the extended K2 mission is not complete. We consider the sample of EBs that is publicly available on the Villanova webpage, in its most recent update (May 2018), and the systems from the Coughlin et al. (2011) paper. The sum of the two catalogs contains 2925 systems, with orbital periods ranging from 0.05 to 1087 days. This global catalog is not in a strict sense a catalog of EBs, as about 600 ellipsoidal binaries are included. Ellipsoidal binaries are noneclipsing tight binaries where the stellar oblateness (ellipsoidal shape) introduces a periodic modulation of the light curve. We focus on γ Dor, δ Sct, and RG oscillators, but we also list all other types of pulsators that we are able to identify, in particular tidal pulsators, and a handful of possible SPBs or white dwarf (WD) pulsators.
Our goal is to decipher whether or not oscillations are detected in the Kepler EBs. We classify the oscillators according to their types, but we do not model each eclipse light curve or attempt to identify the oscillation modes that we detect. The objective is to motivate future in-depth studies of the most interesting cases, which would involve complementary radial velocity (RV) measurements. The paper is organized as follows. We review the main properties of the Villanova EBs and we describe the methods employed to disentangle the stellar pulsations from the eclipse signal in Sect. 2. We then present the detection of 303 systems where pulsations are detected in Sect. 3, including 163 that have been discovered in the present study, before discussing the results in Sect. 4.

Source of information
Most of the paper is based on the Villanova EB database (e.g., Prša et al. 2011;Slawson et al. 2011;Matijevič et al. 2012;Conroy et al. 2014;Kirk et al. 2016), which contains information about 2909 systems that were extracted with the "eclipsing binaries via artificial intelligence" (EBAI) pipeline (Prša et al. 2008). EBAI is a pipeline based on an artificial neural network that automatically extracts the sum of relative radii (R 1 + R 2 )/a, temperature ratio T 2 /T 1 , orbital argument of the periastron ω, eccentricity e and inclination angle i of a binary system. From the database, it is possible to download estimates of the orbital period, the date of primary eclipse, primary and secondary (if any) eclipse depths, durations (widths), relative separation, eclipse morphology, Kepler magnitude, and effective temperature from up to three different sources (Kepler Input Catalog, Pinsonneault et al. 2012, andCasagrande et al. 2011 catalogs). Regarding data availability and history, the database also indicates whether short cadence (1 min) or long cadence (30 min) data are available, and whether or not publications mentioning a given system exist. In addition, the detection of ETVs -an indicator of interacting triple systemsis flagged. We optimized our data processing tools for longcadence data and do not consider possible short-cadence data in this paper.
In Figs. 1 and 2, we present the distribution of the EB sample as a function of orbital period and effective temperature. The shortest orbital period is 0.05 days, the maximum is 1087 days. Ten percent of the stars display a period shorter than 0.33 days A106, page 2 of 17 (tenth percentile), and 10% longer than 43.3 days. The median orbital period is 2.30 days (Fig. 3). Out of the 2625 systems with an estimated effective temperature, the histogram reveals that 545 are in the [4500, 5300] K range, in which we expect RG oscillators -together with K dwarfs -, and 785 lie in the [6000-7500] K range where both δ Sct and γ Dor stars are expected. The eclipse light curves were classified by Matijevič et al. (2012) who introduce the "morphology" parameter c, which is a measure of "detachedness" of the system. This parameter is built upon a grid of simulated light curves that range from well-detached binaries to "overcontact" stars. According to these latter authors, all systems with c < 0.5 are predominantly detached. The range of c for semidetached systems broadly lies in the 0.5 < c < 0.7 range. Overcontact systems dominate the 0.7 < c < 0.8 region, after which a mixture of ellipsoidal variables and systems with uncertain classification sets in, including many HB systems. Figure 3 shows the histogram of the typology of the systems as a function of orbital period based on this rough classification. We note that among the 2909 systems, 175 are not classified. Out of the 2734 systems that are classified, 1431 (52.3%) are detached systems, 413 are semidetached (15.1%), 290 (10.6%) contact binaries, and 599 (21.9%) are classified as ellipsoidal or miscellaneous. It is very likely that most ellipsoidal binaries with periods shorter than one day are either contact or semidetached systems. It is also worth noting that detached systems dominate the sample starting from the median period (P ≥ 2.30 day).
Even though most of the paper is based on the Villanova database, a handful of systems that were published by Coughlin et al. (2011) have never been included in the Villanova database. We count 16 of these systems and include them in our analysis. We also add the RG in an eclipsing HB system recently studied by Kuszlewicz Fig. 3. Eclipse light curve classification according to Matijevič et al. (2012) as a function of orbital period. The size of each bloc is adjusted to include 10% of the total number of systems that were classified. The total number of classified systems is 2734 and each chunk consists of 273 ± 1 systems. The median orbital period is 2.30 days, while the 10 and 90 percentiles are 0.33 and 43.3 days, respectively. The labels "D", "SD", "OC", and "ELL" refer to detached, semidetached, overcontact and ellipsoidal systems. A total of 1431 systems are D, 413 are SD, 290 are OC, and 599 are ELL, i.e., a total of 2135 systems are classified as EBs.
in the Villanova catalog. We performed our study with Kepler public light curves available on the Mikulski Archive for Space Telescopes (MAST) 2 . We work with both the Simple Aperture Photometry (SAP, i.e., raw data) and the Pre-search Data Conditioning Simple Aperture Photometry (PDCSAP) light curves. The latter consist of time series that were corrected for discontinuities, systematic errors and excess flux due to aperture crowding (Twicken et al. 2010). We note that we do not seek complementary information from the Gaia data release 2 catalog because binary stars are not included, and if some are, their parameters may be biased.

Disentangling eclipses, surface activity, and oscillations
Searching for pulsations in EBs entails removing the photometric variability caused by binarity: mainly eclipses and phase effects. The methods we use are described in Gaulme et al. (2013Gaulme et al. ( , 2016 and here we provide a summary. The major challenge in concatenating light curves is to ensure photometric continuity before and after each interruption. The main cause of photometric jumps from quarter to quarter is the fact that the Kepler telescope rotated four times a year, which implied that a given star would fall on four different chips. However, the pointing was fine enough that a star repeatedly covered the same group of pixels every four quarters. Light curves are obtained by adding the pixels of the masks that were designed for every star of the field of view. For a given star, a mask was designed for each of Kepler's positions. Because of the photo response nonuniformity (PRNU) of the pixels and the changing size of the masks, the recorded flux changes. Both PRNU and varying mask areas lead to flux discontinuity that should be adjusted in a multiplicative way. The first correction we apply is therefore a normalization that turns the photoelectric counts into relative flux, by dividing the light curve of each quarter by its average. A median is actually more appropriate than a mean as outliers and large photometric jumps can bias the mean. If photometric variations would only be generated by PRNU and masks, this process should be enough. We note that this is only true for systems where no stellar activity is measurable if we exclude the effect of the differential velocity aberration.
Issues arise with the systems that display strong pseudoperiodic luminosity fluctuations. For those, the average (or median) over a quarter is biased by the fact that the number of pseudo periods is too small to be averaged out. Therefore, the median is not a perfect estimator of the mean photometry. This is an intrinsic limitation of the light curve photometric accuracy. In such cases, jumps remain with amplitudes within a few percent. Given that the remaining jumps are caused by a biased normalization, the second layer of adjustment to be applied should still be done in a multiplicative way. However, this is not possible in practice because none of the quarters can be considered as an absolute reference. The only corrections we may apply are additive, to ensure a smooth aspect of the light curve and to minimize their effects in the Fourier domain.
We employ two ways to smooth the remaining discontinuities once quarters are divided by their median. When a gap is short with respect to the photometric variability timescale, each side of the gap is adjusted accordingly. When a gap is longer than the variability period, we simply adjust the photometry with the difference of the means of each chunk surrounding the gap. Once the complete time series is leveled and concatenated, a linear fit is subtracted from it to take into account the decreasing instrumental sensitivity. Finally, when working with SAP light curves, we compensate for the differential velocity aberration -the motion of the target across a fixed aperture smaller than the point spread function -caused by the pixel scale breathing along the orbit of the satellite (372.5 days), whose peak-to-peak amplitude ranges from 0.5% to ∼10%. This is done by subtracting from each light curve a 372.5-day period sine fitting and a first harmonic, which is enough to reduce its amplitude to less than 0.5%.
We search for stellar pulsations in the power density spectra of the light curves. To minimize the effects of the incomplete duty cycle, we perform gap fillings and make use of the fast Fourier transform. All short gaps (only several missing points) are interpolated with a second-order polynomial estimated from the nearby data points. Long gaps are filled with zeros. To reduce the impact of abrupt discontinuities around long gaps, the edges of each section in between gaps are apodized with a cosine function. This is particularly important when significant variability is detected.
For well-detached systems, in which the time spent during eclipses is less than 20% of the time, we remove the data corresponding to the eclipses and bridge them with a secondorder polynomial constrained by the surrounding data. In cases where time spent during eclipses is more than 20%, we prefer to fold the light curve on the orbital period and subtract the average signal from each orbital period. Another option could have been to model the eclipse light curve with a fitting routine, such as PHOEBE (Prša 2018) or JKTEBOP . However, the advantage of detrending the light curves with a mean folded light curve is that the signal does not have to be modeled, which represents a huge amount of work for a sample of 3000 targets. In addition, residuals of a light curve minus a model are often significant enough to alter the signal. Still there is a big drawback when opting to subtract the folded light curve: because of pointing jitter and PRNU, the amplitude of eclipses may vary from quarter to quarter, or even during a quarter. Consequently, there are significant residuals in the light curves, and a signature of the eclipse signal is still quite visible. In such cases, we then remove all harmonics of the orbital period from the Fourier spectra. The latter method is not very sophisticated but can still help in detecting classical pulsator oscillations, as illustrated in Fig. 5 (middle panel, gray vs. black curve).

Data inspection tool
We developed a quick-look data-inspection tool (DIT) that helps the user to classify a system at a glance ( Fig. 4) 3 . For each star, it is composed of a page divided into a series of ten plots. First, the original Kepler light curves (SAP and PDCSAP) are displayed to look for any major issue regarding the data (panel a). This is useful in the case where a quarter or two are outliers with respect to the others and need to be manually removed before reprocessing the light curve. In a second panel b, the three types of light curves described in the previous section (variability, eclipse, oscillation types) are overplotted to make sure that the disentangling is correct. Specifically, the user checks whether discontinuities are still present in the oscillation-optimized light curve. A fold of the complete light curve is also a good indicator of any misestimate of the orbital period, or of the presence of a third body eclipsing or causing ETV (panel c). A zoom on each eclipse is also displayed for refining this analysis for long orbits, where the eclipses represent less than 10% of the period (panels d and e). We also plot a zoom on a relatively short range (20 days) at a random location, to help the user visually identify slow pulsations in the time domain, such as γ-Dor or tidally induced (panel f).
The last four plots are the power spectral density (PSD) in log-log scale to identify surface activity (typical of solar-like stars) and outstanding peaks indicating oscillations (panel g). It is also useful to check the quality of the background fitting that is performed to whiten the PSD. The background fitting is done following the prescription of Kallinger et al. (2014) for solar-like oscillators. The next plot is the square root of the previous one, that is, the Fourier transform module, which is more appropriate to identify pulsators with a large variety of oscillation amplitudes, such as δ-Scuti (panel h). The last two plots are more focused on the search for solar-like oscillators, as they display the envelope of the autocorrelation of the time series (EACF; see Mosser & Appourchaux 2009), which is usually considered the best performing tool to detect solar-like oscillations and measure the large frequency spacing ∆ν, even in low signal-to-noise-ratio (S/N) conditions (panel i). The last panel is an échelle diagram of the PSD based on the ∆ν deduced from the EACF computation (panel j).
In addition to the plots, the title includes the orbital period, the data sampling (long or short cadence), the light-curve type (SAP, PDCSAP), and the effective temperature. The latter is particularly useful when oscillations are detected. For example, an effective temperature of 4900 K is compatible with the detection of RG oscillations, and a temperature of 7000 K with a δ Scuti. We are aware that published temperatures may be inaccurate, or often not available, but it is valuable support information.

Results
We ran the DIT on all of the targets from the 2018 update of the Villanova catalog, except for the RGs in EB or HB systems published by Hekker et al. (2010), Gaulme et al. (2013Gaulme et al. ( , 2016, Beck et al. (2014Beck et al. ( , 2015, and the forthcoming paper by Benbakoura et al. (in prep.). Indeed, these systems were already well characterized, and first author Gaulme had already processed them. Regarding the classical pulsators, we were not sure about the exact number of systems already known and we did not exclude them a priori from the analysis. The total amount of stars for which we applied the DIT pipeline and that we visually inspected is 2875 out of a total of 2925. . Data inspection tool applied to the EB system KIC 3851949, which includes a solar-like oscillator with ν max ≈ 120 µHz and depleted = 1 modes. Left column, from top to bottom, panel a: Kepler light curves as a function of time (days), where "raw" stands for SAP and "cor" for PDCSAP. Panel b: light curves expressed in relative flux, where the blue line contains the stellar activity and oscillations (eclipses clipped out), the green line is the eclipse signal (activity filtered out), and the red curve is optimized for oscillation searches (activity and eclipses filtered out). Panel c: eclipse light curve folded on the orbital period. Panels d and e: zooms of the folded light curves around the eclipses. Right column, from top to bottom, panel f: zoom of the light curve over 20 days. Panel g: log-log scale display of the power spectral density of the time series expressed in ppm 2 µHz −1 as a function of frequency (µHz). Panel h: amplitude Fourier spectrum of the time series as a function of frequency. Panel i: envelope of the autocorrelation function (EACF) as a function of frequency and time. Panel j: échelle diagram associated with the large frequency spacing automatically determined from the EACF plot. The x-axis is the frequency modulo the large frequency spacing (i.e., from 0 to ∆ν), and the y-axis is frequency.  We first inspected every light curve to check whether orbital parameters were off. In 100 cases, we had to modify the ephemeris published on the Villanova database; many were small inaccuracies regarding eclipse timing and duration. Also, for 36 systems, only primary eclipses were detected. Thanks to the visual inspection of the light curves folded on the orbital period, we identified secondary eclipses for 36 new systems.
Even though it can sometimes be difficult to distinguish the detection of oscillations from imperfect light curve cleaning, we iteratively converged to the detection of pulsations in 303 out A106, page 5 of 17 of the 2925 light curves, that is, in about 10.4% of the cases. To do so, first author Gaulme produced a first screening by inspecting all of the DIT files, and selected about 350 of them. He attempted a preliminary classification according to pulsator type. Co-author Guzik then reviewed all of them and commented on each case with no instruction from Gaulme. Subsequently, Gaulme led a systematic search of every target on a web search engine by entering "KIC" and the KIC number of each star, one by one, with quotes and without quotes. This ended up being more efficient than using publication search tools such as NASA ADS or SIMBAD. The result of this search is certainly not fully exhaustive but is relatively complete. In total, we found 187 systems to have been studied already in one way or another in a peer-review paper or a conference proceeding. The level of study was highly variable, from simple notes in large tables containing many Kepler targets to dedicated papers about single binary systems. Having 187 systems previously cited in the literature does not mean that all were characterized both as binary and pulsators. Indeed, many were studied as double or triple systems but there was no mention of any stellar pulsation. Only 140 systems had already been reported as being part of a multiple-star system and containing a stellar pulsator. The main result of the present paper is that we have identified 163 new pulsating stars in multiple star systems. This represents 54% of the 303 stellar pulsators identified in the Kepler EB catalogs, and an increase of 116% of the known stellar pulsators in EBs among the Kepler data.
Not all of the pulsators we list are actual pulsators in multiple-star systems. Many are false positives, that is, pulsators that are either aligned (within the point spread function or pixel) with an EB or bright stellar pulsators whose light leaks into the considered EB mask. Disentangling those is a huge work that goes well beyond the scope of the present work. To do so, one should first check the Kepler imagettes and check whether the maximum intensity matches the maximum amplitude of the oscillations and the maximum depth of the eclipses. Such an approach was used by Gaulme et al. (2013) for 70 RG candidates in EBs. The second step would consist of consulting existing RV measurements or acquiring new ones. However, we can get a very rough idea of the likelihood that a pulsator belongs to a given binary system. For example, if we detect very clear oscillation peaks in the Fourier domain, while the photometric variability associated with the eclipses is significantly less than 1%, we can assume that the pulsator is not a member of the EB. Nevertheless, it could still be a triple system. Another example, a 10 R RG associated with a binary system that orbits in less than 5 days is very unlikely. Indeed, in such systems, stars are synchronized and the rotation rate of the star could make it disintegrate, depending on the orbital parameters. For RGs specifically, Gaulme et al. (2013) studied the question and used this type of consideration to detect many false positive RGs in EBs.
Our results are displayed in Table A.1, which includes the 303 systems where we are confident to have detected oscillations, plus the system identified by Kuszlewicz et al. (2019). The table is sorted by increasing KIC number; it is optimized to occupy as little room as possible and cannot include all the existing information. We therefore refer to the Villanova catalog and the Coughlin et al. (2011) paper for complete details of orbital parameters. We only list the main properties an observer would need to decide which target to study and observe: effective temperature, Kepler magnitude, orbital period, deepest eclipse depth, phase separation in between primary and secondary eclipses, ratio of eclipse durations, sum of eclipse durations relative to orbital period, and pulsator type. We indicate the number of eclipses per orbital period, which dictates whether it is a true EB or an ellipsoidal-variation binary, or even a HB star. We also list the references we have identified and systems marked with "Y" are those for which this is the first they have been identified as both a pulsator and a binary. Some notes were added to the last column to point out random relevant information.

Discussion
In this section, we highlight the main findings of the current study and list actions that could be taken based on this sample for future studies.

Classical pulsators
Classical pulsators represent the majority of the pulsators in EBs, as we expect from the effective temperature histogram (see Sect. 2.1). A total of 190 γ Dor, δ Sct, or hybrids are identified, including 122 that were previously unknown. There is more uncertainty regarding the detection of γ Dor oscillations with respect to δ Sct because of the frequency range that often overlaps harmonics of the orbital frequency. Due to the large amount of data, we did not look for regular period spacings or whatever type of pattern which would indicate a γ Dor pulsator. We simply flagged as γ Dor pulsators any stars that display pulsation frequencies of less than ∼50 µHz (∼5 c/d) and possibly typical γ Dor features in the time domain (e.g., Kurtz et al. 2015a). In addition, any information about effective temperature helped us to assess γ Dor pulsators. We flagged δ Sct stars whose oscillation spectra resembled those of typical δ Sct (e.g., Baglin et al. 1973;Balona et al. 2015), with frequencies larger than 50 µHz. However, in some cases, oscillation spectra are sparse (few peaks), and there could be some doubt surrounding tidally excited oscillations when peaks are regularly spaced. In some other cases, oscillation spectra closely resemble those of a δ Sct but the stellar effective temperatures are slightly too hot. In these cases we wrote a note in the table, but we keep in mind that effective temperatures from the KIC may be inaccurate, especially for close binary systems. We display three examples of γ Dor and δ Sct pulsators in Fig. 5.
A rather surprising finding of this study is the relatively large number of γ Dor and/or δ Sct pulsators in very-short-period systems. A total of 13 systems with orbital periods of less than 0.5 days display clear oscillations. Some could either be false positives or triple systems. Indeed, three of them are triple systems as ETVs were measured and published in previous articles (Gaulme et al. 2013;Borkovits et al. 2016). However, a more detailed investigation would be worthwhile to decipher whether the ten remaining systems are false positives, triples, or genuine EBs.
The first panel of Fig. 5 is an example of our short-period classical pulsators with the 0.32-day OC binary KIC 8330092. The oscillation spectrum is very clear, with a low-frequency part which we assume to be γ Dor and higher frequencies that are typical of δ Sct. In addition, the crowded aspect of the low-frequency part could be partially composed of Rossby modes as proposed by Saio et al. (2018). In the specific case of KIC 8330092, since the system is an OC, its rotation period is equal to its orbital period. From the light curve, the system seems to be composed of two similar stars (similar eclipse depth and shape), which we assume to be A0-to F5-type stars, that is, with radii and masses of about 1.4 R and 1.4 M . Within such an assumption (M 1 = M 2 , R 1 = R 2 ), the rotational velocity at the equator is V eq ≈ 220 km s −1 , whereas the escape velocity from the gravitational field is V esc ≈ 730 km s −1 , which means that such a star is far from disintegrating due to fast rotation. Moreover, still within the same assumptions, Kepler's third law indicates that the semi-major axis of the system would be a ≈ 2.8 R . These numbers are very rough estimates, but they indicate that a star such as a δ Sct can exist in such a tight contact binary.
Discovering classical pulsators in short-period binary systems is not something new per se. Indeed, several papers have reported detections of δ Sct or γ Dor pulsators in tight binary systems. Here, we mention a few examples. Aerts et al. (2002) report the detection of δ Sct pulsations in the 1.15-day orbit ellipsoidal binary XX Pyx. Dal & Sipahi (2013) report the detection of δ Sct oscillations in the 0.69-day orbit ellipsoidal binary V1464 Aql, and Sipahi & Dal (2014) the detection of γ Dor oscillations in the 0.93-day orbit close binary systems. Zhang et al. (2015) reported the detection of δ Sct oscillations in a component of the 0.65-day orbit near-contact binary V392 Ori. From the present study, the ten systems displaying γ Dor and/or δ Sct with P orb < 0.5 days are anomalous as to our knowledge no oscillations have been detected so far in such short-period systems. In other words, what is new here is the relative amount of short-period binaries among the oscillators (7% with P orb < 0.5 d) but it is also the first time that we identify pulsators in systems with periods of less than 0.65 days. Although some may be false positives, it is still very likely that we have identified the δ Sct in a binary system with the shortest orbit ever recorded.
Beyond short-period systems and thanks to the large amount of δ Sct that we have, we can test observations that have been made in previous studies. As presented first in Soydugan et al. (2006) then confirmed on a larger sample by Liakos & Niarchos (2017), there exists a correlation between the pulsation frequency and the orbital period. For systems with P orb less than a threshold located in between 13 and 40 days, the dominant δ Sct pulsation tends to increase with decreasing orbital period (see Fig. 4 in Liakos & Niarchos 2017). In Fig. 6, we test whether we find the same trend from the 149 possible δ Sct that we identify. For all systems flagged as possible δ Sct pulsators, we measured the frequency and height of the largest peak belonging to the δ Sct domain, that is, by excluding the γ Dor region for the hybrids. There is no obvious trend visible from our sample. In the most populated part of the diagram, we note a small trend however: for 0.6 < P orb < 1.4 d the most likely pulsation frequency is 20 c/d, for 1.4 < P orb < 2.1 d this is 16 c/d, and for 2.1 < P orb < 3.2 d it is 13.5 c/d. Several reasons may be responsible for not confirming the correlation between the pulsation frequency and the orbital period. Firstly, we have no way of knowing which systems are false-positive δ Sct binaries, while Liakos & Niarchos (2017) used systems that were well characterized. Secondly, we have not looked at the possible short-cadence data for part of our systems. Therefore, our maximum frequency available is the Kepler long-cadence Nyquist frequency ( 283 µHz). This means that a fraction of the frequencies that we present are likely to be aliases of higher frequencies. To emphasize this bias and help the comparison, we kept the y-axis boundaries to be the same as that in Liakos & Niarchos (2017) and we plotted the logarithm in base 10 of the pulsation and orbital periods. Thirdly, the Kepler sample lacks long-period systems, which prevents us from having a view on long-orbit systems (8 δ Sct with P orb > 30 days). Once the sample we list here is fully characterized, in particular with the help of complementary RV measurements, it will be possible to lead such a study.
In addition to the pulsation-dominant frequency, we tried to see whether binarity suppresses pulsation of A/F stars, as was   Fig. 6. Pulsation properties of the 148 δ Sct EB candidates identified in this work. Top: pulsation periods of the highest peaks of the oscillation spectra P puls as a function of orbital periods P orb . Bottom and left x-and y-axes are in the same units as in Liakos & Niarchos (2017). As in this latter paper, the vertical line indicates a 13-day orbital period, while the horizontal lines indicate pulsation frequencies of 3 and 80 cycles per day (c/d), respectively. Top and right x-and y-axes are the same in days for P orb and in frequency for pulsations. Bottom: representation of the height of the highest peak per δ Sct spectrum as a function of P orb . In both panels the color scale indicates the number of systems per square, to help the reader look for trends, such as possible clustering. The two vertical dotted lines indicate orbital periods from 0.6 to 3.2 days.
observed by  for RGs in close binary systems. We looked for a dependence of the pulsation amplitude as a function of orbital frequency. The bottom panel of Fig. 6 does not seem to show any correlation. More generally it would be interesting to see whether the fraction of binaries that do not pulsate is different from the fraction of stars in the instability strip that do not pulsate. Answering this question is tricky as the boundaries of the A/F pulsators are not particularly clear due to a lack of precision in the stellar parameters.

Solar-like oscillators
As indicated in the introduction, all solar-like oscillators are RGs. We detect 85 RG oscillation spectra among all the binaries listed in the Villanova catalog and the Coughlin et al. (2011) Fig. 7. Bona fide RGs/EBs discovered in the present work that are possibly SB2. Top panels: Kepler light curves folded over the orbital period (black markers), with a rebinned version of it in 30-min steps (red line). Bottom panels: PSD as a function of frequency of the Kepler light curvesafter eclipse removal -centered around the RG oscillation modes. Power spectral density was smoothed with a triangular weighted moving average (three bins wide at half maximum) to highlight the oscillation mode visibility. positives or triple systems. Even though only RV measurements and close-up imaging could allow us to verify their nature, we already have a general idea of the nature of these systems (false positive, binary, triple).
First of all, by following the conclusions of Gaulme et al. (2013), we are suspicious of RGs associated with a system with an orbital period of less than 10 days, especially if it is oscillating . As an example, so far, the shortest oscillating RG in a binary system is the EB KIC 8702921 with a 19.38-day orbit . Second, the shape of eclipses gives a strong indication of the nature of the system. Since the red giant branch (RGB) represents a short period compared to the overall star lifetime, both components are relatively unlikely to be at the same evolutionary stage. Therefore, as observed in practice, the large majority of RGs that belong to binary systems are on the RGB and are composed of a MS star and an evolved star (e.g., Gaulme et al. 2013;Beck et al. 2014). Therefore, one of the two stars is much smaller than the other, which means that one eclipse displays a flat bottom, while the other has a round shape caused by the limb darkening law of the RG star. Systems composed of a white dwarf and a MS star produce similar light curves to those composed of an M-dwarf and a subgiant; occasionally a hot Jupiter transiting a faint star can also produce this kind of light curve as well, provided the stellar light reflected by the planet is enough to cause a secondary transit. The system effective temperature, even if possibly biased, is complementary information that can help confirm that a system hosts an RG.
Beyond eclipsing systems, many "false positive" RG/EBs are actually HT systems, that is, systems where a close EB orbits an RG. In rare cases, the close EB also eclipses the RG (Derekas et al. 2011), but in general this is not the case (e.g., Gaulme et al. 2013). In this kind of system, the RG tends to cause ETVs. The detection of ETVs is a strong indication that the RG belongs to the triple system. Overall, when no ETVs are measured and the orbital period is too short to be an EB including an RG component, we classify a target as "likely FP or triple". We note that even if we do not detect ETVs, we cannot totally exclude a triple system, as ETVs may exist in wide hierarchical triple systems at a level that is not detectable (low ETV amplitude, slow variations).
Based on these criteria, we estimate that 15 out of the 27 systems flagged as displaying RG oscillations are bona fide pulsators in multiple-star systems. Specifically, ten are EBs, one is an HB with no eclipse, and four are HTs. Another system is a possible RG/EB (KIC 10858117), but the oscillation S/N is too poor for a robust conclusion on whether or not there are indeed RG oscillations in the Fourier spectrum of the time series.
As mentioned above, solar-like oscillators in EBs are unique targets for testing and calibrating asteroseismology. To this aim, individual masses and radii are needed, which implies that the EB systems must be double-lined spectroscopic binaries (SB2) to be considered as asteroseismic test benches. Given that most systems are composed of an MS star with an RG, the flux contrast in between both components is large; the companion usually accounts for a maximum of 10% of the total flux. In practice, for a system to be SB2, the companion star must be hotter than the RG, that is, a G-or F-type dwarf (e.g., Gaulme et al. 2016;Hełminiak et al. 2016Hełminiak et al. , 2017a. With the present work, we can estimate how many RG/EB/SB2 are present in total in the Kepler sample. So far, 11 RB/EB/SB2s have been identified and studied with complementary RV data Rawls et al. 2016;Gaulme et al. 2016;Hełminiak et al. 2015Hełminiak et al. , 2016Brogaard et al. 2018;Themeßl et al. 2018). In the forthcoming paper by Benbakoura et al. (in prep.), three more systems are identified. Furthermore, among those three, one shows no secondary eclipses due to a large eccentricity and a grazing primary, and only masses can be estimated. We are therefore at 14 RG/EB/SB2. In addition, one system from the Gaulme et al. (2016) paper (KIC 8054233) which has a 1058-day orbit was classified as SB1, but it could be an SB2 -higher-S/N observations are needed to make a robust conclusion. Hitherto, the total amount is thus 15 at best.
Among the new systems identified in this paper, only KICs 8308347 and 10920813, with orbital periods of 165 and 54 days, respectively, show evidence that the temperature of the companion is higher than that of the giant ("T 2 > T 1 " in Table A.1). Another case, KIC 8129189 (54-day period), displays deep partial eclipses and a clear oscillation signal with ν max ≈ 250 µHz, which indicates that it could be composed of a small RG (for which oscillations are detected) and a smaller RG or a subgiant star. This latter could also be an SB2. The case of KIC 10491544 displays the same kind of eclipses, namely partial; with T 2 ∼ T 1 , an RG of about 6.6 ± 0.3 R , and an orbit of 23 days could be a bona fide RG. However, their depths vary from quarter to quarter in between 2 and 5% and the S/N on oscillation measurements is very high, whereas it has been observed that modes are greatly depleted in such short-orbit systems ). This latter could be a false positive.
Overall, it appears that the number of RG/EB/SB2s from the A106, page 8 of 17  Fig. 8. Example of pulsators that either mimic tidal pulsators (KICs 8153568 and 6531496) or are actual tidal pulsators (KIC 11572363). In the cases of KIC 8153568 and 6531496 we do not plot a folded light curve but an extract of it over a little more than an orbital period (black) as the oscillations are not exactly in phase with the orbital period. The red curve is the original light curve minus the folded signal. In the case of KIC 8153568, the oscillations occur on the star that is hidden during the deeper eclipses, as the modulation is suppressed during them. Right panel: light curve and Fourier spectrum of KIC 11572363, a 19-day HB system with tidally excited oscillations.
Kepler data will be at best 19, which is a small statistical ensemble to test a tool as widely used as asteroseismology.
To complement these systems, triple systems displaying ETVs may be a good option. Borkovits et al. (2016) showed that it is possible to estimate the mass of the RG and the total mass of the tight binary component of a HT system using only ETV measurements. However, the precision of such estimates is rather poor (e.g., 1.5 ± 0.5 M for the very clear ETV KIC 7955301). To improve the precision, complementary RV measurements even in the case of SB1 systems should drastically reduce the error bars on masses -to the percent level that is required to calibrate asteroseismology. We note that no radii measurements arise from HT systems that are not triply eclipsing. The total number of HT systems for which ETVs have been detected is 13, including 4 from the present study. By considering both EBs/SB2s and HTs, the total number of systems that could be used to calibrate asteroseismic measurements of RG stars is about 30 (maximum 32) from the original Kepler mission.
Overall, to extend the sample of asteroseismic calibrators, we could use noneclipsing binaries that are both astrometric and visual binaries. Marcadon et al. (2018) studied such a system hosting a MS solar-like oscillator. In this specific case though, the S/N of the RV measurement was not high enough to estimate the mass of the oscillating star to better than 5%, which was not sufficient. Similarly to visual binaries, bright EBs could be resolved with interferometric measurements and be able to measure the proper motion of each companion, leading to the same result as for visual ones. So far, typical magnitude limitation is about eight in the visible. A few bright noneclipsing binaries may be present in the Kepler sample, and should be revealed by the ESA Gaia mission in the coming years (data release 3 or 4 planned in the early 2020s).
In addition, more data from other space missions can be considered, but none are comparable to those of Kepler. The extended Kepler mission K2 was limited to 90 days, most TESS fields of view are limited to 27 days, and CoRoT fields lasted 180 days at most. Given that  showed that RG/EBs with orbits shorter than 120 days tend to show oscillation suppression, these three other options (K2, TESS, CoRoT) are less promising than Kepler, even though it will be necessary to consider them too. The ESA PLATO mission, whose launch is planned in late 2026, will surely provide new interesting targets, given its large field of view and long exposures (2 years).

Other types of pulsators
Beyond γ Dor, δ Sct, and RGs we have identified 62 systems where other types of oscillations are possibly detected, most of which (59) are tidally excited. The Kepler mission has led to the discovery of a large number of tidal pulsators and has attracted much interest regarding these (e.g., Welsh et al. 2011;Dong et al. 2013;Maceroni et al. 2014;Smullen & Kobulnicky 2015;Hambleton et al. 2013Hambleton et al. , 2016Hambleton et al. , 2018Lee et al. 2016aLee et al. ,b, 2017Shporer et al. 2016;Dimitrov et al. 2017). Not surprisingly, among the 59 possible tidal pulsators that we list in Table A.1 about half (30) have already been identified as such. We note that many (35) of them are not only tidal pulsators but are also γ Dor or δ Sct, or both. For example, KIC 4544587 was studied in detail by Hambleton et al. (2013): these authors noticed that there were self-excited pressure and gravity modes (δ Sct types and γ Dor respectively), but also tidally excited modes and tidally influenced p-modes, which were identified because their frequencies were harmonics of the orbital period. This latter study perfectly illustrates how complex the identification of tidally excited oscillations can be: several types of oscillation modes can be present and it is hard to identify and verify the tidally excited nature of some of them. In addition, harmonics of the orbital frequency can be caused by insufficient removal of the eclipse signal. This difficulty of identifying tidal oscillations was already pointed out by Aerts et al. (2010), prior to the first Kepler results.
As an example, Fig. 8 shows three examples of pulsators whose oscillation spectra resemble tidal ones, but two are not. The left panel shows an extract of the light curve and Fourier spectrum of KIC 8153568, a detached EB orbiting in 3.607 days which appears to be a typical tidal pulsator where oscillations of constant amplitude are observed all along the orbit. We note that the oscillation amplitude is suppressed during the deeper eclipses, which means that only one of the two components of the binary system oscillates. However, when folding the light curve over the orbital period, the periodic modulation averages out. The only peak that dominates the Fourier spectrum has a frequency of ν peak = 140.511 µHz, which is not an integer multiple of the orbital frequency (factor 43.79). Tidal oscillations are expected at integer multiples of the orbital frequency because excited modes do not oscillate at their natural frequencies, but rather at tidal forcing frequencies, i.e., orbital. When a pulsation is observed at not an exact integer multiple of the orbital frequency, it is almost A106, page 9 of 17 A&A 630, A106 (2019) certainly not a tidally excited pulsation. Also, tidally excited pulsations are almost always found in eccentric systems, which is not the case here. The effective temperature of KIC 8153568 (6800 K) makes it compatible with a δ Sct, even though the oscillation spectrum does not resemble that of a δ Sct, with only one peak surrounded by a few aliases caused by the amplitude modulation of the pulsations during the primary eclipses. Figure 9 displays the time series (where the eclipse signal is removed) folded over 44 × 1/ν peak (i.e., 3.624 days) and shows how the pulsation of the star is constant and coherent.
There are a handful of case similar to KIC 8153568: these also are "false" tidal oscillators, but for different reasons. For example, KIC 6531496 -a detached system with a 14.32-day orbit (Fig. 9, middle) whose light curve resembles that of KIC 8153568 -is of a different nature after further investigation. The time series displays a photometric modulation along the orbit, which appears as a series of peaks in the power spectrum. The dominant peak corresponds to a period of 0.5129 days, which is almost a twentyeighth of the orbital period. When the light curve is folded over 28 × 0.5129 days (14.36 days), a very coherent periodic signal is observed, similar to that of KIC 8153568 (Fig. 9, middle panel). However, the 0.5129-day modulation is not a sine curve but a typical OC binary light curve, which means that KIC 6531496 has no tidal oscillations but is instead either a quadruple system or two binary systems blending into each other. One could then argue that KIC 8153568 is not a binary system including a stellar pulsator, but another quadruple or another pair of blended binaries. However, the pulsation period of KIC 8153568 is 0.08 days which seems a little short to be a binary, even though the Villanova catalog presents several systems with periods down to 0.05 days. Besides, the amplitude of the photometric modulation of KIC 8153568 disappears during one of the two eclipses, which proves that it is intrinsic to one of the two stars.
Another system where the origin of a photometric modulation is unclear is the detached binary KIC 6805146 with a 13.78day orbit, which shows a sine modulation at approximately (but not exactly) P orb /7 (Fig. 9, bottom panel). Again, this cannot be a tidally excited mode because it is not an integer multiple of the orbital period, so the most likely option is that it is a contami-nating ellipsoidal binary with a period of 1.909 days. The final example is KIC 4677321 (described in Table A.1 but not displayed in any figure), where we observe two prominent peaks at frequencies of (18 ± 1/4)ν orb , which suggest a rotational splitting where the rotation period would be ≈4mP orb where m is the azimuthal order of the excited mode. Assessment of the nature of these apparently but unlikely tidal pulsations requires further study involving high-resolution spectroscopy.
Despite the series of examples of misleading cases that resemble tidally excited pulsators, most of the systems we flagged as hosting tidally excited pulsations are HB stars. The right panel of Fig. 8 displays an HB system (KIC 11572363) for which tidally excited modes are reported for the first time. There is a distinction between the HB signal and tidally excited oscillations. The HB signal is the result of strong gravitational distortions and heating during periastron passage and does not last the whole orbital period. Tidally excited modes are oscillations driven by the tidal force onto the internal structure of the star. It has been observed that some HB stars display tidal oscillations, which are often at exact multiples of the orbital frequency (e.g., Fig. 8, right panel), and some others do not. The latter case is often seen with HB systems including an RG star, as observed by Beck et al. (2014). We note that as observed by Thompson et al. (2012), it may be sometimes hard to determine what fraction of the peaks is due to the HB shape and what fraction is due to stellar pulsation.
Finally, a handful of systems display different types of oscillations. The KICs 5217733 and 6806632 display oscillations that resemble those of γ Dor/δ Sct and δ Sct, respectively, but their effective temperatures of ≈9200 K according to the KIC are a little large for this kind of pulsator (Fig. 10). They could be either δ Sct stars with overestimated temperatures or SPBs with underestimated T eff , which we expect to be about 11 000 K. The KICs 6889235 and 8223868, also known as KOI 74 and 81, were identified as likely white dwarfs orbiting an A and B star, respectively . In both cases, oscillations are visible and are likely of tidal nature. The KIC 7749504 is an ellipsoidal binary orbiting in 0.57 days and displays clear peaks (last panel in Fig. 10) and T eff ≈ 11 000 K. This system was part of the search by  for Maia variables that are postulated between δ Sct and SPB, and the author concluded that KIC 7749504 was a rotational variable. We confirm the presence of a series of multiple peaks at more than twice and four times the orbital period, without being able to completely discard the hypothesis of nonrotational peaks. Finally, KIC 11179657 is the only sdB pulsator observed among the EB catalog and was classified as such by Pablo et al. (2012).

Conclusion
In summary, we led the first systematic search for stellar pulsators in the likely-to-be complete Kepler EB catalog, which includes 2925 systems among which about 600 are actually ellipsoidal binaries. To this aim, we developed a dedicated data inspection tool that automatically processes the Kepler light curves, provided elementary information about EB properties is fed in (orbital period, epochs, and duration of eclipses). We focused on three main classes of stellar oscillators: classical δ Sct and γ Dor pulsators, solar-like oscillators which happened to be all RGs, and tidal pulsators.
We first inspected the output figures produced by the DIT (e.g., Fig. 4) and made our own classification. Then, based on a manual search over the internet for possible earlier studies regarding these systems, we established a list of 303 pulsators associated with an EB, including 163 that were reported for the A106, page 10 of 17 first time as pulsator and binary candidates. A total of 149 stars are flagged as δ Sct (100 from this paper), 115 stars as γ Dor (69 new), 85 stars as RGs (27 new), and 59 as tidally excited oscillators (29 new). There is some overlap among these groups, as some display δ Sct, γ Dor, and tidal oscillations altogether. Many of these systems are likely to be false positives, that is, when an EB light curve is blended with a pulsator, and only deeper studies -often involving complementary RV measurementswould allow us to confirm the EB nature of these candidate pulsators. False positives are expected to be more common among the shallow-eclipse and short-period systems.
Of particular interest, some very short-orbit systems (≤0.5 days) display stellar oscillations (δ Sct and γ Dor). Some of them are likely to be false positives, but we know from the literature that δ Sct pulsators have been identified in almost contact binaries, but not at periods of under 0.6 days. Some ground-based support will be led to determine their nature. The opportunity to detect oscillations of contact binaries would provide unique insight into these types of systems, which are indeed very common. Regarding RG oscillators, one of the current major goals is to identify a large sample of oscillators in EBs to test and calibrate asteroseismic inference on stellar masses. We were able to estimate that less than 20 RG/EB/SB2s are present in the whole Kepler dataset, and that hierarchical triple systems may help to bring the sample up to about 30. We also identified new tidally eccentric binaries (HBs) with tidal oscillations, as well as detached binaries that display regular coherent pulsations that resemble tidal modes, but with periods that are not integer fractions of orbital periods.
In the context of imagining the future of asteroseismology, the present discovery of stellar pulsators in EB candidates constitutes a valuable sample that merits further detailed study, especially with the help of complementary observations. Priority should be given to the brightest targets for which a fine characterization is obtained more quickly. Priority should also be given to systems that are unique, such as the three likely RG/EB/SB2s, or some short-period classical pulsators. The current authors are involved in high-resolution spectroscopic observations with the 3.5m telescope of the Apache Point Observatory. However, given the large number of systems, such a survey will not be performed with only one telescope. We invite scientists that are interested to contact us in order to coordinate possible future observations. Notes. Systems are sorted by increasing KIC number (first column). Second and third columns are the effective temperatures, and Kepler magnitudes from the KIC. The "Type" column indicates the type of binary systems: D for detached, SD for semi-detached, OC for over-contact, and ELL for ellipsoidal. The "Ecl" column indicates the number of eclipses that are visible per orbital period (0 for ELL systems, 1 or 2 for the others). Column 6 displays the orbital period P orb expressed in days rounded at two digits. Columns 7-10 indicate the eclipse depth of the deepest eclipse (primary) expressed in percent, the orbital phase separation "sep" in between the secondary and the primary eclipse, the eclipse duration ratio W 2 /W 1 , and the sum of eclipse duration relative to orbital period (W 2 + W 1 )/P orb . The pulsator type "Puls type" is either δ-Sct for δ Scuti, γ Dor for γ Doradus, "Rossby", "RG" for RG solar-like oscillator, "tidal", SPB for slowly pulsating B-type stars, "WD" for white dwarf, or sdB for subdwarf B. The type "rot" stands for surface stellar "rotation" that may mimic stellar pulsations. The literature references are compacted (see below). "New PB" indicates a target identified both as a pulsator and a binary. Systems identified for the first time are marked with a "Y". Notes are random relevant specific information regarding each system. "HB" stands for heartbeat, FP for false positive, KOI for Kepler object of interest, P orb for orbital period, "Ecl." for eclipse, "harm." for harmonics, T 1 and T 2 for star 1 star 2 effective temperatures, T eff refers to the published KIC effective temperature, "SB1" and "SB2" for single or double-lined spectroscopic binary. The label "P orb = P orb /2" means that the actual orbital period is half of what Villanova database reported.