The primary Kepler Mission provided nearly continuous monitoring of ∼200,000 objects with unprecedented photometric precision. We present the final catalog of eclipsing binary systems within the 105 deg 2 Kepler field of view. This release incorporates the full extent of the data from the primary mission (Q0-Q17 Data Release). As a result, new systems have been added, additional false positives have been removed, ephemerides and principal parameters have been recomputed, classifications have been revised to rely on analytical models, and eclipse timing variations have been computed for each system. We identify several classes of systems including those that exhibit tertiary eclipse events, systems that show clear evidence of additional bodies, heartbeat systems, systems with changing eclipse depths, and systems exhibiting only one eclipse event over the duration of the mission. We have updated the period and galactic latitude distribution diagrams and included a catalog completeness evaluation. The total number of identified eclipsing and ellipsoidal binary systems in the Kepler field of view has increased to 2878, 1.3% of all observed Kepler targets. An online version of this catalog with downloadable content and visualization tools is maintained at http://keplerEBs.villanova.edu.

We have discovered a class of eccentric binary systems within the Kepler data archive that have dynamic tidal distortions and tidally-induced pulsations. Each has a uniquely shaped light curve that is characterized by periodic brightening or variability at time scales of 4-20 days, frequently accompanied by shorter period oscillations. We can explain the dominant features of the entire class with orbitally-varying tidal forces that occur in close, eccentric binary systems. The large variety of light curve shapes arises from viewing systems at different angles. This hypothesis is supported by spectroscopic radial velocity measurements for five systems, each showing evidence of being in an eccentric binary system. Prior to the discovery of these 17 new systems, only four stars, where KOI-54 is the best example, were known to have evidence of these dynamic tides and tidally-induced oscillations. We perform preliminary fits to the light curves and radial velocity data, present the overall properties of this class and discuss the work required to accurately model these systems.

