Contact binary stars have been known to have a peculiar and somewhat mysterious hydro- and thermodynamical structure since their discovery, which directly affects the radiation distribution in their atmospheres. Over the past several decades, however, observational data of contact binaries have been modeled through a simplified approach, involving the artificial concatenation of the two components of the contact envelope and populating their respective surfaces with either blackbody atmospheres or plane-parallel model atmospheres of single stars. We show the implications this approach has on the reliability of the system parameter values and propose a method to overcome these issues with a new generalized radiative transfer code, COBAIN (COntact Binary Atmospheres with INterpolation). The basic principles of COBAIN are outlined and their application to different geometries and polytropic stellar structures is discussed. We present initial tests on single non-rotating, uniformly rotating and differentially rotating stars, as well as on simplified polytropic structural models of contact binaries. We briefly discuss the final goal of this ambitious project, which is the computation of model atmosphere tables under the correct assumptions for contact binary stars, to be used in modern binary star analysis codes.

We present the Kepler Object of Interest (KOI) catalog of transiting exoplanets based on searching four years of Kepler time series photometry (Data Release 25, Q1-Q17). The catalog contains 8054 KOIs of which 4034 are planet candidates with periods between 0.25 and 632 days. Of these candidates, 219 are new and include two in multi-planet systems (KOI-82.06 and KOI-2926.05), and ten high-reliability, terrestrial-size, habitable zone candidates. This catalog was created using a tool called the Robovetter which automatically vets the DR25 Threshold Crossing Events (TCEs, Twicken et al. 2016). The Robovetter also vetted simulated data sets and measured how well it was able to separate TCEs caused by noise from those caused by low signal-to-noise transits. We discusses the Robovetter and the metrics it uses to sort TCEs. For orbital periods less than 100 days the Robovetter completeness (the fraction of simulated transits that are determined to be planet candidates) across all observed stars is greater than 85%. For the same period range, the catalog reliability (the fraction of candidates that are not due to instrumental or stellar noise) is greater than 98%. However, for low signal-to-noise candidates between 200 and 500 days around FGK dwarf stars, the Robovetter is 76.7% complete and the catalog is 50.5% reliable. The KOI catalog, the transit fits and all of the simulated data used to characterize this catalog are available at the NASA Exoplanet Archive.

It has long been recognized that the finite speed of light can affect the observed time of an event. For example, as a source moves radially toward or away from an observer, the path length and therefore the light travel time to the observer decreases or increases, causing the event to appear earlier or later than otherwise expected, respectively. This light travel time effect (LTTE) has been applied to transits and eclipses for a variety of purposes, including studies of eclipse timing variations (ETVs) and transit timing variations (TTVs) that reveal the presence of additional bodies in the system. Here we highlight another non-relativistic effect on eclipse or transit times arising from the finite speed of light---caused by an asymmetry in the transverse velocity of the two eclipsing objects, relative to the observer. This asymmetry can be due to a non-unity mass ratio or to the presence of external barycentric motion. Although usually constant, this barycentric and asymmetric transverse velocities (BATV) effect can vary between sequential eclipses if either the path length between the two objects or the barycentric transverse velocity varies in time. We discuss this BATV effect and estimate its magnitude for both time-dependent and time-independent cases. For the time-dependent cases, we consider binaries that experience a change in orbital inclination, eccentric systems with and without apsidal motion, and hierarchical triple systems. We also consider the time-independent case which, by affecting the primary and secondary eclipses differently, can influence the inferred system parameters, such as the orbital eccentricity.

In this work we quantify the performance of $LSST$ on the detection of eclipsing binaries. We use $Kepler$ observed binaries to create a large sample of simulated pseudo-$LSST$ binary light curves. From these light curves, we attempt to recover the known binary signal. The success rate of period recovery from the pseudo-$LSST$ light curves is indicative of $LSST$'s expected performance. Using an off-the-shelf Analysis of Variance (AoV) routine, we successfully recover 71% of the targets in our sample. We examine how the binary period impacts recovery success and see that for periods longer than 10~days the chance of successful binary recovery drops below 50%.

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.

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 21 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 d 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.

