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 estimation, optimization, and sampling algorithms. The presented framework supports multiple forward-models, including the robust model built into PHOEBE itself.

Traditionally, the effects of interstellar extinction on binary star light curves have been treated as a uniform reduction in the observed brightness of the system that is independent of orbital phase. However, unless the orbital plane of the system coincides with the plane of the sky, or if the two stars are completely identical and present with minimal mutual irradiation and tidal/rotational distortions, then this is unlikely to be an accurate representation of the effect of interstellar extinction. Here, we present an updated treatment of interstellar extinction as incorporated in the PHOEBE 2.2 release (publicly available from http://phoebe-project.org) and assess the importance of using such an approach in the modeling of different types of binary systems. We also present the incorporation of PHOENIX model atmospheres into the PHOEBE 2.2 release, providing increased fidelity on computed observables down to lower temperatures than previously available. The importance of these new code developments is then highlighted via an extincted toy model of the eclipsing white-dwarf-subdwarf binary SDSS~J235524.29+044855.7 -- demonstrating that, in the age of LSST as well as complementary space-based photometric missions, a proper accounting for extinction and as well as the use of realistic model atmospheres will be essential in deriving accurate binary parameters.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

Theoretical scaling relations for solar-like oscillators and red giants are widely used to estimate fundamental stellar parameters. The accuracy and precision of these relations have often been questioned in the literature, with studies often utilising binarity for model-independent validation. However, it has not been tested if the photometric effects of binarity introduce a systematic effect on the extraction of the seismic properties of the pulsating component(s). In this paper, we present an estimation of the impact of a contaminating photometric signal with a distinct background profile on the global asteroseismic parameter $\nu_{\text{max}}$ through the analysis of synthetic red-giant binary light curves. We generate the pulsational and granulation parameters for single red giants with different masses, radii and effective temperatures from theoretical scaling relations and use them to simulate single red-giant light curves with the characteristics of Kepler long-cadence photometric data. These are subsequently blended together according to their light ratio to generate binary red-giant light curves of various configurations. We then perform a differential analysis to characterise the systematic effects of binarity on the extraction of $\nu_{\text{max}}$. We find that the $\nu_{\text{max}}$ extraction for red-giant power spectra featuring overlapping power excesses is unreliable if unconstrained priors are used. Outside of this scenario, we obtain results that are nearly identical to single-star case. We conclude that i) the photometric effects of binarity on the extraction of $\nu_{\text{max}}$ are largely negligible as long as the power excesses of the individual components do not overlap, and that ii) there is minimal advantage to using more than two super-Lorentzian components to model the granulation signal of a binary red-giant.

The Maunakea Spectroscopic Explorer (MSE) is a planned 11.25-m aperture facility with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. A rebirth of the 3.6m Canada-France-Hawaii Telescope on Maunakea, MSE will use 4332 fibers operating at three different resolving powers (R ~ 2500, 6000, 40000) across a wavelength range of 0.36-1.8mum, with dynamical fiber positioning that allows fibers to match the exposure times of individual objects. MSE will enable spectroscopic surveys with unprecedented scale and sensitivity by collecting millions of spectra per year down to limiting magnitudes of g ~ 20-24 mag, with a nominal velocity precision of ~100 m/s in high-resolution mode. This white paper describes science cases for stellar astrophysics and exoplanet science using MSE, including the discovery and atmospheric characterization of exoplanets and substellar objects, stellar physics with star clusters, asteroseismology of solar-like oscillators and opacity-driven pulsators, studies of stellar rotation, activity, and multiplicity, as well as the chemical characterization of AGB and extremely metal-poor stars.

Light curves of astrophysical objects frequently contain strictly periodic signals. In those cases we can use that property to aid the detrending algorithm to fully disentangle an unknown periodic signal and an unknown baseline signal with no power at that period. The periodic signal is modeled as a discrete probability distribution function (pdf), while the baseline signal is modeled as a residual timeseries. Those two components are disentangled by minimizing the length of the residual timeseries w.r.t. the per-bin pdf fluxes. We demonstrate the use of the algorithm on a synthetic case, on the eclipsing binary KIC 3953981 and on the eccentric ellipsoidal variable KIC 3547874. We further discuss the parameters and the limitations of the algorithm and speculate on the two most common use cases: detrending the periodic signal of interest and measuring the dependence of instrumental response on controlled instrumental variables. A more sophisticated version of the algorithm is released as open source on github and available via pip.

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.

The advent of large scale multi-epoch surveys raises the need for automated light curve (LC) processing. This is particularly true for eclipsing binaries (EBs), which form one of the most populated types of variable objects. The Gaia mission, launched at the end of 2013, is expected to detect of the order of few million EBs over a 5-year mission. We present an automated procedure to characterize EBs based on the geometric morphology of their LCs with two aims: first to study an ensemble of EBs on a statistical ground without the need to model the binary system, and second to enable the automated identification of EBs that display atypical LCs. We model the folded LC geometry of EBs using up to two Gaussian functions for the eclipses and a cosine function for any ellipsoidal-like variability that may be present between the eclipses. The procedure is applied to the OGLE-III data set of EBs in the Large Magellanic Cloud (LMC) as a proof of concept. The bayesian information criterion is used to select the best model among models containing various combinations of those components, as well as to estimate the significance of the components. Based on the two-Gaussian models, EBs with atypical LC geometries are successfully identified in two diagrams, using the Abbe values of the original and residual folded LCs, and the reduced $\chi^2$. Cleaning the data set from the atypical cases and further filtering out LCs that contain non-significant eclipse candidates, the ensemble of EBs can be studied on a statistical ground using the two-Gaussian model parameters. For illustration purposes, we present the distribution of projected eccentricities as a function of orbital period for the OGLE-III set of EBs in the LMC, as well as the distribution of their primary versus secondary eclipse widths.

