Finding binaries among Kepler pulsating stars from phase modulation of their pulsations

We present a method for finding binaries among pulsating stars that were observed by the Kepler Mission. We use entire 4 yr light curves to accurately measure the frequencies of the strongest pulsation modes, and then track the pulsation phases at those frequencies in 10-d segments. This produces a series of time-delay measurements in which binarity is apparent as a periodic modulation whose amplitude gives the projected light travel time across the orbit. Fourier analysis of this time-delay curve provides the parameters of the orbit, including the period, eccentricity, angle of ascending node, and time of periastron passage. Differentiating the time-delay curve yields the full radial-velocity curve directly from the Kepler photometry, without the need for spectroscopy. We show examples with delta scuti stars having large numbers of pulsation modes, including one system in which both components of the binary are pulsating. The method is straightforward to automate, thus radial velocity curves can be derived for hundreds of non-eclipsing binary stars from Kepler photometry alone.