Heartbeat stars are eccentric (e > 0.2) ellipsoidal variables whose light curves resemble a cardiogram. We present the observations and corresponding model of KIC 3749404, a highly eccentric (e = 0.66), short period (P = 20.3 d) heartbeat star with tidally induced pulsations. A binary star model was created using PHOEBE, which we modified to include tidally induced pulsations and Doppler boosting. The morphology of the photometric periastron variation (heartbeat) depends strongly on the eccentricity, inclination and argument of periastron. We show that the inclusion of tidally induced pulsations in the model significantly changes the parameter values, specifically the inclination and those parameters dependent on it. Furthermore, we determine the rate of apsidal advance by modelling the periastron variation at the beginning and end of the 4-yr Kepler data set and dividing by the elapsed time. We compare the model with the theoretical expectations for classical and general relativistic apsidal motion and find the observed rate to be two orders of magnitude greater than the theoretical rate. We find that the observed rate cannot be explained by tidally induced pulsations alone and consequently hypothesise the presence of a third body in the system.

In this brief communication we provide the rationale for, and the outcome of the International Astronomical Union (IAU) resolution vote at the XXIX-th General Assembly in Honolulu, Hawaii, in 2015, on recommended nominal conversion constants for selected solar and planetary properties. The problem addressed by the resolution is a lack of established conversion constants between solar and planetary values and SI units: a missing standard has caused a proliferation of solar values (e.g., solar radius, solar irradiance, solar luminosity, solar effective temperature and solar mass parameter) in the literature, with cited solar values typically based on best estimates at the time of paper writing. As precision of observations increases, a set of consistent values becomes increasingly important. To address this, an IAU Working Group on Nominal Units for Stellar and Planetary Astronomy formed in 2011, uniting experts from the solar, stellar, planetary, exoplanetary and fundamental astronomy, as well as from general standards fields to converge on optimal values for nominal conversion constants. The effort resulted in the IAU 2015 Resolution B3, passed at the IAU General Assembly by a large majority. The resolution recommends the use of nominal solar and planetary values, which are by definition exact and are expressed in SI units. These nominal values should be understood as conversion factors only, not as the true solar/planetary properties or current best estimates. Authors and journal editors are urged to join in using the standard values set forth by this resolution in future work and publications to help minimize further confusion.

The Kepler Mission has provided unprecedented, nearly continuous photometric data of $\sim$200,000 objects in the $\sim$105 deg$^{2}$ field of view from the beginning of science operations in May of 2009 until the loss of the second reaction wheel in May of 2013. The Kepler Eclipsing Binary Catalog contains information including but not limited to ephemerides, stellar parameters and analytical approximation fits for every known eclipsing binary system in the Kepler Field of View. Using Target Pixel level data collected from Kepler in conjunction with the Kepler Eclipsing Binary Catalog, we identify false positives among eclipsing binaries, i.e. targets that are not eclipsing binaries themselves, but are instead contaminated by eclipsing binary sources nearby on the sky and show eclipsing binary signatures in their light curves. We present methods for identifying these false positives and for extracting new light curves for the true source of the observed binary signal. For each source, we extract three separate light curves for each quarter of available data by optimizing the signal-to-noise ratio, the relative percent eclipse depth and the flux eclipse depth. We present 289 new eclipsing binaries in the Kepler Field of View that were not targets for observation, and these have been added to the Catalog. An online version of this Catalog with downloadable content and visualization tools is maintained at http://keplerEBs.villanova.edu.

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 square degree 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.

The original {\it Kepler} mission observed and characterized over 2400 eclipsing binaries in addition to its prolific exoplanet detections. Despite the mechanical malfunction and subsequent non-recovery of two reaction wheels used to stabilize the instrument, the {\it Kepler} satellite continues collecting data in its repurposed {\it K2} mission surveying a series of fields along the ecliptic plane. Here we present an analysis of the first full baseline {\it K2} data release: the Campaign 0 data-set. In the 7761 light curves, we have identified a total of 207 eclipsing binaries. Of these, 97 are new discoveries that were not previously identified. Our pixel-level analysis of these objects has also resulted in identification of several false positives (observed targets contaminated by neighboring eclipsing binaries), as well as the serendipitous discovery of two short period exoplanet candidates. We provide catalog cross-matched source identifications, orbital periods, morphologies and ephemerides for these eclipsing systems. We also describe the incorporation of the K2 sample into the Kepler Eclipsing Binary Catalog\footnote{\url{keplerebs.villanova.edu/k2}}, present spectroscopic follow-up observations for a limited selection of nine systems, and discuss prospects for upcoming {\it K2} campaigns.