In the new era of large-scale astronomical surveys, automated methods of analysis and classification of bulk data are a fundamental tool for fast and efficient production of deliverables. This becomes ever more imminent as we enter the Gaia era. We investigate the potential detectability of eclipsing binaries with Gaia using a data set of all Kepler eclipsing binaries sampled with Gaia cadence and folded with the Kepler period. The performance of fitting methods is evaluated with comparison to real Kepler data parameters and a classification scheme is proposed for the potentially detectable sources based on the geometry of the light curve fits. The polynomial chain (polyfit) and two-Gaussian models are used for light curve fitting of the data set. Classification is performed with a combination of the t-SNE (t-distrubuted Stochastic Neighbor Embedding) and DBSCAN (Density-Based Spatial Clustering of Applications with Noise) algorithms. We find that approximately 68% of Kepler Eclipsing Binary sources are potentially detectable by Gaia when folded with the Kepler period and propose a classification scheme of the detectable sources based on the morphological type indicative of the light curve, with subclasses that reflect the properties of the fitted model (presence and visibility of eclipses, their width, depth, etc.).

The ESA Gaia mission provides a unique time-domain survey for more than one billion sources brighter than G=20.7 mag. Gaia offers the unprecedented opportunity to study variability phenomena in the Universe thanks to multi-epoch G-magnitude photometry in addition to astrometry, blue and red spectro-photometry, and spectroscopy. Within the Gaia Consortium, Coordination Unit 7 has the responsibility to detect variable objects, classify them, derive characteristic parameters for specific variability classes, and provide global descriptions of variable phenomena. We describe the variability processing and analysis that we plan to apply to the successive data releases, and we present its application to the G-band photometry results of the first 14 months of Gaia operations that comprises 28 days of Ecliptic Pole Scanning Law and 13 months of Nominal Scanning Law. Out of the 694 million, all-sky, sources that have calibrated G-band photometry in this first stage of the mission, about 2.3 million sources that have at least 20 observations are located within 38 degrees from the South Ecliptic Pole. We detect about 14% of them as variable candidates, among which the automated classification identified 9347 Cepheid and RR Lyrae candidates. Additional visual inspections and selection criteria led to the publication of 3194 Cepheid and RR Lyrae stars, described in Clementini et al. (2016). Under the restrictive conditions for DR1, the completenesses of Cepheids and RR Lyrae stars are estimated at 67% and 58%, respectively, numbers that will significantly increase with subsequent Gaia data releases. Data processing within the Gaia Consortium is iterative, the quality of the data and the results being improved at each iteration. The results presented in this article show a glimpse of the exceptional harvest that is to be expected from the Gaia mission for variability phenomena. [abridged]

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 XXIXth IAU General Assembly in Honolulu adopted IAU 2015 Resolution B2 on recommended zero points for the absolute and apparent bolometric magnitude scales. The resolution was proposed by the IAU Inter-Division A-G Working Group on Nominal Units for Stellar and Planetary Astronomy after consulting with a broad spectrum of researchers from the astronomical community. Resolution B2 resolves the long-standing absence of an internationally-adopted zero point for the absolute and apparent bolometric magnitude scales. Resolution B2 defines the zero point of the absolute bolometric magnitude scale such that a radiation source with $M_{\rm Bol}$ = 0 has luminosity L$_{\circ}$ = 3.0128e28 W. The zero point of the apparent bolometric magnitude scale ($m_{\rm Bol}$ = 0) corresponds to irradiance $f_{\circ}$ = 2.518021002e-8 W/m$^2$. The zero points were chosen so that the nominal solar luminosity (3.828e26 W) adopted by IAU 2015 Resolution B3 corresponds approximately to $M_{\rm Bol}$(Sun) = 4.74, the value most commonly adopted in recent literature. The nominal total solar irradiance (1361 W/m$^2$) adopted in IAU 2015 Resolution B3 corresponds approximately to apparent bolometric magnitude $m_{\rm bol}$(Sun) = -26.832. Implicit in the IAU 2015 Resolution B2 definition of the apparent bolometric magnitude scale is an exact definition for the parsec (648000/$\pi$ au) based on the IAU 2012 Resolution B2 definition of the astronomical unit.

Astronomers commonly quote the properties of celestial objects in units of parameters for the Sun, Jupiter, or the Earth. The resolution presented here was proposed by the IAU Inter-Division Working Group on Nominal Units for Stellar and Planetary Astronomy and passed by the XXIXth IAU General Assembly in Honolulu. IAU 2015 Resolution B3 adopts a set of nominal solar, terrestrial, and jovian conversion constants for stellar and (exo)planetary astronomy which are defined to be exact SI values. While the nominal constants are based on current best estimates (CBEs; which have uncertainties, are not secularly constant, and are updated regularly using new observations), they should be interpreted as standard values and not as CBEs. IAU 2015 Resolution B3 adopts five solar conversion constants (nominal solar radius, nominal total solar irradiance, nominal solar luminosity, nominal solar effective temperature, and nominal solar mass parameter) and six planetary conversion constants (nominal terrestrial equatorial radius, nominal terrestrial polar radius, nominal jovian equatorial radius, nominal jovian polar radius, nominal terrestrial mass parameter, and nominal jovian mass parameter).