The precision of photometric and spectroscopic observations has been systematically improved in the last decade, mostly thanks to space-borne photometric missions and ground-based spectrographs dedicated to finding exoplanets. The field of eclipsing binary stars strongly benefited from this development. Eclipsing binaries serve as critical tools for determining fundamental stellar properties (masses, radii, temperatures and luminosities), yet the models are not capable of reproducing observed data well either because of the missing physics or because of insufficient precision. This led to a predicament where radiative and dynamical effects, insofar buried in noise, started showing up routinely in the data, but were not accounted for in the models. PHOEBE (PHysics Of Eclipsing BinariEs;http://phoebe-project.org) is an open source modeling code for computing theoretical light and radial velocity curves that addresses both problems by incorporating missing physics and by increasing the computational fidelity. In particular, we discuss triangulation as a superior surface discretization algorithm, meshing of rotating single stars, light time travel effect, advanced phase computation, volume conservation in eccentric orbits, and improved computation of local intensity across the stellar surfaces that includes photon-weighted mode, enhanced limb darkening treatment, better reflection treatment and Doppler boosting. Here we present the concepts on which PHOEBE is built on and proofs of concept that demonstrate the increased model fidelity.

We present Kepler photometry and ground based spectroscopy of KIC 4544587, a short-period eccentric eclipsing binary system with self-excited pressure and gravity modes, tidally excited modes, tidally influenced p modes, and rapid apsidal motion of 182 y per cycle. The primary and secondary components of KIC 4544587 reside within the δ Scuti and γ Dor instability region of the Hurtzsprung-Russell diagram, respectively. By applying the binary modelling software phoebe to prewhitened Kepler photometric data and radial velocity data obtained using the William Herschel Telescope and 4-m Mayall telescope at KPNO, the fundamental parameters of this important system have been determined, including the stellar masses, 1.98 ±0.07 M ⊙ and 1.60 ± 0.06 M ⊙ , and radii, 1.76 ± 0.03 R ⊙ and 1.42 ± 0.02 R ⊙ , for the primary and secondary components, respectively. Frequency analysis of the residual data revealed 31 modes, 14 in the gravity mode region and 17 in the pressure mode region. Of the 14 gravity modes 8 are orbital harmonics: a signature of tidal resonance. While the measured amplitude of these modes may be partially attributed to residual signal from binary model subtraction, we demonstrate through consideration of the folded light curve that these frequencies do infact correspond to tidally excited pulsations. Furthermore, we present an echelle diagram of the pressure mode frequency region (modulo the orbital frequency) and demonstrate that the tides are also influencing the p modes. A first look at asteroseismology hints that the secondary component is responsible for the p modes, which is contrary to our expectation that the hotter star should pulsate in higher radial overtone, higher frequency p modes.

Context. The unparalleled photometric data obtained by NASA's Kepler Space Telescope has led to improved understanding of red giant stars and binary stars. Seismology allows us to constrain the properties of red giants. In addition to eclipsing binaries, eccentric non-eclipsing binaries that exhibit ellipsoidal modulations have been detected with Kepler. Aims. We aim to study the properties of eccentric binary systems containing a red giant star and to derive the parameters of the primary giant component. Methods. We applied asteroseismic techniques to determine the masses and radii of the primary component of each system. For a selected target, light and radial velocity curve modelling techniques were applied to extract the parameters of the system and its primary component. Stellar evolution and its effects on the evolution of the binary system were studied from theoretical models. Results. The paper presents the asteroseismic analysis of 18 pulsating red giants in eccentric binary systems, for which masses and radii were constrained. The orbital periods of these systems range from 20 to 440 days. The results of our ongoing radial velocity monitoring programme with the Hermes spectrograph reveal an eccentricity range of e = 0.2 to 0.76. As a case study we present a detailed analysis of KIC 5006817, whose rich oscillation spectrum allows for detailed seismic analysis. From seismology we constrain the rotational period of the envelope to be at least 165 d, which is roughly twice the orbital period. The stellar core rotates 13 times faster than the surface. From the spectrum and radial velocities we expect that the Doppler beaming signal should have a maximum amplitude of 300 ppm in the light curve. Fixing the mass and radius to the asteroseismically determined values, we find from our binary modelling a value of the gravity darkening exponent that is significantly larger than expected. Through binary modelling, we determine the mass of the secondary component to be 0.29 ± 0.03 M . Conclusions. For KIC 5006817 we exclude pseudo-synchronous rotation of the red giant with the orbit. The comparison of the results from seismology and modelling of the light curve shows a possible alignment of the rotational and orbital axis at the 2σ level. Red giant eccentric systems could be progenitors of cataclysmic variables and hot subdwarf B stars.

PHOEBE 2 is a Python package for modeling the observables of eclipsing star systems, but until now has focused entirely on the forward-model -that is, generating a synthetic model given fixed values of a large number of parameters describing the system and the observations. The inverse problem, obtaining orbital and stellar parameters given observational data, is more complicated and computationally expensive as it requires generating a large set of forward-models to determine which set of parameters and uncertainties best represent the available observational data. The process of determining the best solution and also of obtaining reliable and robust uncertainties on those parameters often requires the use of multiple algorithms, including both optimizers and samplers. Furthermore, the forward-model of PHOEBE has been designed to be as physically robust as possible, but is computationally expensive compared to other codes. It is useful, therefore, to use whichever code is most efficient given the reasonable assumptions for a specific system, but learning the intricacies of multiple codes presents a barrier to doing this in practice. Here we present the 2.3 release of PHOEBE (publicly available from http://phoebe-project.org) which introduces a general framework for defining and handling distributions on parameters, and utilizing multiple different

Heartbeat stars (HB stars) are a class of eccentric binary stars with close periastron passages. The characteristic photometric HB signal evident in their light curves is produced by a combination of tidal distortion, heating, and Doppler boosting near orbital periastron. Many HB stars continue to oscillate after periastron and along the entire orbit, indicative of the tidal excitation of oscillation modes within one or both stars. These systems are among the most eccentric binaries known, and they constitute astrophysical laboratories for the study of tidal effects. We have undertaken a radial velocity (RV) monitoring campaign of Kepler HB stars in order to measure their orbits. We present our first results here, including a sample of 22 Kepler HB systems, where for 19 of them we obtained the Keplerian orbit and for 3 other systems we did not detect a statistically significant RV variability. Results presented here are based on 218 spectra obtained with the Keck/HIRES spectrograph during the 2015 Kepler observing season, and they have allowed us to obtain the largest sample of HB stars with orbits measured using a single instrument, which roughly doubles the number of HB stars with an RV measured orbit. The 19 systems measured here have orbital periods from 7 to 90 days and eccentricities from 0.2 to 0.9. We show that HB stars draw the upper envelope of the eccentricity-period distribution. Therefore, HB stars likely represent a population of stars currently undergoing high eccentricity migration via tidal orbital circularization, and they will allow for new tests of high eccentricity migration theories.

Binary systems where the axis of rotation (spin) of one or both components is tilted w.r.t. the axis of revolution are called misaligned binary systems. The angle of misalignment, obliquity, has been measured for a handful of stars and extrasolar planets to date. Here we present a mathematical framework for a complete and rigorous treatment of misalignment and introduce an extension to the public PHOEBE code that implements this framework. We discuss misalignment for the Roche geometry and introduce methods for computing stellar shapes, equilibrium (generalized Lagrange) points of the potential and minimal requirements for lobe existence. Efficient parameterization of misalignment is proposed in the plane-of-sky coordinates and implementation details in PHOEBE are given alongside the proof-of-concept toy model, comparison with a known misaligned binary DI Her, and comparison with a misaligned planetary system Kepler-13. We provide important mathematical details of the model in the Appendix. This paper accompanies the release of PHOEBE 2.1, which is available from its website http://phoebe-project.org.