KIC 2835289 is a triple stellar system that consists of an inner ellipsoidal variable with an orbital period of ~0.86 days and an outer star that eclipses the inner pair every ~750 days. Two eclipse events were observed by the Kepler mission, but they do not fully constrain our photodynamical models. The next eclipse event will occur on May 14 UT, and we solicit community involvement to follow it up from the ground. All details are available in the attached call. Please contact the authors to join the follow-up campaign.

KIC 8569819 is an eclipsing binary star with an early F primary and G secondary in a 20.85-d eccentric orbit. The primary is a delta Sct - gamma Dor star pulsating in both p modes and g modes. Using 4 years of Kepler Mission photometric data we independently model the light curve using the traditional technique with the modelling code PHOEBE, and we study the orbital characteristics using the new frequency modulation (FM) technique. We show that both methods provide the equivalent orbital period, eccentricity and argument of periastron, thus illustrating and validating the FM technique. In the amplitude spectrum of the p-mode pulsations we also discovered an FM signal compatible with a third body in the system, a low mass M dwarf in an 861-d orbit around the primary pair. However, the eclipses show no timing variations, indicating that the FM signal is a consequence of the intrinsic change in pulsation frequency, thus providing a cautionary tale. Our analysis shows the potential of the FM technique using {\it Kepler} data, and we discuss the prospects to detect planets and brown dwarfs in Kepler data for A and F stars even in the absence of transits and with no spectroscopic radial velocity curves. This opens the possibility of finding planets orbiting hotter stars that cannot be found by traditional techniques.

Over 2500 eclipsing binaries were identified and characterized from the ultra-precise photometric data provided by the Kepler space telescope. Kepler is now beginning its second mission, K2, which is proving to again provide ultra-precise photometry for a large sample of eclipsing binary stars. In the 1951 light curves covering 12 days in the K2 engineering data-set, we have identified and determined the ephemerides for 31 eclipsing binaries that demonstrate the capabilities for eclipsing binary science in the upcoming campaigns in K2. Of those, 20 are new discoveries. We describe both manual and automated approaches to harvesting the complete set of eclipsing binaries in the K2 data, provide identifications and details for the full set of eclipsing binaries present in the engineering data-set, and discuss the prospects for application of eclipsing binary searches in the K2 mission.

We report on the masses, sizes, and orbits of the planets orbiting 22 Kepler stars. There are 49 planet candidates around these stars, including 42 detected through transits and 7 revealed by precise Doppler measurements of the host stars. Based on an analysis of the Kepler brightness measurements, along with high-resolution imaging and spectroscopy, Doppler spectroscopy, and (for 11 stars) asteroseismology, we establish low false-positive probabilities for all of the transiting planets (41 of 42 have a false-positive probability under 1%), and we constrain their sizes and masses. Most of the transiting planets are smaller than 3X the size of Earth. For 16 planets, the Doppler signal was securely detected, providing a direct measurement of the planet's mass. For the other 26 planets we provide either marginal mass measurements or upper limits to their masses and densities; in many cases we can rule out a rocky composition. We identify 6 planets with densities above 5 g/cc, suggesting a mostly rocky interior for them. Indeed, the only planets that are compatible with a purely rocky composition are smaller than ~2 R_earth. Larger planets evidently contain a larger fraction of low-density material (H, He, and H2O).

We present a catalog of precise eclipse times and analysis of third body signals among 1279 close binaries in the latest Kepler Eclipsing Binary Catalog. For these short period binaries, Kepler's 30 minute exposure time causes significant smearing of light curves. In addition, common astrophysical phenomena such as chromospheric activity, as well as imperfections in the light curve detrending process, can create systematic artifacts that may produce fictitious signals in the eclipse timings. We present a method to measure precise eclipse times in the presence of distorted light curves, such as in contact and near-contact binaries which exhibit continuously changing light levels in and out of eclipse. 236 systems for which we find a timing variation signal compatible with the presence of a third body are identified. These are modeled for the light time travel effect and the basic properties of the third body are derived. This study complements Orosz et al. (2013; in prep), which focuses on eclipse timing variations of longer period binaries with flat out-of-eclipse regions. Together, these two papers provide comprehensive eclipse timings for all binaries in the Kepler Eclipsing Binary Catalog, as an ongoing resource freely accessible online to the community.

