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

TL;DR

This study systematically searched the Kepler eclipsing binary catalog for stellar pulsators, discovering 303 systems with oscillations, including many new detections, to enhance understanding of stellar evolution.

Contribution

It is the first comprehensive search for pulsators within Kepler eclipsing binaries, introducing an automated tool and identifying numerous new pulsating systems.

Findings

Identified 303 pulsating eclipsing binary systems in Kepler data.

Discovered 163 new pulsators among the identified systems.

Classified stars into δ Scuti, γ Doradus, red giant, and tidally excited pulsators.

Abstract

Eclipsing binaries (EBs) are unique targets for measuring precise stellar properties and constrain stellar evolution models. In particular, it is possible to measure at the percent level masses and radii of both components of a double-lined spectroscopic EB. 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 low-mass star studies. The Kepler mission has led to the discovery of thousands of classical pulsators such as $\delta$ 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 eclipsing binary catalog. The focus is mainly aimed at discovering $\delta$ Scuti, $\gamma$ 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 eclipsing binary database and identified 303 systems whose light curves display oscillations, including 163 new discoveries. A total of 149 stars are flagged as $\delta$ Scuti (100 from this paper), 115 stars as $\gamma$ Doradus (69 new), 85 stars as red giants (27 new), 59 as tidally excited oscillators (29 new). There is some overlap among these groups, as some display several types of oscillations. Despite many of these systems are likely to be false positives, i. e., 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.

Full text

1 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 to the red giant phase, 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 of comparing the oscillation global properties of a given star to the Sun’s, 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 main-sequence (MS) to red-giant (RG) stars (Chaplin & Miglio 2013; Chiappini et al. 2015). It also will 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 $\displaystyle\gamma$ Doradus ($\displaystyle\gamma$ Dor), $\displaystyle\delta$ Scuti ($\displaystyle\delta$ Sct), $\displaystyle\beta$ Cephei, slowly pulsating B-type stars (SPB), pulsating B stars exhibiting emission lines (Be), or rapidly oscillating Ap stars are becoming of prime importance too. These stars are more massive and hotter than the Sun and display pressure modes (e.g., $\displaystyle\delta$ Sct), gravity modes (e.g., $\displaystyle\gamma$ Dor), or both (hybrid). So far, asteroseismology for $\displaystyle\gamma$ Dor and $\displaystyle\delta$ Sct stars has been hindered because of difficulty of mode identification, 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 $\displaystyle\gamma$ Dor $\displaystyle 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 $\displaystyle\gamma$ Dor stars as a seismic diagnostic. García Hernández et al. (2016) report progress in developing a method to use low-order $\displaystyle p$ modes to determine the mean density of $\displaystyle\delta$ Sct stars, and they apply it successfully to a CoRoT star, and Mirouh et al. (2019) report theoretical studies about period spacing in fast-rotating $\displaystyle\delta$-Sct stars. Kurtz et al. (2015b) show that complex frequency spectra of $\displaystyle\gamma$ Dor, SPB, and Be stars can be explained by just a few $\displaystyle 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 main-sequence pulsators and developing asteroseismic techniques to derive interior structure properties.

Given the importance of asteroseismology for both solar-like stars and classical pulsators, it is fundamental to identify a set of benchmarks to refine the stellar models on. Such benchmarks are stars whose main physical properties are known with high precision, especially mass, radius, metallicity and temperature. In the past decades, eclipsing binaries (EBs) have become very popular as benchmarks for stellar physics, as they provide accurate ways to measure masses, radii, and distances. It is possible to determine the mass and radius of each component of a double-line spectroscopic binary (SB2) from measurements of eclipse photometry and radial velocities. It is also possible to measure the mass of stars belonging to highly eccentric binary systems nicknamed “heartbeat” (HB) stars because their lightcurves recall electrocardiograms (e.g., Welsh et al. 2011), hierarchical triple systems (HTs, e.g., Borkovits et al. 2016), and visual binaries (e.g., Marcadon et al. 2018).