Gamma Doradus and delta Scuti pulsators cover the transition region between low mass and massive main-sequence stars, and as such, are critical for testing stellar models. When they reside in binary systems, we can combine two independent methods to derive critical information, such as precise fundamental parameters to aid asteroseismic modelling. In the Kepler light curve of KIC10080943, clear signatures of gravity- and pressure-mode pulsations have been found. Ground-based spectroscopy revealed this target to be a double-lined binary system. We present the analysis of four years of Kepler photometry and high-resolution spectroscopy to derive observational constraints with which to evaluate theoretical predictions of the stellar structure and evolution for intermediate-mass stars. We used the method of spectral disentangling to determine atmospheric parameters for both components and derive the orbital elements. With PHOEBE, we modelled the ellipsoidal variation and reflection signal of the binary in the light curve and used classical Fourier techniques to analyse the pulsation modes. We show that the eccentric binary system KIC10080943 contains two hybrid pulsators with masses $M_1=2.0\pm0.1~M_\odot$ and $M_2=1.9\pm0.1~M_\odot$, with radii $R_1=2.9\pm0.1~R_\odot$ and $R_2=2.1\pm0.2~R_\odot$. We detect rotational splitting in the g modes and p modes for both stars and use them to determine a first rough estimate of the core-to-surface rotation rates for the two components, which will be improved by future detailed seismic modelling.

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.

We report the discovery of two new low-mass, thermally bloated, hot white dwarfs among the Kepler sample of eclipsing binaries. These are KIC 9164561 and KIC 10727668 with orbital periods of 1.2670 and 2.3058 days, respectively. The current primary in both systems is an A star of about 2 Msun. This brings the number of similar binaries among the Kepler sample to six, and the two new systems have the shortest orbital periods among them. The white dwarf in KIC 9164561 has the largest thermal bloating, compared to its cold degenerate radius, of about a factor of 14. We utilize RV measurements of the A star in KIC 9164561 to determine the white dwarf mass rather accurately: 0.197 +/- 0.005 Msun. The mass of the white dwarf in KIC 10727668 is based on the Doppler boosting signal in the Kepler photometry, and is less accurately determined to be 0.266 +/- 0.035 Msun. Based on the inferred radii and effective temperatures of these two white dwarfs we are able to make an independent theoretical estimate of their masses to within ~0.01 Msun based on evolutionary models of their cooling history after they lose their hydrogen-rich envelopes. We also present evidence that there is a third body in the KIC 9164561 system with an orbital period of 8-14 years.

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.

KIC 8560861 (HD 183648) is a marginally eccentric (e=0.05) eclipsing binary with an orbital period of P_orb=31.973d, exhibiting mmag amplitude pulsations on time scales of a few days. We present the results of the complex analysis of high and medium-resolution spectroscopic data and Kepler Q0 -- Q16 long cadence photometry. The iterative combination of spectral disentangling, atmospheric analysis, radial velocity and eclipse timing variation studies, separation of pulsational features of the light curve, and binary light curve analysis led to the accurate determination of the fundamental stellar parameters. We found that the binary is composed of two main sequence stars with an age of 0.9\+-0.2 Gyr, having masses, radii and temperatures of M_1=1.93+-0.12 M_sun, R_1=3.30+-0.07 R_sun, T_eff1=7650+-100 K for the primary, and M_2=1.06+-0.08 M_sun, R_2=1.11+-0.03 R_sun, T_eff2=6450+-100 K for the secondary. After subtracting the binary model, we found three independent frequencies, two of which are separated by twice the orbital frequency. We also found an enigmatic half orbital period sinusoidal variation that we attribute to an anomalous ellipsoidal effect. Both of these observations indicate that tidal effects are strongly influencing the luminosity variations of HD 183648. The analysis of the eclipse timing variations revealed both a parabolic trend, and apsidal motion with a period of (P_apse)_obs=10,400+-3,000 y, which is three times faster than what is theoretically expected. These findings might indicate the presence of a distant, unseen companion.

K2 is the mission concept for a repurposed Kepler mission that uses two reaction wheels to maintain the satellite attitude and provide ~85 days of coverage for ten 105-deg^2 fields along the ecliptic in the first 2.5 years of operation. We examine stellar populations based on the improved Besancon model of the Galaxy, comment on the general properties for the entire ecliptic plane, and provide stellar occurrence rates in the first 6 tentative K2 campaigns grouped by spectral type and luminosity class. For each campaign we distinguish between main sequence stars and giants and provide their density profile as a function of galactic latitude. We introduce the crowding metric that serves for optimized target selection across the campaigns. For all main sequence stars we compute the expected planetary occurrence rates for three planet sizes: 2-4 REarth, 4-8 REarth, and 8-32 REarth with orbital periods up to 50 days. In conjunction with Gaia and the upcoming TESS and Plato missions, K2 will become a gold mine for stellar and planetary astrophysics.

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 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 delta Scuti and gamma 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 MSun and 1.60+/-0.06 MSun, and radii, 1.76+/-0.03 RSun and 1.42+/-0.02 RSun, 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 in fact 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.

We present the detection of five planets -- Kepler-62b, c, d, e, and f -- of size 1.31, 0.54, 1.95, 1.61 and 1.41 Earth radii, orbiting a K2V star at periods of 5.7, 12.4, 18.2, 122.4 and 267.3 days, respectively. The outermost planets (Kepler-62e & -62f) are super-Earth-size (1.25 < planet radius/earth radius < 2.0) planets in the habitable zone (HZ) of their host star, receiving 1.2 +- 0.2 and 0.41 +- 0.05 times the solar flux at Earth's orbit. Theoretical models of Kepler-62e and -62f for a stellar age of ~7 Gyr suggest that both planets could be solid: either with a rocky composition or composed of mostly solid water in their bulk.

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.