We show that the Kepler spacecraft in two-reaction wheel mode of operation is very well suited for the study of eclipsing binary star systems. Continued observations of the Kepler field will provide the most enduring and long-term valuable science. It will enable the discovery and characterization of eclipsing binaries with periods greater than 1 year - these are the most important, yet least understood binaries for habitable-zone planet background considerations. The continued mission will also enable the investigation of hierarchical multiple systems (discovered through eclipse timing variations), and provide drastically improved orbital parameters for circumbinary planetary systems.

We discuss the discovery and characterization of the circumbinary planet Kepler-38b. The stellar binary is single-lined, with a period of 18.8 days, and consists of a moderately evolved main-sequence star (M_A = 0.949 +/- 0.059 solar masses and R_A = 1.757 +/- 0.034 solar radii) paired with a low-mass star (M_B = 0.249 +/- 0.010 solar masses and R_B = 0.2724 +/- 0.0053 solar radii) in a mildly eccentric (e=0.103) orbit. A total of eight transits due to a circumbinary planet crossing the primary star were identified in the Kepler light curve (using Kepler Quarters 1 through 11), from which a planetary period of 105.595 +/- 0.053 days can be established. A photometric dynamical model fit to the radial velocity curve and Kepler light curve yields a planetary radius of 4.35 +/- 0.11 Earth radii, or 1.12 +/- 0.03 Neptune radii. Since the planet is not sufficiently massive to observably alter the orbit of the binary from Keplerian motion, we can only place an upper limit on the mass of the planet of 122 Earth masses (7.11 Neptune masses or 0.384 Jupiter masses) at 95% confidence. This upper limit should decrease as more Kepler data become available.