Recent space missions based on high-precision photometry drastically changed the number of known eclipsing binaries, as well as the sensitivity of measurements (e.g., Prša et al. 2011; Debosscher et al. 2011; Coughlin et al. 2011). The original Kepler mission has discovered 2,925 EBs111according to Villanova’s Kepler EB catalog, updated on May 6, 2018, which includes 2909 systems (http://keplerebs.villanova.edu/), and Coughlin et al. (2011) which lists 82 systems among which 16 are not in the Villanova catalog., including over 800 with periods from 10 to 1100 days, over 50 HB systems (e.g., Thompson et al. 2012; Beck et al. 2014; Gaulme et al. 2014; Shporer et al. 2016), and 222 triples displaying eclipse timing variations (ETVs) (Borkovits et al. 2016). Beyond Kepler, the CoRoT mission has discovered a few thousands of EB systems (Cilia Damiani, priv. comm.), and Kepler’s extended mission, K2, has found almost 700 EBs in the first six fields of view (according to Villanova’s database).

As regards classical pulsators, much research about binary systems including a pulsator has been led for the past thirty years (Szatmary 1990). The majority concerns $\displaystyle\delta$-Sct pulsators. We refer the reader to the comprehensive reviews available in Liakos & Niarchos (2017) and Kahraman Aliçavuş et al. (2017). In brief, there are 199 confirmed cases of binary systems containing at least a $\displaystyle\delta$-Sct pulsator (Liakos & Niarchos 2017), 87 of which being detached or semidetached eclipsing binaries, the other being visual, ellipsoidal variables and spectroscopic binaries. Among the catalog of 199 $\displaystyle\delta$-Sct in binaries, 29 are targets of the NASA original Kepler mission, among which 16 are eclipsing binaries. Debosscher et al. (2011), who were the first to look for oscillators in eclipsing binaries observed by Kepler, detected 14 classical pulsators in EBs, of which five are either $\displaystyle\gamma$ Dor or SPB. In addition, Gaulme & Guzik (2014) reported the identification of eight bona-fide classical pulsators among the Kepler EBs, one including a $\displaystyle\gamma$ Dor and seven $\displaystyle\delta$ Sct. Note that in addition to the eclipsing binaries, the measurement of fine frequency modulation of classical pulsators observed by Kepler led to the discovery of 341 new binary systems (Shibahashi & Kurtz 2012; Murphy et al. 2018). These systems are wide binaries – otherwise no fluctuation of mode frequencies would be measurable – and are likely not to be eclipsing.

As regards solar-like oscillators belonging to EBs, all are red giants (RGs) detected by the Kepler mission (Hekker et al. 2010; Gaulme et al. 2013, 2014; Beck et al. 2014, 2015; Kuszlewicz et al. 2019, Benbakoura et al., in prep.). So far, eleven wide SB2 EBs including an oscillating RG have been fully characterized with the help of radial-velocity ground-based support (Frandsen et al. 2013; Rawls et al. 2016; Gaulme et al. 2016; Brogaard et al. 2018; Themeßl et al. 2018) and three more are part of the upcoming paper of Benbakoura et al. (in prep.). Note that an equivalent number of RG displaying oscillations have been detected in HB systems (Gaulme et al. 2013, 2014; Beck et al. 2014, 2015; Kuszlewicz et al. 2019), but most do not show eclipses and are single-line spectroscopic binaries (SB1s).

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 eclipsing binaries, which is a unique class of stellar benchmarks. This motivated us to lead the first systematic search for any stellar pulsators in eclipsing binaries 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 eclipsing binaries as about 600 ellipsoidal binaries are included. Ellipsoidal binaries are non-eclipsing tight binaries where the stellar oblateness (ellipsoidal shape) introduces a periodic modulation of the light curve. We focus on $\displaystyle\gamma$ Dor, $\displaystyle\delta$ 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 SPB or white dwarf (WD) pulsators.

Our goal is to indicate whether 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. In terms of organization 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 Section 2. We then present the detection of 303 systems where pulsations are detected, including 163 that have been discovered in the present study (Sect 3) before discussing the results (Section 4).

2 Method

2.1 Source of information

Most of the paper is based on the Villanova eclipsing binary 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 $\displaystyle(R_{1}+R_{2})/a$, temperature ratio $\displaystyle T_{2}/T_{1}$, orbital argument of the periastron $\displaystyle\omega$, eccentricity $\displaystyle e$ and inclination angle $\displaystyle 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, and Casagrande et al. 2011 catalogs). Regarding data availability and history, it is also indicated whether short cadence (1 minute) or long cadence (30 minutes) data are available, and publications mentioning a given system exist. Besides, the detection of ETVs – an indicator of interacting triple systems – is flagged too. We optimized our data processing tools for long-cadence 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 day, the maximum is 1087 days. Ten percent of the stars display a period shorter than 0.33 day (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 $\displaystyle[4500,5300]$ K range, in which we expect red giant oscillators – together with K dwarfs –, and 785 lie in the [6000-7500] K range where both $\displaystyle\delta$ Sct and $\displaystyle\gamma$ Dor stars are expected.

The eclipse light curves were classified by Matijevič et al. (2012) who introduce the “morphology” parameter $\displaystyle c$, which is a measure of “detachedness” of the system. According to them, all systems with $\displaystyle c<0.5$ are predominantly detached. This parameter is built upon a grid of simulated light curves that range from well detached binaries to overcontact stars. The range of $\displaystyle c$ for semidetached systems broadly lies in the $\displaystyle 0.5<c<0.7$ range. Overcontact systems dominate the $\displaystyle 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. 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 semi-detached (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 period shorter than one day are either contact or semi-detached systems. It is worth noticing that detached systems dominate the sample starting from the median period ($\displaystyle P\geq 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 et al. (2019), which is not listed in the Villanova catalog. We performed our study with Kepler public light curves available on the Mikulski Archive for Space Telescopes (MAST)222http://archive.stsci.edu/kepler/. 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 corrected from discontinuities, systematic errors and excess flux due to aperture crowding (Twicken et al. 2010). Note that we do not seek for complementary information from the GAIA data release 2 catalog because binary stars are not included, and if some are, their parameters may be biased.

2.2 Disentangling eclipses, surface activity, and oscillations

Searching for pulsations in eclipsing binaries entails removing the photometric variability caused by binarity: eclipses and phase effects mainly. The methods we use are described in Gaulme et al. (2013, 2014, 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 non-uniformity (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 each quarter’s light curve 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. As a matter of fact, this is 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 pseudo-periodic 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 satellite’s orbit (372.5 days), whose peak-to-peak amplitude ranges from 0.5 % to $\displaystyle\sim 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.

In the case of 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 second-order polynomial, constrained by the surrounding data. In the case where time spent during eclipses is more than that, we prefer folding the light curve on the orbital period and subtract the average signal to 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 (Southworth 2013). However the advantage of detrending the light curves with a mean folded light curve avoids having to model the signal, 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. Then, 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 rather gross, but can still help in detecting classical pulsator oscillations, as illustrated in Fig. 5 (middle panel, gray vs. black curve).

2.3 Data inspection tool

We developed a quick look tool that helps the user classify a system at a glance (Fig. 4)333The ultimate goal of the DIT is to be public, but it is still under development and needs more time.. For each star, it is composed of 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). It is useful in the case where a quarter or two are outliers with respect to the others and deserve to be manually removed before reprocessing the light curve. In a second panel (b), the three types of light curves described in the previous subsection (variability, eclipse, oscillation types) are overplotted to make sure that the disentangling is correct. Especially, the user checks whether discontinuities are still present in the oscillation-optimized light curve. Then, a fold of the complete light curve is 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, 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 $\displaystyle\gamma$-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, i.e., the Fourier transform module, which is more appropriate to identify pulsators with a large variety of oscillation amplitude, such as $\displaystyle\delta$-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 $\displaystyle\Delta\nu$, even in low signal-to-noise (SNR) conditions (panel i). The last panel is an échelle diagram of the PSD, based on the $\displaystyle\Delta\nu$ deduced from the EACF computation (panel j).

Besides 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 red giant oscillations, and a temperature of 7000 K with a $\displaystyle\delta$ Scuti. We are aware that published temperatures may be off beyond the fact that they are often not available, but it is valuable support information.

3 Results

We ran the data inspection tool (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. (2013, 2016); Beck et al. (2014, 2015) and the forthcoming paper Benbakoura et al. (in prep). Indeed, these systems were already well characterized, and first author Gaulme already had processed them. As regards 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 2925 total.

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 of the 2925 light curves, i.e., 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. Then co-author Guzik reviewed all of them and commented on each case with no instruction from Gaulme. Then Gaulme led a systematic search of every target on a web search engine by entering, one by one, “KIC” and the KIC number of each star, with quotes and without quotes. Except for pictures of soccer players KICking balls, it ended up being more efficient than only through publication search tools such as NASA ADS or SIMBAD. The result of this search is certainly not fully exhaustive but is pretty much complete. In total, 187 systems were already studied 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. Actually, only 140 systems were already published as both multiple-star system and stellar pulsator. The main result of the present paper is thus to have identified 163 new pulsating stars in multiple star systems. It 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.

Now, not all of the pulsators we list are actual pulsators in multiple-star systems. Many are false positives, i.e., pulsators that are either aligned (within the point spread function or pixel) with an eclipsing binary 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 of intensity matches the maximum of 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 getting existing or acquiring new RV measurements. 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 eclipsing binary. Still, it can be a triple system. Another example, a 10-$\displaystyle R_{\odot}$ 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 RG specifically, Gaulme et al. (2013) studied the question and used this type of consideration to classify many RG apparently in EBs as false positives.

Our results are displayed in Table LABEL:tab:maintable, 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 and is optimized to not take too much room and cannot include all the existing information. We therefore refer to the Villanova catalog and the Coughlin et al. (2011) paper to get all details about 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 tells whether it is a true eclipsing binary or an ellipsoidal-variation binary, or even an HB star. We also list the references we have identified and we indicate with “Y” is this is the first time a system is identified as both a pulsator and a binary. Then we add some notes in the last column to point out random relevant information.

4 Discussion

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

4.1 Classical pulsators

Classical pulsators represent the majority of the pulsators in eclipsing binaries, as we expect from the effective temperature histogram (see Sect. 2.1). A total of 190 $\displaystyle\gamma$ Dor, $\displaystyle\delta$ Sct or hybrid are identified, including 122 that were previously unknown. There is more uncertainty regarding the detection of $\displaystyle\gamma$ Dor oscillations with respect to $\displaystyle\delta$ Sct because of the frequency range that often overlaps harmonics of the orbital frequency. Because of the large amount of data, we did not look for regular period spacings or whatever type of pattern that would for sure indicate that a pulsator is a $\displaystyle\gamma$ Dor. We just flagged as $\displaystyle\gamma$ Dor pulsators stars that display frequencies less than $\displaystyle\sim 50\,\mu$Hz ($\displaystyle\sim 5$ c/d) and possibly typical $\displaystyle\gamma$ Dor features in the time domain (e.g., Kurtz et al. 2015b). In addition, any information about effective temperature helped us assessing $\displaystyle\gamma$ Dor pulsators. We flagged $\displaystyle\delta$ Sct stars whose oscillation spectra looked like typical $\displaystyle\delta$ Sct (e.g., Baglin et al. 1973; Balona et al. 2015), with frequencies larger than 50 $\displaystyle\mu$Hz. However, in some cases, oscillation spectra are sparse (few peaks), and there could be some doubts with tidally excited oscillations when peaks are regularly spaced. In some other cases, oscillation spectra look much like a $\displaystyle\delta$ Sct but the stellar effective temperatures are a little hot. In these cases we write a note in the table, but we keep in mind that effective temperatures from the KIC may be off, and especially for close binary systems. We display three examples of $\displaystyle\gamma$ Dor and $\displaystyle\delta$ Sct pulsators in Fig. 5.

A rather surprising finding of this study is the relatively large number of $\displaystyle\gamma$ Dor and/or $\displaystyle\delta$ Sct pulsators in very short period systems. A total of 13 systems with orbital periods less than 0.5 days display clear oscillations. We naturally can assume them to all be either false positives or triple systems. Actually, three of them are indeed triple systems as ETVs were measured and published in previous articles (Gaulme et al. 2013; Borkovits et al. 2016). However, it would be worth checking more into details whether the ten remaining systems are false positives, triples or genuine eclipsing binaries.

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 $\displaystyle\gamma$ Dor and higher frequencies that are typical of $\displaystyle\delta$ Sct. Besides, 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, i.e., with radii and masses of about $\displaystyle 1.4\ R_{\odot}$ and $\displaystyle 1.4\ M_{\odot}$. Within such an assumption ($\displaystyle M_{1}=M_{2}$, $\displaystyle R_{1}=R_{2}$), the rotational velocity at equator is $\displaystyle V_{\mathrm{eq}}\approx 220$ km s*-1*, whereas the escape velocity from the gravitational field is $\displaystyle V_{\mathrm{esc}}\approx 730$ km s*-1*, which means that such a star is far from disintegrating because of fast rotation. Moreover, still within same assumptions, Kepler’s third law indicates that the semi-major axis of the system would be $\displaystyle a\approx 2.8R_{\odot}$. Of course, these numbers are very rough estimates, but they indicate that a star such as a $\displaystyle\delta$ Sct can exist in such a tight contact binary.

Discovering classical pulsators in short period binary systems is not something new per se: several papers report detections of $\displaystyle\delta$ Sct or $\displaystyle\gamma$ Dor pulsators in tight binary systems. We here mention a few examples. Aerts et al. (2002) report the detection of $\displaystyle\delta$ Sct pulsations in the 1.15-day orbit ellipsoidal binary XX Pyx. Dal & Sipahi (2013) report the detection of $\displaystyle\delta$ Sct oscillations in the 0.69-day orbit ellipsoidal binary V1464 Aql, and Sipahi & Dal (2014) the detection of $\displaystyle\gamma$ Dor oscillations in the 0.93-day orbit close binary systems. Zhang et al. (2015) reported the detection of $\displaystyle\delta$ Sct oscillations in a component of the 0.65-day orbit near-contact binary V392 Ori. From the present study, the ten systems displaying $\displaystyle\gamma$ Dor and/or $\displaystyle\delta$ Sct with $\displaystyle P_{\mathrm{orb}}<0.5$ day are odd as no oscillations have been detected so far, as far as we know, 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 $\displaystyle P_{\mathrm{orb}}<0.5$ d) but it is also the first time that we identify pulsators in systems with periods less than 0.65 days. Although some may be false positives, it is still very likely that we have identified the $\displaystyle\delta$ Sct in a binary system with the shortest orbit ever.

Beyond short period systems and thanks to the large amount of $\displaystyle\delta$ Sct that we have, we can test observations that have been done in previous works. 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 $\displaystyle P_{\mathrm{orb}}$ less than a threshold located in between 13 and 40 days, the dominant $\displaystyle\delta$ 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 $\displaystyle\delta$ Sct that we identify. For all systems flagged as possible $\displaystyle\delta$ Sct pulsator, we measured the frequency and height of the largest peak belonging to the $\displaystyle\delta$ Sct domain, i.e., by excluding the $\displaystyle\gamma$ Dor region for the hybrids. It arises that there is no obvious trend visible from our sample. In the most populated part of the diagram, we note a small trend though: for $\displaystyle 0.6<P_{\mathrm{orb}}<1.4$ d the most likely pulsation frequency is 20 c/d, for $\displaystyle 1.4<P_{\mathrm{orb}}<2.1$ d is is 16 c/d, and for $\displaystyle 2.1<P_{\mathrm{orb}}<3.2$ d it is 13.5 c/d.

Several reasons may be responsible for not confirming that observation. Firstly, we have no way to know which systems are false-positive $\displaystyle\delta$ 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 ($\displaystyle\simeq 283\ \mu$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 $\displaystyle 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 of long-period systems, which prevents us to have a view on long-orbit systems (8 $\displaystyle\delta$ Sct with $\displaystyle P_{\mathrm{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 kind of study.

In addition to the pulsation dominant frequency, we tried to see whether binarity suppresses pulsation of A/F stars, as was observed by Gaulme et al. (2014) 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 than 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 is not very clear because the stellar parameters are not precise enough.

4.2 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) table. Among those, 27 are new to be both identified as oscillating RG and a binary (eclipsing, ellipsoidal or HB). As mentioned earlier, a large fraction of them must be either false positives or triple systems. Even though only RV measurements and close up imaging could allow us to make sure of their nature, we can already get some ideas on the nature of these systems (false positive, binary, triple).

First of all, by following the conclusions of Gaulme et al. (2013), we are suspicious with RGs associated with a system of orbital period less than 10 days, especially if is oscillating (Gaulme et al. 2014). As an example, so far, the shortest oscillating RG in a binary system is the 19.38-day orbit EB KIC 8702921 (Gaulme et al. 2016). Secondly, the shape of eclipses is a strong insight on the nature of the system. Since the RG branch (RGB) represents a short period compared to the overall star lifetime, it is rather unlikely to have both components on the same evolutionary state. Therefore, as observed in practice, the large majority of RGs that belong to binary systems are on the RGB and are composed of a main-sequence 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. Similar light curves can be produced with systems composed of a white dwarf and a main sequence star, with an M-dwarf and a subgiant, or occasionally by a hot Jupiter transiting a faint star provided the stellar light reflected by the planet is enough to cause a secondary transit. The system effective temperature, even if possibly biased, is a complementary information that can help confirm that a system hosts an RG.

Beyond eclipsing systems, many “false positive” RG/EBs are actually HT systems, i.e., where a close eclipsing binary orbits an RG. In rare cases, the close EB also eclipses the RG (Derekas et al. 2011), but in general it is not the case (e.g., Gaulme et al. 2013). In such kind of systems, 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”. Note that even if we do not detect ETVs, we cannot totally exclude a system to be triple, as ETVs may exist in wide hierarchical triple systems at a level that is not detectable (low ETVs 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, 10 are EBs, one is an HB with no eclipse, and 4 are HTs. Another system is a possible RG/EB (KIC 10858117), but the oscillation SNR is so poor that it is hard to be fully convinced that there are indeed RG oscillations in the Fourier spectrum of the time series.

As reminded earlier, solar-like oscillators in EBs are unique targets for testing and calibrating asteroseismology. For that goal, 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 account 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, i.e. a G or F type dwarf (e.g., Gaulme et al. 2016; Hełminiak et al. 2016, 2017b). 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 (Frandsen et al. 2013; Rawls et al. 2016; Gaulme et al. 2016; Hełminiak et al. 2015, 2016; Brogaard et al. 2018; Themeßl et al. 2018). In the forthcoming paper by Benbakoura et al., three more systems are identified. Actually 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 thus 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 turn into an SB2 with higher SNR observations in the future. 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 companion’s temperature is larger than the giant’s (“$\displaystyle T_{2}>T_{1}$” in Table LABEL:tab:maintable). Another case, KIC 8129189 (54-day period), displays deep partial eclipses and a clear oscillation signal with $\displaystyle\nu_{\mathrm{max}}\approx 250\ \mu$Hz, which indicates it could be composed of a small RG (for which oscillations are detected) and a smaller RG or a subgiant star. This latter could be an SB2 too. The case of KIC 10491544 displays the same kind of eclipses – partial –, with $\displaystyle T_{2}\sim T_{1}$, an RG of about $\displaystyle 6.6\pm 0.3R_{\odot}$, 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 oscillations SNR is very high, whereas it has been observed that modes are much depleted in such short-orbit systems (Gaulme et al. 2014). This latter could be a false positive. Overall, it appears that the number of RG/EB/SB2s from the 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, just based on ETV measurements. However, the precision of these estimate is rather poor (e.g., $\displaystyle 1.5\pm 0.5M_{\odot}$ 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. Note that no radii measurements arise from HT systems that are not triply eclipsing. The total amount of HT systems for which ETVs have sure be detected is 13, including 4 from the present study. By considering both EB/SB2s and HTs, the total number of systems that could help calibrating asteroseismic measurements of red giant stars is about 30 (maximum 32) from the original Kepler mission.

Overall, to extend the sample of asteroseismic calibrators, we could use non eclipsing binaries that are both astrometric and visual binaries. Marcadon et al. (2018) studied such a system hosting a main-sequence solar-like oscillator. In this specific case though, the SNR of the RV measurement was not high enough to estimate the oscillating star’s mass to better than 5 %, which was not sufficient. Similarly to visual binaries, bright EB 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 non-eclipsing 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).

Besides, more data from other space missions can be considered too, but none are comparable to Kepler. K2 is limited to 90 days, most TESS fields of view are limited to 27 days, and CoRoT fields were lasting 180 days at maximum. Given that Gaulme et al. (2014) 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 provide for sure new interesting targets, given its large field of view and long exposures (2 years).

4.3 Other types of pulsators

Beyond $\displaystyle\gamma$ Dor, $\displaystyle\delta$ Sct, and RG we have identified 62 systems where other types of oscillations are possibly detected, most of which (59) being tidally excited. The Kepler mission has led to discover a large number of tidal pulsators and attracted much interest regarding it (e.g., Welsh et al. 2011; Dong et al. 2013; Maceroni et al. 2014; Smullen & Kobulnicky 2015; Hambleton et al. 2013, 2016, 2018; Kjurkchieva et al. 2016; Lee et al. 2016a, b; Shporer et al. 2016; Lee et al. 2017; Dimitrov et al. 2017). Not surprisingly, among the 59 possibly tidal pulsators that we list in Table LABEL:tab:maintable about half (30) were already identified as such. Note that many (35) of them are not only tidal pulsators, but also $\displaystyle\gamma$ Dor or $\displaystyle\delta$ Sct, or both. For example, KIC 4544587 was studied in detail by Hambleton et al. (2013): they noticed that there were self-excited pressure and gravity modes ($\displaystyle\delta$ Sct types and $\displaystyle\gamma$ Dor respectively), but also tidally excited modes and tidally influenced p-modes, which were identified because their frequencies were harmonics of the orbital period. Their work perfectly illustrates how complex the identification of tidally excited oscillations might be: several types of oscillation modes can be present and it is hard to identify and make sure of the tidally-excited nature of some of them. Besides, harmonics of the orbital frequency can also 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 look like tidal ones, but where 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 day, which looks 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 $\displaystyle\nu_{\mathrm{peak}}=140.511\ \mu$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 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 $\displaystyle\delta$ Sct, even though the oscillation spectrum does not look like a $\displaystyle\delta$ Sct one, 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 $\displaystyle 44\times 1/\nu_{\mathrm{peak}}$ (i.e., 3.624 days) and it shows how constant and coherent is the pulsation of that star.

A handful of cases look similar to KIC 8153568: they also are “false” tidal oscillators, but for different reasons. For example, KIC 6531496 – a 14.32-day orbit detached system (Fig. 9, middle) – whose light curve looks like 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 close to being a 28th of the orbital period. When we fold the light curve over $\displaystyle 28\times 0.5129$ days (14.36 days), we observe a very coherent periodic signal as with 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, KIC 8153568’s pulsation period is 0.08 days which seems a little short to be a binary, even though the Villanova catalog presents some 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 13.78-day orbit detached binary KIC 6805146, which shows a sine modulation at about but not exactly $\displaystyle P_{\mathrm{orb}}/7$ (Fig. 9, bottom panel). Again, it 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 contaminating ellipsoidal binary with period of 1.909 days. A last example is KIC 4677321 (described in Table LABEL:tab:maintable but not displayed in any figure), where we observe two prominent peaks at frequencies $\displaystyle(18\pm 1/4)\nu_{\mathrm{orb}}$, which suggest a rotational splitting where rotation period would be $\displaystyle\approx 4mP_{\mathrm{orb}}$ where $\displaystyle m$ is the azimuthal order of the excited mode. Assessing the nature of these tidally-looking but tidally-unlikely pulsations requires some more studies, 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 in between the HB signal and tidally excited oscillations. The HB signal is the result of strong gravitational distortions and heating during periastron passage and it 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). 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. KICs 5217733 and 6806632 display oscillations that look like those of $\displaystyle\gamma$ Dor/$\displaystyle\delta$ Sct and $\displaystyle\delta$ Sct respectively but their effective temperatures of $\displaystyle\approx 9200$ K according to the KIC are a little large for such kind of pulsators (Fig. 10). It could be either a $\displaystyle\delta$ Sct with an overestimation of its $\displaystyle T_{\mathrm{eff}}$ or an SPB with an underestimated $\displaystyle T_{\mathrm{eff}}$, which we expect to be about 11,000 K. KICs 6889235 and 8223868, also known as KOI 74 and 81 were identified as likely white dwarfs orbiting an A and B star respectively (Rowe et al. 2010). In both cases, oscillations are visible and are likely of tidal nature. KIC 7749504 is an ellipsoidal binary orbiting in 0.57 days: it displays clear peaks (last panel in Fig. 10) and $\displaystyle T_{\mathrm{eff}}\approx 11000$ K. This system was part of Balona (2015)’s search for Maia variables that are postulated between $\displaystyle\delta$ Sct and SPB, and the author concluded it 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 non-rotational peaks. Finally, KIC 11179657 is the only sdB pulsator observed among the EB catalog and was classified as such by Pablo et al. (2012).

5 Conclusion

To summarize, we led the first systematic search for stellar pulsators in the likely-to-be complete Kepler eclipsing binary catalog, which includes 2925 systems among which about 600 are actually ellipsoidal binaries. For this, we developed a dedicated data inspection tool that automatically processes the Kepler light curves, provided elementary information about EB properties are fed in (orbital period, epochs and duration of eclipses). We focused on three main classes of stellar oscillators: classical $\displaystyle\delta$ Sct and $\displaystyle\gamma$ Dor pulsators, solar-like oscillators which happened to be all red giants, 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 first time as pulsator and binary candidates. A total of 149 stars are flagged as $\displaystyle\delta$ Sct (100 from this paper), 115 stars as $\displaystyle\gamma$ Dor (69 new), 85 stars as RGs (27 new), 59 as tidally excited oscillators (29 new). There is some overlap among these groups, as some display $\displaystyle\delta$ Sct, $\displaystyle\gamma$ Dor and tidal oscillations altogether. Many of these systems are likely to be false positives, i. e., when an EB light curve is blended with a pulsator, and only deeper studies – often involving complementary RV measurements – would allow us to make these candidate confirmed pulsators in EBs. False positives are expected to be relatively more common among the shallow-eclipse and short-period systems.

Among interesting facts are some very short orbit systems ($\displaystyle\leq 0.5$ days) displaying stellar oscillations ($\displaystyle\delta$ Sct and $\displaystyle\gamma$ Dor). Some of them are likely to be false positives, but we know from the literature that $\displaystyle\delta$ Sct pulsators have been identified in almost contact binaries, but not at period under 0.6 days. Some ground based support will be led to determine their nature. The opportunity of detecting oscillations of contact binaries would provide unique views on these types of systems that are very common. As regards red giant oscillators, one of the major current goal 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/SB2 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 (HB) with tidal oscillations, as well as detached binaries that display regular coherent pulsations that look like 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 eclipsing binary candidates constitutes a valuable sample that deserves to be studied in further detail, especially with the help of complementary observations. Priority should be given to the brightest targets for which a fine characterization is obtained more quickly, especially about the rate of false positives. Priority should also be given to systems that are unique, as the three likely RG/EB/SB2s, or some short period classical pulsators. The authors are involved in high-resolution spectroscopic observations with the 3.5-m telescope of the Apache Point Observatory. However, given the large number of systems, such a survey will not be performed only with one telescope. We invite scientists that are interested to contact us in order to coordinate possible future observations.

Acknowledgments

P. Gaulme was supported in part by the German space agency (Deutsches Zentrum für Luft- und Raumfahrt) under PLATO data grant 50OO1501. P. Gaulme acknowledges NASA grant NNX17AF74G for partial support. The authors thank Jason Jackiewicz at New Mexico State University for comments and text editing. P. Gaulme thanks Cilia Damiani at Max Planck Institut für Sonnensystemforschung, Jim Fuller at Caltech, Donald W. Kurtz at University of Central Lancashire, Dominic Bowman at KU Leuven, Kosmas Gazeas at University of Athens, and Giovanni Mirouh at University of Surrey for useful comments on the manuscript. The authors thank the Apache Point Observatory 3.5-meter telescope, which is owned and operated by the Astrophysical Research Consortium, for granting us observing time during quarters 2 and 3, 2019 to support the follow up of the present work.

Appendix A Table of pulsators in eclipsing binaries

Literature abbreviation used in Table LABEL:tab:maintable:

AS14 Kahraman Aliçavuş & Soydugan (2014), Ali17 Aliçavuş & Soydugan (2017), Arm14, Armstrong et al. (2014), Blj15 Balaji et al. (2015), Bal11 Balona et al. (2011), Bal13 Balona et al. (2013), Bal14 Balona (2014), Bal15 Balona (2015), Bal15b Balona et al. (2015), Bck14 Beck et al. (2014), Bed11 Bedding et al. (2011), Bel19 Bell et al. (2019), Bog15 Bognár et al. (2015), Bow16 Bowman et al. (2016), Bie17 Biersteker & Schlichting (2017), Bor11 Borucki et al. (2011), Brk14 Borkovits et al. (2014), Brk16 Borkovits et al. (2016), Bra15 Bradley et al. (2015), Bro18 Brogaard et al. (2018), Cei17 Ceillier et al. (2017), Con14 Conroy et al. (2014), Cou11 Coughlin et al. (2011), Deb11 Debosscher et al. (2011), Deb13 Debosscher et al. (2013), Dem15 Demircan & Bulut (2015), Dim17 Dimitrov et al. (2017), DR12 Dodson-Robinson (2012), Don13 Dong et al. (2013), Fra13 Frandsen et al. (2013), Gau13 Gaulme et al. (2013), Gau14 Gaulme et al. (2014), Gau16 Gaulme et al. (2016), GG14 Gaulme & Guzik (2014), Gie15 Gies et al. (2015), Guo16 Guo et al. (2016), Guo17 Guo et al. (2017b), Guo17b Guo et al. (2017a), Fai15 Faigler et al. (2015), Geo18 George et al. (2018), Hek10 Hekker et al. (2010), Ham13 Hambleton et al. (2013), Ham16 Hambleton et al. (2016), Ham18 Hambleton et al. (2018), Hel16 Hełminiak et al. (2016), Hel17 Hełminiak et al. (2017a), Hel17b Hełminiak et al. (2017b), KN18 Katsova & Nizamov (2018), Kju16 Kjurkchieva et al. (2016), Kju16b Kjurkchieva & Atanasova (2016), Kou18 Kouzuma (2018), Kur15 Kurtz et al. (2015a), Kus19 Kuszlewicz et al. (2019) Lee14 Lee et al. (2014), Lee15 Lee et al. (2015), Lee16a Lee et al. (2016a), Lee16b Lee et al. (2016b), Lee17 Lee et al. (2017), Leh13 Lehmann et al. (2013), LiB14 Lillo-Box et al. (2014), Lur17 Lurie et al. (2017), Lia17 Liakos (2017), Lia18 Liakos (2018), LN17 Liakos & Niarchos (2017), Mac14 Maceroni et al. (2014), MH18 Manuel & Hambleton (2018), Mas17 Masuda (2017), Mat17 Matson et al. (2017), McN12 McNamara et al. (2012), Moy17 Moya et al. (2017), Mur15 Murphy et al. (2015), Mur18 Murphy et al. (2018), Nie15 Niemczura et al. (2015), Nie17 Niemczura et al. (2017), Ozd17 Özdarcan & Dal (2017), Oro15 Orosz (2015), Pab12 Pablo et al. (2012), Ram14 Ramsay et al. (2014), Rap13 Rappaport et al. (2013), Rap15 Rappaport et al. (2015), Raw16 Rawls et al. (2016), Row10 Rowe et al. (2010), Row15 Rowe et al. (2015), Sow17 Sowicka et al. (2017), Sai18 Saio et al. (2018), Snd16 Sandquist et al. (2016), San16 Santerne et al. (2016), Shp16 Shporer et al. (2016), Smu15 Smullen & Kobulnicky (2015), Sou11 Southworth et al. (2011), Sou15 Southworth (2015), The18 Themeßl et al. (2018), Tho12 Thompson et al. (2012), TM16 Turner & Maynard (2016), Tur15 Turner & Holaday (2015), Tur19 Turner (2019), Uyt11 Uytterhoeven et al. (2011), Wel11 Welsh et al. (2011), Yak15 Yakut (2015), Zha17 Zhang et al. (2017a), Zha17b Zhang et al. (2017b), ZXB17 Zhang et al. (2017c), Zho10 Zhou (2010), Zim17 Zimmerman et al. (2017)