We analyze the deviations of transit times from a linear ephemeris for the Kepler Objects of Interest (KOI) through Quarter six (Q6) of science data. We conduct two statistical tests for all KOIs and a related statistical test for all pairs of KOIs in multi-transiting systems. These tests identify several systems which show potentially interesting transit timing variations (TTVs). Strong TTV systems have been valuable for the confirmation of planets and their mass measurements. Many of the systems identified in this study should prove fruitful for detailed TTV studies.

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.

New transiting planet candidates are identified in sixteen months (May 2009 - September 2010) of data from the Kepler spacecraft. Nearly five thousand periodic transit-like signals are vetted against astrophysical and instrumental false positives yielding 1,091 viable new planet candidates, bringing the total count up to over 2,300. Improved vetting metrics are employed, contributing to higher catalog reliability. Most notable is the noise-weighted robust averaging of multi-quarter photo-center offsets derived from difference image analysis which identifies likely background eclipsing binaries. Twenty-two months of photometry are used for the purpose of characterizing each of the new candidates. Ephemerides (transit epoch, T_0, and orbital period, P) are tabulated as well as the products of light curve modeling: reduced radius (Rp/R*), reduced semi-major axis (d/R*), and impact parameter (b). The largest fractional increases are seen for the smallest planet candidates (197% for candidates smaller than 2Re compared to 52% for candidates larger than 2Re) and those at longer orbital periods (123% for candidates outside of 50-day orbits versus 85% for candidates inside of 50-day orbits). The gains are larger than expected from increasing the observing window from thirteen months (Quarter 1-- Quarter 5) to sixteen months (Quarter 1 -- Quarter 6). This demonstrates the benefit of continued development of pipeline analysis software. The fraction of all host stars with multiple candidates has grown from 17% to 20%, and the paucity of short-period giant planets in multiple systems is still evident. The progression toward smaller planets at longer orbital periods with each new catalog release suggests that Earth-size planets in the Habitable Zone are forthcoming if, indeed, such planets are abundant.

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.

We present a hierarchical triple star system (KIC 9140402) where a low mass eclipsing binary orbits a more massive third star. The orbital period of the binary (4.98829 Days) is determined by the eclipse times seen in photometry from NASA's Kepler spacecraft. The periodically changing tidal field, due to the eccentric orbit of the binary about the tertiary, causes a change in the orbital period of the binary. The resulting eclipse timing variations provide insight into the dynamics and architecture of this system and allow the inference of the total mass of the binary ($0.424 \pm 0.017 \text{M}_\odot$) and the orbital parameters of the binary about the central star.

In this paper we estimate the Large Synoptic Survey Telescope (LSST) yield of eclipsing binary stars, which will survey ~20,000 square degrees of the southern sky during the period of 10 years in 6 photometric passbands to r ~ 24.5. We generate a set of 10,000 eclipsing binary light curves sampled to the LSST time cadence across the whole sky, with added noise as a function of apparent magnitude. This set is passed to the Analysis of Variance (AoV) period finder to assess the recoverability rate for the periods, and the successfully phased light curves are passed to the artificial intelligence-based pipeline EBAI to assess the recoverability rate in terms of the eclipsing binaries' physical and geometric parameters. We find that, out of ~24 million eclipsing binaries observed by LSST with S/N>10 in mission life-time, ~28% or 6.7 million can be fully characterized by the pipeline. Of those, ~25% or 1.7 million will be double-lined binaries, a true treasure trove for stellar astrophysics.

We report the distribution of planets as a function of planet radius (R_p), orbital period (P), and stellar effective temperature (Teff) for P < 50 day orbits around GK stars. These results are based on the 1,235 planets (formally "planet candidates") from the Kepler mission that include a nearly complete set of detected planets as small as 2 Earth radii (Re). For each of the 156,000 target stars we assess the detectability of planets as a function of R_p and P. We also correct for the geometric probability of transit, R*/a. We consider first stars within the "solar subset" having Teff = 4100-6100 K, logg = 4.0-4.9, and Kepler magnitude Kp < 15 mag. We include only those stars having noise low enough to permit detection of planets down to 2 Re. We count planets in small domains of R_p and P and divide by the included target stars to calculate planet occurrence in each domain. Occurrence of planets varies by more than three orders of magnitude and increases substantially down to the smallest radius (2 Re) and out to the longest orbital period (50 days, ~0.25 AU) in our study. For P < 50 days, the radius distribution is given by a power law, df/dlogR= k R^\alpha. This rapid increase in planet occurrence with decreasing planet size agrees with core-accretion, but disagrees with population synthesis models. We fit occurrence as a function of P to a power law model with an exponential cutoff below a critical period P_0. For smaller planets, P_0 has larger values, suggesting that the "parking distance" for migrating planets moves outward with decreasing planet size. We also measured planet occurrence over Teff = 3600-7100 K, spanning M0 to F2 dwarfs. The occurrence of 2-4 Re planets in the Kepler field increases with decreasing Teff, making these small planets seven times more abundant around cool stars than the hottest stars in our sample. [abridged]