Most Sun-like stars in the Galaxy reside in gravitationally-bound pairs of stars called "binary stars". While long anticipated, the existence of a "circumbinary planet" orbiting such a pair of normal stars was not definitively established until the discovery of Kepler-16. Incontrovertible evidence was provided by the miniature eclipses ("transits") of the stars by the planet. However, questions remain about the prevalence of circumbinary planets and their range of orbital and physical properties. Here we present two additional transiting circumbinary planets, Kepler-34 and Kepler-35. Each is a low-density gas giant planet on an orbit closely aligned with that of its parent stars. Kepler-34 orbits two Sun-like stars every 289 days, while Kepler-35 orbits a pair of smaller stars (89% and 81% of the Sun's mass) every 131 days. Due to the orbital motion of the stars, the planets experience large multi-periodic variations in incident stellar radiation. The observed rate of circumbinary planets implies > ~1% of close binary stars have giant planets in nearly coplanar orbits, yielding a Galactic population of at least several million.

We present an automated classification of 2165 \textit{Kepler} eclipsing binary (EB) light curves that accompanied the second \textit{Kepler} data release. The light curves are classified using Locally Linear Embedding, a general nonlinear dimensionality reduction tool, into morphology types (detached, semi-detached, overcontact, ellipsoidal). The method, related to a more widely used Principal Component Analysis, produces a lower-dimensional representation of the input data while preserving local geometry and, consequently, the similarity between neighboring data points. We use this property to reduce the dimensionality in a series of steps to a one-dimensional manifold and classify light curves with a single parameter that is a measure of "detachedness" of the system. This fully automated classification correlates well with the manual determination of morphology from the data release, and also efficiently highlights any misclassified objects. Once a lower-dimensional projection space is defined, the classification of additional light curves runs in a negligible time and the method can therefore be used as a fully automated classifier in pipeline structures. The classifier forms a tier of the \textit{Kepler} EB pipeline that pre-processes light curves for the artificial intelligence based parameter estimator.

A search of the time-series photometry from NASA's Kepler spacecraft reveals a transiting planet candidate orbiting the 11th magnitude G5 dwarf KIC 10593626 with a period of 290 days. The characteristics of the host star are well constrained by high-resolution spectroscopy combined with an asteroseismic analysis of the Kepler photometry, leading to an estimated mass and radius of 0.970 +/- 0.060 MSun and 0.979 +/- 0.020 RSun. The depth of 492 +/- 10ppm for the three observed transits yields a radius of 2.38 +/- 0.13 REarth for the planet. The system passes a battery of tests for false positives, including reconnaissance spectroscopy, high-resolution imaging, and centroid motion. A full BLENDER analysis provides further validation of the planet interpretation by showing that contamination of the target by an eclipsing system would rarely mimic the observed shape of the transits. The final validation of the planet is provided by 16 radial velocities obtained with HIRES on Keck 1 over a one year span. Although the velocities do not lead to a reliable orbit and mass determination, they are able to constrain the mass to a 3{\sigma} upper limit of 124 MEarth, safely in the regime of planetary masses, thus earning the designation Kepler-22b. The radiative equilibrium temperature is 262K for a planet in Kepler-22b's orbit. Although there is no evidence that Kepler-22b is a rocky planet, it is the first confirmed planet with a measured radius to orbit in the Habitable Zone of any star other than the Sun.

Planet Hunters is a new citizen science project, designed to engage the public in an exoplanet search using NASA Kepler public release data. In the first month after launch, users identified two new planet candidates which survived our checks for false- positives. The follow-up effort included analysis of Keck HIRES spectra of the host stars, analysis of pixel centroid offsets in the Kepler data and adaptive optics imaging at Keck using NIRC2. Spectral synthesis modeling coupled with stellar evolutionary models yields a stellar density distribution, which is used to model the transit orbit. The orbital periods of the planet candidates are 9.8844 \pm0.0087 days (KIC 10905746) and 49.7696 \pm0.00039 (KIC 6185331) days and the modeled planet radii are 2.65 and 8.05 R\oplus. The involvement of citizen scientists as part of Planet Hunters is therefore shown to be a valuable and reliable tool in exoplanet detection.

Kepler-16 is an eccentric low-mass eclipsing binary with a circumbinary transiting planet. Here we investigate the angular momentum of the primary star, based on Kepler photometry and Keck spectroscopy. The primary star's rotation period is 35.1 +/- 1.0 days, and its projected obliquity with respect to the stellar binary orbit is 1.6 +/- 2.4 degrees. Therefore the three largest sources of angular momentum---the stellar orbit, the planetary orbit, and the primary's rotation---are all closely aligned. This finding supports a formation scenario involving accretion from a single disk. Alternatively, tides may have realigned the stars despite their relatively wide separation (0.2 AU), a hypothesis that is supported by the agreement between the measured rotation period and the "pseudosynchronous" period of tidal evolution theory. The rotation period, chromospheric activity level, and fractional light variations suggest a main-sequence age of 2-4 Gyr. Evolutionary models of low-mass stars can match the observed masses and radii of the primary and secondary stars to within about 3%.

We report the detection of a planet whose orbit surrounds a pair of low-mass stars. Data from the Kepler spacecraft reveal transits of the planet across both stars, in addition to the mutual eclipses of the stars, giving precise constraints on the absolute dimensions of all three bodies. The planet is comparable to Saturn in mass and size, and is on a nearly circular 229-day orbit around its two parent stars. The eclipsing stars are 20% and 69% as massive as the sun, and have an eccentric 41-day orbit. The motions of all three bodies are confined to within 0.5 degree of a single plane, suggesting that the planet formed within a circumbinary disk.

Kepler observations of the star HD 187091 (KID 8112039, hereafter KOI-54) revealed a remarkable light curve exhibiting sharp periodic brightening events every 41.8 days with a superimposed set of oscillations forming a beating pattern in phase with the brightenings. Spectroscopic observations revealed that this is a binary star with a highly eccentric orbit, e=0.83. We are able to match the Kepler light curve and radial velocities with a nearly face-on (i=5.5 degree) binary star model in which the brightening events are caused by tidal distortion and irradiation of nearly identical A stars during their close periastron passage. The two dominant oscillations in the light curve, responsible for the beating pattern, have frequencies that are the 91st and 90th harmonic of the orbital frequency. The power spectrum of the light curve, after removing the binary star brightening component, reveals a large number of pulsations, 30 of which have a signal-to-noise ratio > 7. Nearly all of these pulsations have frequencies that are either integer multiples of the orbital frequency or are tidally-split multiples of the orbital frequency. This pattern of frequencies unambiguously establishes the pulsations as resonances between the dynamic tides at periastron and the free oscillation modes of one of the stars. KOI-54 is only the 4th star to show such a phenomenon, and is by far the richest in terms of excited modes.