On 1 February 2011 the Kepler Mission released data for 156,453 stars observed from the beginning of the science observations on 2 May through 16 September 2009. There are 1235 planetary candidates with transit like signatures detected in this period. These are associated with 997 host stars. Distributions of the characteristics of the planetary candidates are separated into five class-sizes; 68 candidates of approximately Earth-size (radius < 1.25 Earth radii), 288 super-Earth size (1.25 Earth radii < radius < 2 Earth radii), 662 Neptune-size (2 Earth radii < radius < 6 Earth radii), 165 Jupiter-size (6 Earth radii < radius < 15 Earth radii), and 19 up to twice the size of Jupiter (15 Earth radii < radius < 22 Earth radii). In the temperature range appropriate for the habitable zone, 54 candidates are found with sizes ranging from Earth-size to larger than that of Jupiter. Five are less than twice the size of the Earth. Over 74% of the planetary candidates are smaller than Neptune. The observed number versus size distribution of planetary candidates increases to a peak at two to three times Earth-size and then declines inversely proportional to area of the candidate. Our current best estimates of the intrinsic frequencies of planetary candidates, after correcting for geometric and sensitivity biases, are 6% for Earth-size candidates, 7% for super-Earth size candidates, 17% for Neptune-size candidates, and 4% for Jupiter-size candidates. Multi-candidate, transiting systems are frequent; 17% of the host stars have multi-candidate systems, and 33.9% of all the candidates are part of multi-candidate systems.

The Kepler Mission provides nearly continuous monitoring of ~156 000 objects with unprecedented photometric precision. Coincident with the first data release, we presented a catalog of 1879 eclipsing binary systems identified within the 115 square degree Kepler FOV. Here, we provide an updated catalog augmented with the second Kepler data release which increases the baseline nearly 4-fold to 125 days. 386 new systems have been added, ephemerides and principle parameters have been recomputed. We have removed 42 previously cataloged systems that are now clearly recognized as short-period pulsating variables and another 58 blended systems where we have determined that the Kepler target object is not itself the eclipsing binary. A number of interesting objects are identified. We present several exemplary cases: 4 EBs that exhibit extra (tertiary) eclipse events; and 8 systems that show clear eclipse timing variations indicative of the presence of additional bodies bound in the system. We have updated the period and galactic latitude distribution diagrams. With these changes, the total number of identified eclipsing binary systems in the Kepler field-of-view has increased to 2165, 1.4% of the Kepler target stars.

We present Kepler satellite photometry of KIC 10661783, a short-period binary star system which shows total eclipses and multi-periodic delta Scuti pulsations. A frequency analysis of the eclipse-subtracted light curve reveals at least 68 frequencies of which 55 or more can be attributed to pulsation modes. The main limitation on this analysis is the frequency resolution within the 27-day short-cadence light curve. Most of the variability signal lies in the frequency range 18 to 31 c/d, with amplitudes between 0.1 and 4 mmag. One harmonic term (2.f) and a few combination frequencies (f_i+f_j) have been detected. From a plot of the residuals versus orbital phase we assign the pulsations to the primary star in the system. The pulsations were removed from the short-cadence data and the light curve was modelled using the Wilson-Devinney code. We are unable to get a perfect fit due to the residual effects of pulsations and also to the treatment of reflection and reprocessing in the light curve model. A model where the secondary star fills its Roche lobe is favoured, which means that KIC 10661783 can be classified as an oEA system. Further photometric and spectroscopic observations will allow the masses and radii of the two stars to be measured to high precision and hundreds of delta Scuti pulsation frequencies to be resolved. This could lead to unique constraints on theoretical models of delta Scuti stars, if the evolutionary history of KIC 10661783 can be accounted for.

The Kepler space mission is devoted to finding Earth-size planets in habitable zones orbiting other stars. Its large, 105-deg field-of-view features over 156,000 stars that are observed continuously to detect and characterize planet transits. Yet this high-precision instrument holds great promise for other types of objects as well. Here we present a comprehensive catalog of eclipsing binary stars observed by Kepler in the first 44 days of operation, the data which are publicly available through MAST as of 6/15/2010. The catalog contains 1879 unique objects. For each object we provide its Kepler ID (KID), ephemeris (BJD0, P0), morphology type, physical parameters (Teff, log g, E(B-V), crowding), and principal parameters (T2/T1, q, fillout factor and sin i for overcontacts, and T2/T1, (R1+R2)/a, e sin(w), e cos(w), and sin i for detached binaries). We present statistics based on the determined periods and measure an average occurence rate of eclipsing binaries to be ~1.2% across the Kepler field. We further discuss the distribution of binaries as function of galactic latitude, and thoroughly explain the application of artificial intelligence to obtain principal parameters in a matter of seconds for the whole sample. The catalog was envisioned to serve as a bridge between the now public Kepler data and the scientific community interested in eclipsing binary stars.

Light curves from the Kepler Mission contain valuable information on the nature of the phenomena producing the transit-like signals. To assist in exploring the possibility that they are due to an astrophysical false positive, we describe a procedure (BLENDER) to model the photometry in terms of a "blend" rather than a planet orbiting a star. A blend may consist of a background or foreground eclipsing binary (or star-planet pair) whose eclipses are attenuated by the light of the candidate and possibly other stars within the photometric aperture. We apply BLENDER to the case of Kepler-9, a target harboring two previously confirmed Saturn-size planets (Kepler-9b and Kepler-9c) showing transit timing variations, and an additional shallower signal with a 1.59-day period suggesting the presence of a super-Earth-size planet. Using BLENDER together with constraints from other follow-up observations we are able to rule out all blends for the two deeper signals, and provide independent validation of their planetary nature. For the shallower signal we rule out a large fraction of the false positives that might mimic the transits. The false alarm rate for remaining blends depends in part (and inversely) on the unknown frequency of small-size planets. Based on several realistic estimates of this frequency we conclude with very high confidence that this small signal is due to a super-Earth-size planet (Kepler-9d) in a multiple system, rather than a false positive. The radius is determined to be 1.64 (+0.19/-0.14) R(Earth), and current spectroscopic observations are as yet insufficient to establish its mass.

Ten days of commissioning data (Quarter 0) and thirty-three days of science data (Quarter 1) yield instrumental flux timeseries of ~150,000 stars that were combed for transit events, termed Threshold Crossing Events (TCE), each having a total detection statistic above 7.1-sigma. TCE light curves are modeled as star+planet systems. Those returning a companion radius smaller than 2R_J are assigned a KOI (Kepler Object of Interest) number. The raw flux, pixel flux, and flux-weighted centroids of every KOI are scrutinized to assess the likelihood of being an astrophysical false-positive versus the likelihood of a being a planetary companion. This vetting using Kepler data is referred to as data validation. Herein, we describe the data validation metrics and graphics used to identify viable planet candidates amongst the KOIs. Light curve modeling tests for a) the difference in depth of the odd- versus even-numbered transits, b) evidence of ellipsoidal variations, and c) evidence of a secondary eclipse event at phase=0.5. Flux-weighted centroids are used to test for signals correlated with transit events with a magnitude and direction indicative of a background eclipsing binary. Centroid timeseries are complimented by analysis of images taken in-transit versus out-of-transit, the difference often revealing the pixel contributing the most to the flux change during transit. Examples are shown to illustrate each test. Candidates passing data validation are submitted to ground-based observers for further false-positive elimination or confirmation/characterization.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

Accurate photometric CoRoT space observations of a secondary seismological target, HD 174884, led to the discovery that this star is an astrophysically important double-lined eclipsing spectroscopic binary in an eccentric orbit (e of about 0.3), unusual for its short (3.65705d) orbital period. The high eccentricity, coupled with the orientation of the binary orbit in space, explains the very unusual observed light curve with strongly unequal primary and secondary eclipses having the depth ratio of 1-to-100 in the CoRoT 'seismo' passband. Without the high accuracy of the CoRoT photometry, the secondary eclipse, 1.5 mmag deep, would have gone unnoticed. A spectroscopic follow-up program provided 45 high dispersion spectra. The analysis of the CoRoT light curve was performed with an adapted version of PHOEBE that supports CoRoT passbands. The final solution was obtained by simultaneous fitting of the light and the radial velocity curves. Individual star spectra were derived by spectrum disentangling. The uncertainties of the fit were derived by bootstrap resampling and the solution uniqueness was tested by heuristic scanning. The results provide a consistent picture of the system composed of two late B stars. The Fourier analysis of the light curve fit residuals yields two components, with orbital frequency multiples and an amplitude of about 0.1 mmag, which are tentatively interpreted as tidally induced pulsations. An extensive comparison with theoretical models is carried out by means of the Levenberg-Marquardt minimization technique and the discrepancy between models and the derived parameters is discussed. The best fitting models yield a young system age of 125 million years which is consistent with the eccentric orbit and synchronous component rotation at periastron.

Analyses of very accurate CoRoT space photometry, past Johnson V photoelectric photometry and high-resolution \'echelle spectra led to the determination of improved and consistent fundamental stellar properties of both components of AU Mon. We derived new, accurate ephemerides for both the orbital motion (with a period of 11.113d) and the long-term, overall brightness variation (with a period of 416.9d) of this strongly interacting Be + G semi-detached binary. It is shown that this long-term variation must be due to attenuation of the total light by some variable circumbinary material. We derived the binary mass ratio $M_{\rm G}/M_{\rm B}$ = 0.17\p0.03 based on the assumption that the G-type secondary fills its Roche lobe and rotates synchronously. Using this value of the mass ratio as well as the radial velocities of the G-star, we obtained a consistent light curve model and improved estimates of the stellar masses, radii, luminosities and effective temperatures. We demonstrate that the observed lines of the B-type primary may not be of photospheric origin. We also discover rapid and periodic light changes visible in the high-quality residual CoRoT light curves. AU Mon is put into perspective by a comparison with known binaries exhibiting long-term cyclic light changes.

We present the JHKs light curves for the double-lined eclipsing binary 2MASS J05352184-0546085, in which both components are brown dwarfs. We analyze these light curves with the published Ic-band light curve and radial velocities to provide refined measurements of the system's physical parameters. The component masses and radii are here determined with an accuracy of ~6.5% and ~1.5%, respectively. We confirm the previous surprising finding that the primary brown dwarf has a cooler effective temperature than its companion. Next, we perform a detailed study of the variations in the out-of-eclipse phases of the light curves to ascertain the properties of any inhomogeneities on the surfaces of the brown dwarfs. Our analysis reveals two low-amplitude periodic signals, one attributable to the rotation of the primary (with a period of 3.293+/-0.001 d) and the other to that of the secondary (14.05+/-0.05 d). Finally, we explore the effects on the derived physical parameters of the system when spots are included in the modeling. The observed low-amplitude rotational modulations are well fit by cool spots covering a small fraction of their surfaces. To mimic the observed ~200 K suppression of the primary's temperature, our model requires that the primary possess a very large spot coverage fraction of ~65%. Altogether, a spot configuration in which the primary is heavily spotted while the secondary is lightly spotted can explain the apparent temperature reversal and can bring the temperatures of the brown dwarfs into agreement with the predictions of theoretical models.

Observational astronomy has changed drastically in the last decade: manually driven target-by-target instruments have been replaced by fully automated robotic telescopes. Data acquisition methods have advanced to the point that terabytes of data are flowing in and being stored on a daily basis. At the same time, the vast majority of analysis tools in stellar astrophysics still rely on manual expert interaction. To bridge this gap, we foresee that the next decade will witness a fundamental shift in the approaches to data analysis: case-by-case methods will be replaced by fully automated pipelines that will process the data from their reduction stage, through analysis, to storage. While major effort has been invested in data reduction automation, automated data analysis has mostly been neglected despite the urgent need. Scientific data mining will face serious challenges to identify, understand and eliminate the sources of systematic errors that will arise from this automation. As a special case, we present an artificial intelligence (AI) driven pipeline that is prototyped in the domain of stellar astrophysics (eclipsing binaries in particular), current results and the challenges still ahead.

The next decade of survey astronomy has the potential to transform our knowledge of variable stars. Stellar variability underpins our knowledge of the cosmological distance ladder, and provides direct tests of stellar formation and evolution theory. Variable stars can also be used to probe the fundamental physics of gravity and degenerate material in ways that are otherwise impossible in the laboratory. The computational and engineering advances of the past decade have made large-scale, time-domain surveys an immediate reality. Some surveys proposed for the next decade promise to gather more data than in the prior cumulative history of astronomy. The actual implementation of these surveys will have broad implications for the types of science that will be enabled. We examine the design considerations for an optimal time-domain photometric survey dedicated to variable star science, including: observing cadence, wavelength coverage, photometric and astrometric accuracy, single-epoch and cumulative depth, overall sky coverage, and data access by the broader astronomical community. The best surveys must combine aspects from each of these considerations to fully realize the potential for the next decade of time-domain science.

Eclipsing binary stars provide highly accurate measurements of the fundamental physical properties of stars. They therefore serve as stringent tests of the predictions of evolutionary models upon which most stellar age determinations are based. Models generally perform very well in predicting coeval ages for eclipsing binaries with main-sequence components more massive than ~1.2 Msun; relative ages are good to ~5% or better in this mass regime. Low-mass main-sequence stars (M < 0.8 Msun) reveal large discrepancies in the model predicted ages, primarily due to magnetic activity in the observed stars that appears to inhibit convection and likely causes the radii to be 10-20% larger than predicted. In mass-radius diagrams these stars thus appear 50-90% older or younger than they really are. Aside from these activity-related effects, low-mass pre--main-sequence stars at ages ~1 Myr can also show non-coevality of ~30% due to star formation effects, however these effects are largely erased after ~10 Myr.

During the OGLE-2 operation, Soszynski et al. (2003) found 3 LMC candidates for an RR Lyr-type component in an eclipsing binary system. Two of those have orbital periods that are too short to be physically plausible and hence have to be optical blends. For the third, OGLE052218.07-692827.4, we developed a model of the binary that could host the observed RR Lyr star. After being granted HST/WFPC2 time, however, we were able to resolve 5 distinct sources within a 1.3" region that is typical of OGLE resolution, proving that OGLE052218.07-692827.4 is also an optical blend. Moreover, the putative eclipsing binary signature found in the OGLE data does not seem to correspond to a physically plausible system; the source is likely another background RR Lyr star. There are still no RR Lyr stars discovered so far in an eclipsing binary system.

Achieving maximum scientific results from the overwhelming volume of astronomical data to be acquired over the next few decades will demand novel, fully automatic methods of data analysis. Artificial intelligence approaches hold great promise in contributing to this goal. Here we apply neural network learning technology to the specific domain of eclipsing binary (EB) stars, of which only some hundreds have been rigorously analyzed, but whose numbers will reach millions in a decade. Well-analyzed EBs are a prime source of astrophysical information whose growth rate is at present limited by the need for human interaction with each EB data-set, principally in determining a starting solution for subsequent rigorous analysis. We describe the artificial neural network (ANN) approach which is able to surmount this human bottleneck and permit EB-based astrophysical information to keep pace with future data rates. The ANN, following training on a sample of 33,235 model light curves, outputs a set of approximate model parameters (T2/T1, (R1+R2)/a, e sin(omega), e cos(omega), and sin i) for each input light curve data-set. The whole sample is processed in just a few seconds on a single 2GHz CPU. The obtained parameters can then be readily passed to sophisticated modeling engines. We also describe a novel method polyfit for pre-processing observational light curves before inputting their data to the ANN and present the results and analysis of testing the approach on synthetic data and on real data including fifty binaries from the Catalog and Atlas of Eclipsing Binaries (CALEB) database and 2580 light curves from OGLE survey data. [abridged]

One of the most important changes in observational astronomy of the 21st Century is a rapid shift from classical object-by-object observations to extensive automatic surveys. As CCD detectors are getting better and their prices are getting lower, more and more small and medium-size observatories are refocusing their attention to detection of stellar variability through systematic sky-scanning missions. This trend is aditionally powered by the success of pioneering surveys such as ASAS, DENIS, OGLE, TASS, their space counterpart Hipparcos and others. Such surveys produce massive amounts of data and it is not at all clear how these data are to be reduced and analysed. This is especially striking in the eclipsing binary (EB) field, where most frequently used tools are optimized for object-by-object analysis. A clear need for thorough, reliable and fully automated approaches to modeling and analysis of EB data is thus obvious. This task is very difficult because of limited data quality, non-uniform phase coverage and solution degeneracy. This paper reviews recent advancements in putting together semi-automatic and fully automatic pipelines for EB data processing. Automatic procedures have already been used to process Hipparcos data, LMC/SMC observations, OGLE and ASAS catalogs etc. We discuss the advantages and shortcomings of these procedures.

With recent observational advancements, substantial amounts of photometric and spectroscopic eclipsing binary data have been acquired. As part of an ongoing effort to assemble a reliable pipeline for fully automatic data analysis, we put Powell's direction set method to the test. The method does not depend on numerical derivatives, only on function evaluations, and as such it cannot diverge. Compared to differential corrections (DC) and Nelder & Mead's downhill simplex (NMS) method, Powell's method proves to be more efficient in terms of solution determination and the required number of iterations. However, its application is still not optimal in terms of time cost. Causes for this deficiency are identified and two steps toward the solution are proposed: non-ortogonality of the parameter set should be removed and better initial directions should be determined before the minimization is initiated. Once these setbacks are worked out, Powell's method will probably replace DC and NMS as the default minimizing algorithm in PHOEBE modeling package.

Eclipsing binary stars are gratifying objects because of their unique geometrical properties upon which all important physical parameters such as masses, radii, temperatures, luminosities and distance may be obtained in absolute scale. This poses strict demand on the model to be free of systematic effects that would influence the results later used for calibrations, catalogs and evolution theory. We present an objective scheme of obtaining individual temperatures of both binary system components by means of color-index constraining, with the only requirement that the observational data-set is acquired in a standard photometric system. We show that for a modest case of two similar main-sequence components the erroneous approach of assuming the temperature of the primary star from the color index yields temperatures which are systematically wrong by ~100K.

PHOEBE (PHysics Of Eclipsing BinariEs) is a modeling package for eclipsing binary stars, built on top of the widely used WD program (Wilson & Devinney 1971). This introductory paper overviews most important scientific extensions (incorporating observational spectra of eclipsing binaries into the solution-seeking process, extracting individual temperatures from observed color indices, main-sequence constraining and proper treatment of the reddening), numerical innovations (suggested improvements to WD's Differential Corrections method, the new Nelder & Mead's downhill Simplex method) and technical aspects (back-end scripter structure, graphical user interface). While PHOEBE retains 100% WD compatibility, its add-ons are a powerful way to enhance WD by encompassing even more physics and solution reliability. The operability of all these extensions is demonstrated on a synthetic main-sequence test binary; applications to real data will be published in follow-up papers. PHOEBE is released under the GNU General Public License, which guarranties it to be free, open to anyone interested to join in on future development.

Eclipsing binaries are extremely attractive objects because absolute physical parameters (masses, luminosities, radii) of both components may be determined from observations. Since most efforts to extract these parameters were based on dedicated observing programs, existing modeling code is based on interactivity. Gaia will make a revolutionary advance in shear number of observed eclipsing binaries and new methods for automatic handling must be introduced and thoroughly tested. This paper focuses on Nelder & Mead's downhill simplex method applied to a synthetically created test binary as it will be observed by Gaia.

Interstellar and atmospheric extinctions redden the observational photometric data and they should be handled rigorously. This paper simulates the effect of reddening for the modest case of two main sequence T1 = 6500K and T2 = 5500K components of a detached eclipsing binary system. It is shown that simply subtracting a constant from measured magnitudes (the approach often used in the field of eclipsing binaries) to account for reddening should be avoided. Simplified treatment of the reddening introduces systematics that reaches \~0.01mag for the simulated case, but can be as high as ~0.2mag for e.g. B8V--K4III systems. With rigorous treatment, it is possible to uniquely determine the color excess value E(B-V) from multi-color photometric light curves of eclipsing binaries.

The orbits and physical parameters of three close, double-lined G0 eclipsing binaries have been derived combining H_P, V_T, B_T photometry from the Hipparcos/Tycho mission with 8480-8740 \AA ground-based spectroscopy. The setup is mimicking the photometric and spectroscopic observations that should be obtained by GAIA. The binaries considered here are all of G0 spectral type, but each with its own complications: V781 Tau is an overcontact system with components of unequal temperature, UV Leo shows occasional surface spots and GK Dra contains a delta-Scuti variable. Such peculiarities will be common among binaries to be discovered by GAIA. We find that the values of masses, radii and temperatures for such stars can be derived with a 1-2% accuracy using the adopted GAIA-like observing mode.

Because of GAIA's estimated harvest of cca. 10^5 eclipsing binaries (Munari et al. 2001) automative procedures for extracting physical parameters from observations must be introduced. We present preliminary results of automative scripting applied to 5 eclipsing binaries, for which photometric and radial velocity observations were taken from literature. Although the results are encouraging, extensive testing on a wider sample has to be performed.

The orbits and physical parameters of three detached, double-lined A-F eclipsing binaries have been derived combining H_P, V_T, B_T photometry from the Hipparcos/Tycho mission with 8500-8750 Ang ground-based spectroscopy, mimicking the photometric+spectroscopic observations that should be obtained by GAIA, the approved Cornerstone 6 mission by ESA. This study has two main objectives, namely (a) to derive reasonable orbits for a number of new eclipsing binaries and (b) to evaluate the expected performances by GAIA on eclipsing binaries and the accuracy achievable on the determination of fundamental stellar parameters like masses and radii. It is shown that a 1% precision in the basic stellar parameters can be achieved by GAIA on well observed detached eclipsing binaries provided that the spectroscopic observations are performed at high enough resolution. Other types of eclipsing binaries (including semi-detached and contact types) and different spectral types will be investigated in following papers along this series.