MNRAS 000 ,

1-37 (2019)

## Modeling Kepler Eclipsing Binaries: Homogeneous Inference of Orbital &amp; Stellar Properties

D. Windemuth, 1 ? E. Agol, 1 2 ; A. Ali, 1 F. Kiefer 3 4 ;

- 1 Department of Astronomy, University of Washington, Seattle WA 98195, USA
- 2 Guggenheim Fellow
- 3 Institut d'Astrophysique de Paris, 75014 Paris, France
- 4 CNES Fellow

Accepted 2019 July 30. Received 2019 July 26; in original form 2018 August 17

## ABSTRACT

We report on the properties of eclipsing binaries from the Kepler mission with a newly developed photometric modeling code, which uses the light curve, spectral energy distribution of each binary, and stellar evolution models to infer stellar masses without the need for radial velocity measurements. We present solutions and posteriors to orbital and stellar parameters for 728 systems, forming the largest homogeneous catalogue of full Kepler binary parameter estimates to date. Using comparisons to published radial velocity measurements, we demonstrate that the inferred properties (e.g., masses) are reliable for well-detached main-sequence binaries, which make up the majority of our sample. The fidelity of our inferred parameters degrades for a subset of systems not well described by input isochrones, such as short-period binaries that have undergone interactions, or binaries with post-main sequence components. Additionally, we identify 35 new systems which show evidence of eclipse timing variations, perhaps from apsidal motion due to binary tides or tertiary companions. We plan to subsequently use these models to search for and constrain the presence of circumbinary planets in Kepler eclipsing binary systems.

Key words: binaries: eclipsing

## 1 INTRODUCTION

A significant fraction of stars live in binaries or higher-order multiple systems (Abt 1979; Duchˆne &amp; Kraus 2013); despite e this, the effects of binarity on a system's formation, stellar evolution, and dynamical evolution, including planetary system evolution, are not well understood. These uncertainties propagate to other fields of astronomy. For example, they have profound implications on interpretations of stellar populations, the foundation of galactic and extragalactic astronomy (Eldridge et al. 2017; Izzard &amp; Halabi 2018). Constraining binary processes requires knowledge of their fundamental parameters, such as radii, temperatures, mass ratio, and age. Fortunately, the geometry of eclipsing binary (EB) stars allows for direct constraints on such parameters. This makes EBs excellent astrophysical tools not only to calibrate single star stellar evolution models (e.g., Stassun et al. 2009; Claret, &amp; Torres 2018; del Burgo &amp; Allende Prieto 2018), but also to test models of binary tidal evolution (Birkby et al. 2012; Lurie et al. 2017; Fleming et al. 2019), along with their planetary systems (Martin et al. 2015; Fleming et al. 2018). However, direct measurements of binary absolute dimensions rely on the fortuitous condition that the binary is an eclipsing, double-lined spectroscopic system, and necessitate multi-band time-series and radial velocity (RV) information, which is observationally expensive to obtain for a large ensemble of systems. Making meaningful inferences about populations of binaries, therefore, may benefit from a different approach.

The goal of this work is to estimate absolute stellar and orbital dimensions for a large ensemble of GLYPH&lt;24&gt; detached eclipsing binaries in the Kepler field, relying only on photometry . Rather than use time-intensive measurements of a small number of bright EBs to constrain stellar models, we leverage theoretical stellar isochrones to infer properties for a much larger sample

? E-mail: windemut@uw.edu

## 2 D. Windemuth et al.

of systems. The basis for our study begins with the Villanova Kepler EB Catalogue (VKEBC; Prˇ sa et al. 2011; Kirk et al. 2016), which contains over 2800 eclipsing and ellipsoidal binaries from the Kepler mission. The majority of these objects have only crude phenomenological estimates of eccentricity, temperature ratios, and sum of radii relative to orbits, without a full characterization of the binary parameters. Detailed, full parameter studies have typically focused on individual systems or smaller samples of binaries for e.g., asteroseismic (Guo et al. 2016; Gaulme et al. 2016), eccentricity-tidal (Kjurkchieva et al. 2016), and RV analysis (Matson et al. 2017). Here, we statistically combine the constraining power of both light curves and SEDs to derive full system solutions and posteriors to 728 Kepler EBs, creating a large homogeneous catalogue.

Such a catalogue has yet to be developed for the binaries in the Kepler field, although comparable analyses have been carried out for other catalogues of eclipsing binaries. The most similar study was applied to the Trans-Atlantic Exoplanet Survey (TrES; a pre-cursor survey to Kepler ) by Devor et al. (2008) using light curves collected to discover transiting exoplanets from the ground, supplemented by multi-band photometry. The advantage of the Kepler field is that it forms a homogeneous dataset with high precision, continuous time-series photometry over a long baseline, and that it is extremely well studied in the characterization of single stars and exoplanets.

The original impetus for our catalogue is to characterize a large population of binaries to a sufficient degree to enable a search for transiting circumbinary planets (CBPs) and to quantify the frequency of CBPs as a function of stellar binary properties (Windemuth et al., in prep). In addition to this, our analysis can provide the basis for additional studies of binary stellar populations. Time-intensive radial velocity follow-up of Kepler EBs would benefit from informative priors on masses, radii, temperatures, and ages of a large EB sample to select the preferred targets to observe. Our homogeneous parameter estimates could also permit more detailed statistical testing of physical processes, such as tidal circularization (Duquennoy &amp; Mayor 1991; Van Eylen et al. 2016) and rotational synchronization (Lurie et al. 2017), and could facilitate detection of additional bodies in the system, such as tertiary stars (e.g., Borkovits et al. 2015, 2016). Our methodology may be applied to current and upcoming surveys in the era of large-scale precision photometry such as OGLE, TESS , and LSST, which have all-sky coverage, opening up an avenue for binary star galactic archaeology.

With these goals in mind, here we present a Bayesian characterization of 728 well- to semi-detached EBs using Kepler time-series photometry and archival multi-band imaging, as well as Gaia parallaxes when available; that is, we derive full system parameters with only photometry and astrometry, and no RV information. Simultaneous modeling of the light curve (LC) and spectral energy distribution (SED) provides joint constraints on EB orbital (e.g., period, eccentricity, argument of periastron, inclination) and stellar (e.g., mass, radius, temperature, age) properties, subject to the assumption that the two stars fall on a single isochrone. The paper is organized as follows. In section 2, we describe our sample selection and the simultaneous LC + SED model. We present the results of our analysis in section 3, and discuss our results in the context of stellar evolution, systematic biases, previous EB radial velocity and population studies in section 4. In section 5, we conclude and summarize our key findings.

## 2 METHODS

We developed a Python code, dubbed 'KEBLAT', to simultaneously fit the light curves and SEDs of 728 Kepler EBs to obtain their orbital and stellar parameters. The open-source code is available on Github . In this section we describe the 1 target selection process ( § 2.1), the datasets gathered ( 2.2), § the light curve modeling ( 2.3), § spectral energy distribution modeling ( 2.4), and joint LC+SED modeling of each binary ( 2.5 and 2.6). To test the accuracy of our model inference, we § § also extract RVs for 3 EBs in our sample that have not been studied previously ( § 2.7); we add these to published EBs with RV mass solutions to compare to our inferred masses below ( 4.1.1). Throughout the paper, we refer to properties of the primary § and secondary components with subscripts 1 and 2, respectively, and define the primary component as the eclipsed star that exhibits the deeper light curve minimum.

## 2.1 Sample selection

We selected our binary systems from the Villanova Kepler EB Catalogue or VKEBC (Prˇ sa et al. 2011; Kirk et al. 2016). This sample, totaling 2877 targets, has been compiled from the entire Kepler prime mission and includes binaries which are not eclipsing (e.g., ellipsoidal variables). We made selection cuts to reduce grazing (i.e., non-constraining) geometries and the effects of close binarity that are not included in our model (e.g., tidal distortions and Doppler boosting of the binary components).

First, to diminish the presence of background, grazing, or non- eclipsing binaries, we required primary and secondary eclipse depths to be &gt; 5% and &gt; 0.1% of the normalized flux, respectively. Next, we excluded targets which exhibit strong distortions in the light curve by selecting for binaries with morphology parameter morph &lt; 0.6. This morphology parameter,

1 https://github.com/savvytruffle/keblat/

provided by VKEBC, is based on a locally linear embedding scheme that classifies light curve shape, such that a low value correlates with well-detached systems exhibiting clearly separated eclipses, and higher morphology value correlates with overcontact (OC) systems with sinusoidal variations (Matijeviˇ et al. 2012). c

In total, we investigated 728 GLYPH&lt;24&gt; detached eclipsing binary systems. The large reduction in the number of systems from the full VKEBC is related to the fact that short-period binaries have the highest eclipse probability, but also have the largest tidal distortions; our sample emphasizes the longer period eclipsing binaries.

## 2.2 Data Acquisition

## 2.2.1 Kepler photometric time series data

We extracted light curves from Kepler long-cadence ( GLYPH&lt;24&gt; 30 min) simple aperture photometry (SAP) based on Data Release 24, spanning GLYPH&lt;24&gt; 4 years or 17 quarters. To correct for slight variations in instrumental sensitivity across quarterly rotations, whereby light from a target may fall upon different detector and pixels, we normalized the raw SAP flux by the quarterly median value. To lessen the impact of outliers in the data, we then de-weighted data points with quality flags &gt; 8 (signifying cosmic ray hits, space craft re-pointing, etc) by inflating their error bars 10-fold. We did not fit the entire Kepler light curve for our final analysis. Instead, we clipped the light curve around each eclipse with a window 1.5-2.0 times eclipse durations, which were initially taken from VKEBC, and then iteratively refined during the optimization process. This window around each eclipse allowed us to fit a polynomial to the variable component of each binary which is not included in our physical model, while minimizing the model evaluation and computational time.

## 2.2.2 Spectral energy distribution data

As we did not have radial velocity data for most of our targets, we required additional constraints to infer the absolute parameters of each system. We collected literature photometric data on each target, and modeled these with stellar evolution models (see § 2.4). We identified extant photometry for each target by cross-matching our sample against sources from the Howell-Everett (HE, Everett et al. 2012), Sloan Digitial Sky Survey (SDSS, Ahn et al. 2012), Two-Micron All-sky Survey (2MASS, Skrutskie et al. 2006), and Wide-field Infrared Survey Explorer (WISE, Cutri &amp; et al. 2014) catalogues on Vizier using a 2 00 search cone. We gathered together the available archival photometry to construct SEDs, excluding targets which had fewer than five photometric data points.

Having multiple measurements spanning a broad wavelength range is crucial to leverage the shape of stellar SEDs to constrain the stellar properties. The U and B flux, in particular, distinguish young and hot stars from cooler, older stars. While the HE and 2MASS campaigns have nearly complete coverage of the Kepler field, the WISE footprint on our sample covers only GLYPH&lt;24&gt; 50%, and the SDSS footprint covers only GLYPH&lt;24&gt; 10%. Therefore, we supplemented SDSS photometry with g 0 r 0 i 0 z 0 estimates from the Kepler Input Catalogue (KIC), which were designed to mimic SDSS bandpasses. The typical agreement between KIC and SDSS griz measurements for our sample is of order GLYPH&lt;24&gt; 0.05 mag. Furthermore, we used the measured ephemerides to identify that about 10-15% of the photometric observations had been taken during an eclipse; we inflated their error bars to 3 mags to diminish their influence during fitting.

## 2.3 Light curve model

The LC module couples a Keplerian orbit solver to the analytic Mandel &amp; Agol (2002) transit model with a quadratic limbdarkening law to fit the observed light curve. The Keplerian orbital model determines the instantaneous relative positions of the binary stars, which are used, in conjunction with radius ratio and a limb darkening prescription, to compute the obscured-to-unobscured system flux. We apply a triangular parameterization of the quadratic limb darkening coefficients to allow uniform prior sampling (Kipping 2013); these coefficients are allowed to float and are not explicitly coupled to model atmospheres.

Using this model, KEBLAT constructs a template light curve with high temporal fidelity at GLYPH&lt;24&gt; 1min, incrementally summed to produce a GLYPH&lt;24&gt; 30 min integration, i.e., the Kepler long-cadence exposure time. The code then uses this high fidelity template to linearly interpolate the flux at each observed cadence. This method speeds up the computation time and accounts for distortions in the light curve due to finite integration times (Kipping 2010). For stellar eclipses, which are much higher signalto-noise than exoplanet transits, this step is especially important; the effect of finite integration time is strongest during ingress and egress, which in turn can affect inference of eccentricity, inclination, radii, and limb darkening coefficients. We test this interpolation scheme against direct computation of the light curve and found that the maximum error is GLYPH&lt;24&gt; 10 GLYPH&lt;0&gt; 6 , at least 2-3 orders of magnitude smaller than typical in-eclipse residuals due to dynamical and stellar variability. We apply quarterly contamination values from Kepler to model (varying) third-light contributions, and allow each quarter's contamination value to float during optimization.

MNRAS 000 , 1-37 (2019)

## 4 D. Windemuth et al.

In total, the light curve model has 13+ N quarters free parameters describing stellar properties ( M sum , R 2 R 1 , R sum , F 2 F 1 ), orbital elements ( e sin ! , e cos ! , b , P, t PE ), limb darkening ( q 11 ; q 12 ; q 21 ; q 22 ), a systematic error term GLYPH&lt;27&gt; sys to account for underestimated observational uncertainties, and crowding values for quarter j ( c j ). Note that we define the impact parameter with respect to the primary, such that b = a cos i GLYPH&lt;157&gt; R 1 . Transformed parameters such as e sin ! , e cos ! , and limb darkening coefficients have finite bounds. For non-transformed variables, we place broad bounds; see Table 1 for a detailed summary of each parameter.

During each model evaluation, we simultaneously fit a second order local polynomial around each eclipse to marginalize the background continuum flux, which may be influenced by instrument variations and/or stellar activity. We used a window of 1.5 - 2 times the primary eclipse duration to ensure our polynomial choice is appropriate, and found that a second order polynomial was typically sufficient. Because the eclipses are masked during polynomial marginalization, higher order polynomials may overfit the data and potentially introduce spurious signals.

## 2.4 Stellar evolution and photometric models

The SED model assumes a co-eval binary specified by age GLYPH&lt;28&gt; , metallicity z , and masses M 1 and M 2 . We utilize PARSEC (Bressan et al. 2012) stellar evolution isochrones, a large grid of stellar models featuring properties such as luminosity, log g , effective temperature, and flux in various bandpasses. We adopt a fine grid spacing, with z grid = [0.001, 0.002, 0.004, 0.008, 0.01, 0.015, 0.02, 0.03, 0.04, 0.05, 0.06], GLYPH&lt;1&gt; GLYPH&lt;28&gt; grid = 0.02, and median GLYPH&lt;1&gt; M grid = 0.006. For a specified set of mass, age, and metallicity, KEBLAT performs a bi-linear interpolation over the PARSEC grid to determine the predicted temperatures and absolute magnitudes in UBV griz JHK W W 1 2 bandpasses, i.e., the stellar SED. We calculate the effective radius of each star from the bolometric luminosity via the Stefan-Boltzmann law; these radii values are necessary to couple to the LC module.

To fit the archival data, we sum the predicted SEDs of both stellar components, compute the distance modulus, and correct for dust extinction along the line of sight assuming an exponential dust distribution with scale-height h 0 = 119 pc (Kruse &amp; Agol 2014):

$$\quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \quad \quad \quad$$

where b G is a target's galactic latitude, d is distance in pc, and E B ' GLYPH&lt;0&gt; V ' and A GLYPH&lt;21&gt; are reddening and band specific extinction computed from a Milky Way extinction law with R = 3 1 : (Fitzpatrick 1999). The integrated absolute magnitude of the binary in a given bandpass is

$$\text{Mag} _ { \lambda, \text{binary} } = \text{Mag} _ { \lambda, 1 } - 2. 5 \log _ { 1 0 } \left ( 1 + 1 0 ^ { - 0. 4 \text{(Mag} _ { \lambda, 2 } - \text{Mag} _ { \lambda, 1 } } \right ) \,.$$

We use cross-matched Gaia distances derived from Bailer-Jones et al. (2018) and Schlafly &amp; Finkbeiner (2011) dust maps results, when available, to place Gaussian priors on d and E B ' GLYPH&lt;0&gt; V ' in our model. Accurate distances from parallax may be used to place better constraints on EB masses since mass correlates tightly with luminosity on the main sequence. Well-calibrated eclipsing binaries can be used as standard candles to calibrate parallaxes (Southworth et al. 2005), and the converse should also be true for systems with well-constrained geometries (Stassun &amp; Torres 2016). In particular, accurate distances may better constrain masses for binaries with non-total or non-annular eclipses by precisely determining total system luminosities. However, rather than using reported uncertainties as fixed GLYPH&lt;27&gt; d , GLYPH&lt;27&gt; E B ' GLYPH&lt;0&gt; V ' , we allow the widths of these Gaussian priors to float to tolerate inaccuracies in the dust map or Gaia data due to source confusion or presence of tertiary companions, which has an occurrence rate of GLYPH&lt;24&gt; 15 GLYPH&lt;0&gt; 20 % in the Kepler field (Gies et al. 2012; Rappaport et al. 2013; Conroy et al. 2014; Orosz 2015). Frequent eclipses or nearby long-period binaries may also deleteriously affect the accuracy of Gaia astrometry.

As with the LC model, the SED model fits for a systematic error, GLYPH&lt;27&gt; SED , to account for underestimated observational errors and uncertainties in the isochrone models. This systematic error also encapsulates any contaminant flux to the binary by additional components, as we do not fit for a third star given the limited number of available SED measurements. As a result, we expect this systematic term to be relatively large compared to reported observational uncertainties.

In practice, we use the sum and ratio of masses ( M sum = M 1 + M 2 , Q = M 2 GLYPH&lt;157&gt; M 1 ), instead of individual masses ( M 1 , M 2 ), to optimize the SED fits. The SED module has in total 10 possible free parameters describing the fundamental stellar properties ( M sum , Q , age, metallicity), dust scaleheight h 0 , GLYPH&lt;27&gt; SED , and distance d and reddening E B ' GLYPH&lt;0&gt; V ' , with their associated uncertainties GLYPH&lt;27&gt; d , GLYPH&lt;27&gt; E B ' GLYPH&lt;0&gt; V ' ; for this work, we fix h 0 = 119 pc. We give descriptions of each parameter and associated bounds in Table 1.

## 2.5 Joint SED and LC model

We couple the isochrone interpolator to the light curve solver for simultaneous SED + LC fitting. During a joint fit, the SED model provides binary radii and Kepler bandpass fluxes, which are used as inputs to the LC module; this effectively reduces the number of free parameters from 13+ N quarters +10 to 20+ N quarters . The full set of parameters is listed in Table 1. We express

the joint likelihood function as the product of likelihoods given all sets of data (e.g., LC, SED, extinction, and distance), assuming that each data point is uncorrelated and uncertainties are Gaussian:

$$\log L & = - \frac { 1 } { 2 } \sum \left ( \frac { \Delta L C _ { i } ^ { 2 } } { ( \sigma _ { L C, i } ^ { 2 } + \sigma _ { L C, s y s } ^ { 2 } ) } + \log ( \sigma _ { L C, i } ^ { 2 } + \sigma _ { L C, s y s } ^ { 2 } ) \right ) - \frac { 1 } { 2 } \sum \left ( \frac { \Delta S E D _ { j } ^ { 2 } } { ( \sigma _ { S E D, j } ^ { 2 } + \sigma _ { S E D, s y s } ^ { 2 } ) } + \log ( \sigma _ { S E D, j } ^ { 2 } + \sigma _ { S E D, s y s } ^$$

Here, GLYPH&lt;1&gt; LC, GLYPH&lt;1&gt; SED, GLYPH&lt;1&gt; d , and GLYPH&lt;1&gt; E(B-V) are the light curve, spectral energy distribution, distance, and reddening fit residuals, respectively. The distance and E(B-V) uncertainties GLYPH&lt;27&gt; d and GLYPH&lt;27&gt; E B ' GLYPH&lt;0&gt; V ' are allowed to float. The systematic error terms for LC and SED are added in quadrature to the observed uncertainties, and importantly the log quantities penalize inflating the error bars. This prevents the model from converging to a poor fit compensated by over-inflating the error bars, and allows the model to weight the LC or SED differently, if the data prefers it. The best-fit solution corresponds to the maximum likelihood model.

## 2.6 Model initialization, optimization, and Bayesian parameter estimation

Due to the non-linear nature of the SED and light curve models, we fit the time-series and multi-band photometry in three stages:

- (i) A fit to only the Kepler light curve.
- (ii) A fit to only the spectral energy distribution.
- (iii) A combined fit to the Kepler light curve and spectral energy distribution, using the prior stages to initialize the fit.

For each stage, we optimize the solutions via lmfit (Newville et al. 2014), a non-linear least squares minimization algorithm based on the Levenberg-Marquardt algorithm (Mor´ 1978). We find that the optimization routine is highly sensitive e to initial parameters; for initial solutions far away from local minimum, the optimization scheme may become stuck, or tend toward boundary values. Consequently, we initialize the optimizer with a combination of gridded parameters and parameter guesses from observables. We describe our fitting procedures in practice below.

For LC fitting, we use observed eclipse depths, durations, and times to estimate initial parameters, and pad these guesses with a grid in radius ratio, e cos ! , e sin ! , and impact parameter. We then optimize each set of initial parameters and save the solution with the smallest chi-square value. For SED fitting, we use a brute-force grid method to find the best-fit SED model(s) with stellar parameters consistent with the light curve solution: we initialize each set of SED fits in a grid of M sum 2 f 0 5 1 0 1 5 2 0 2 5 4 0 : ; : ; : ; : ; : ; : g M GLYPH&lt;12&gt; , Q 2 f 0 1 0 3 0 5 0 9 : ; : ; : ; : g , z value from KIC, E B ' GLYPH&lt;0&gt; V ' value from Schlafly &amp; Finkbeiner (2011) and d from Bailer-Jones et al. (2018). Then, for each set of M 1 ; M 2 ; z , we determine the upper age limit given the initial mass in the PARSEC model, and evenly sample age values in log space. For each nested subgrid, we optimize the fit via lmfit and save solutions with the lowest GLYPH&lt;31&gt; 2 value.

The best-fit parameters from the separate SED and LC solutions seed the joint-fit optimization. When radii values from SED optimization are in tension with LC-derived ones, we pad the separately derived best-fit parameters with a grid of values in mass sum, mass ratio, and age, and step through a grid of these initial parameters for simultaneous SED+LC fitting.

To quantify uncertainties and degeneracies in the model, we use the best-fit solution from the joint SED+LC optimization to seed Monte Carlo Markov Chain (MCMC) simulations. We use an affine-invariant MCMC ensemble sampler, emcee (Foreman-Mackey et al. 2013), with 128 walkers to sample the posterior distribution of parameters for each binary. We place log-uniform priors on the systematic error terms, Gaussian priors on d and E(B-V), and flat priors on all other parameters; see Table 1 for details on parameter bounds. In addition, uniform sampling of age and eccentricity, under parameter transformation to GLYPH&lt;28&gt; = log 10 (age [yr]) and ( e sin !; e cos ! ), requires placing a prior of exp ' GLYPH&lt;28&gt; ln 10 'GLYPH&lt;157&gt; ln 10 and 1 GLYPH&lt;157&gt; e , respectively. As the model already contains a large number of free parameters, we fix the crowding parameters during posterior sampling to the values found from the optimization stage.

We generate Markov chains with 128 walkers for N iter = 100,000 - 700,000 iterations, discarding chains before burn-in for posterior analysis. We define the burn-in period as when the chains cross 5 times the median log-likelihood value of the entire ensemble population, adapting the methodology of Tegmark et al. (2004). Because the chains are initialized from a highdimensional Gaussian distribution centered on the lmfit optimized solution, burn-in is typically low ( . 10 000 ; iterations), and a long burn-in period may signify that the optimized solution was located at a local minimum. For systems exhibiting this behaviour, we restart the Markov chains using the previous run's maximum-likelihood solution.

We use multiple diagnostics to assess convergence. First, we check that the acceptance fraction is between 0.01 and 0.4. Next, we estimate the integrated autocorrelation time, GLYPH&lt;28&gt; acf , of the ensemble and verify that it is appropriate for the chain length. Finally, we visually inspect that the chains appear well-mixed in trace plots. Longer autocorrelation times typically

## 6 D. Windemuth et al.

correlate with smaller acceptance fractions, and large GLYPH&lt;28&gt; acf values may indicate a parameter is not well constrained by the data, the parameter space is multi-modal, or that the observed binary cannot be fully captured by co-eval model isochrones. We note that for shallowly eclipsing systems or those exhibiting tidal distortions, the autocorrelation times for b , age, metallicity, and mass may be long, with effective chain length N iter GLYPH&lt;157&gt; GLYPH&lt;28&gt; acf = 15 - 20, even after &gt; 500,000 generations of 128 walkers. This is because grazing geometries do not well constrain the light curve model, and isochrone models do not well predict stars modified by additional physics. Rather than excluding these EBs from our final sample, for completeness, we report their parameter estimation and discuss at length the factors that limit those results in § 4.1.1 and § 4.2. In contrast, flat-bottomed light curves indicative of total or annular eclipses completely break any degeneracy with inclination. For non-grazing, welldetached systems, our median chain lengths are &gt; 50 times GLYPH&lt;28&gt; acf . We find that period, times of eclipse, and e cos ! converge well for all EBs.

## 2.7 Radial Velocities

In order to validate our method, we compare our inferred stellar masses against RV masses from existing literature studies (see § 4). To supplement the number of RV comparison targets, we mine spectroscopic data from The Apache Point Observatory Galactic Evolution Experiment (APOGEE; Majewski et al. 2017), and extract RVs for three EBs in our sample: KICs 5284133, 5460835, and 6610219. They compose a small overlap of double-lined spectroscopic systems (SB2s) with GLYPH&lt;21&gt; 6 number of APOGEE visits and were not being studied by the EB Working Group.

The radial velocities of the SB2 spectra are derived using a multi-order version of the Two-Dimensional Cross-Correlation algorithm called TODMOR (Zucker &amp; Mazeh 1994; Zucker et al. 2003; Kiefer et al. 2018). We refer the reader to Kiefer et al. (2016) for a full description of the method and summarize the technique in brief here. First, we normalize each APOGEE spectrum by a percentile rank-filtered version of itself (Hodgson et al. 1985; Faigler et al. 2015; Halbwachs et al. 2016). Empirically, we fix the percentile at the 75%-level and the window size at 50 ˚ , a compromise to suppress undesirable deformation of the A pseudo-continuum but conserving the narrower stellar lines. For each target, the spectra with the largest separation of the components are then matched by GLYPH&lt;31&gt; 2 -optimization with 2 synthetic spectra from the PHOENIX library (Hauschildt, Allard, &amp; Baron 1999). The synthetic spectra are also normalized with the same procedure, and broadened to match the line spread function of APOGEE. Averaging over all selected spectra, we obtain stellar parameters for both components of each binary system, including effective temperature, surface gravity, metallicity, flux ratio at 1.6 GLYPH&lt;22&gt; m, and v sin i . Finally, we calculate the two-dimensional cross-correlation function of each observed spectrum and the two synthetic spectra, as explained in Zucker et al. (2003), to derive the RVs of both binary components simultaneously.

## 3 RESULTS

We present the results of our simultaneous SED + LC modeling of 728 Kepler EBs. We find that the majority of our models are good fits to the data. In Tables 2 and 3, we list orbital and stellar solutions for a small subset of EBs in our sample, respectively; these inferred binary parameters correspond to the 50th percentile values of the MCMC chains after burn-in, with error bars corresponding to 16th and 84th percentile values.

The formal uncertainties on mass, age, and metallicity, likely underestimate true, physical uncertainties, as they are reliant on the input stellar isochrone models. Kepler light curves, with their continuous, high-precision 30 min cadence data, offer exquisite constraints on the timing and shapes of eclipses, which map onto the precision of retrieved orbital elements. The inferred geometric elements are particularly precise if the system exhibit flat-bottomed primary and/or secondary eclipse, as total or annular eclipses break degeneracies associated with inclination. For low signal-to-noise or shallow eclipses, however, there is less constraining power in parameters that influence the duration and shapes of the eclipses, such as e sin ! vector, impact parameter (inclination), radius ratio, and flux ratio. This degeneracy is a general problem in eclipse modeling, and may be ameliorated by additional spectroscopic (e.g., flux ratios if systems are SB2s) or astrometric constraints. We note that about 10% of our total EB sample have shallow eclipse depths (PE and SE depths &lt; 0.05), which are susceptible to degeneracies associated with inclination, and thus may not have unique solutions.

Furthermore, the 50th percentile values for each parameter frequently correspond to a model which is a poor fit to the data, because the probability distributions are highly asymmetric; therefore, we also report parameter values corresponding to the maximum likelihood (ML) model in Table 4. We show the ML fits to the SED and light curves, for 5 targets in Figures 1 to 5; these exhibit a range of orbital and stellar parameters, and demonstrate both the versatility and pitfalls of our modeling technique. The full suite of data products for all EBs modeled, including maximum likelihood fits and posterior parameter distributions plots, are available in the Github repository.

## 3.1 Goodness of fit

To quantify the goodness of fit for our EB sample using our technique, we measure median absolute scatter between our models and LC+SED data. Figure 6 shows the distribution of residuals for the time series during eclipse (top) and multiband (bottom) photometric fits, in log space. For a majority (75%) of EBs in our sample, the median residuals are . 1 parts per thousand and . 0.1 mag for light curve and SED, respectively, indicating that in general model solutions well describe the data (see, e.g., Figures 1, 2, and 5). We quantify light curve jitter or 'white' noise by computing the median absolute difference (MAD) in out-of-eclipse flux, and overplot the MAD distribution (red dashed line) against the residual distribution for comparison. Note that the MAD noise level computed here typically underestimates true noise in the Kepler light curves, because it does not account for correlated noise due to e.g., spot modulation, Doppler beaming, and third body perturbation (see e.g., Figure 3). Nevertheless, the majority of systems have model residuals comparable to the white noise, while about 25% of modeled systems appear to have residuals larger than measured jitter. Systems with morphology parameter &gt; 0.5 (shown in orange) systematically exhibit larger residuals in the light curve and in the SED; these are typically short period ellipsoidal or interacting binaries (see, e.g., Figure 4). We further discuss sources of noise that may cause additional scatter during eclipse, as well as sources for the long tail of large residual distributions in § 4.2.

## 3.2 Orbital properties

The observed sky-projected orbits of binaries can be described by their period P , eccentricity e , argument of periastron ! , inclination i , and semi-major axis a . We generate population statistics and histograms using the 50th percentile MCMC posterior values. These distributions provide consistency checks and insight into the physical processes that shaped the binary systems' dynamical states and histories, in the context of observation and selection biases. To fit each distribution, we compute the 50th percentile x-axis values in each bin, and assign Poisson error to bin counts. We adopt this method to generate histograms over convolving posterior probability densities for simplicity, and confirm that this approximation is valid by checking that the 1 GLYPH&lt;27&gt; uncertainties for each system are much smaller than the histogram bin width. Figure 7 shows the resulting distributions for relative orbital size, eccentricity, inclination, and period.

The binary period distribution of our sample (upper left panel of Figure 7) peaks around P GLYPH&lt;24&gt; 4 d and follows a gamma distribution with GLYPH&lt;11&gt; = 7 6 : ; GLYPH&lt;12&gt; = 4 7 : . The shape reflects our target selection process; for reference, the period distribution of the entire VKEBC is shown in dark orange and given half weight to keep the y-axis within range. The decline in the number of P &lt; 3 d EBs in our sample relative to the full catalogue is a direct consequence of selecting for morph &lt; 0.6 binaries, i.e., at least semi-detached (SD) systems, as over-contact systems are more likely found around binaries with small separations; by contrast, the underlying Kepler EB distribution peaks around P GLYPH&lt;24&gt; 0 5 : d.

The inclination distribution is shown in the top right panel, where we fold the values about 90 GLYPH&lt;14&gt; . The distribution peaks around GLYPH&lt;24&gt; edge-on, as expected for eclipsing geometries, and can be roughly fit by an exponential distribution with GLYPH&lt;21&gt; = 0 28 : . The tail of the distribution is somewhat heavy because a large portion of the EBs in our sample have short periods ( P GLYPH&lt;25&gt; 1 GLYPH&lt;0&gt; 2 d), which allows for smaller inclination values (less edge-on) to satisfy the condition for eclipse: j a cos i j GLYPH&lt;20&gt; R 1 + R 2 for circular orbits. We expect this tail may be even larger if we relaxed the morph &lt; 0 6 : criteria, as higher morphology values correlate strongly with shorter periods and thus higher eclipse probability.

The eccentricity distribution (bottom left panel) cannot be described by a single population, as evidenced by the bimodal distribution in log e with peaks around GLYPH&lt;0&gt; 1 6 : and GLYPH&lt;0&gt; 0 5 : ( e GLYPH&lt;25&gt; 0 03 0 3 : ; : , respectively). The excess of GLYPH&lt;24&gt; circular binaries is expected, given the predominance of shorter period P &lt; 10 d binaries in the sample, which are typically tidally circularized; this peak near circular orbit would be even stronger if we included the full VKEBC sample, which contain more short period EBs. The lower right panel of Figure 7 shows the distribution of binary orbital separations as the sum of stellar radii relative to the semi-major axis, R 1 + R 2 a . This is to first order a proxy for light curve morphology; smaller values indicate the stars take up a small fraction of the average orbital separations, while larger values mean the sizes of the stars are comparable to their orbital separation. In the latter case, the stars would experience large tidal forces, resulting in orbital circularization and ellipsoidal deformation of stellar surfaces. Indeed, EBs in our sample with higher morphology values ( morph &gt; 0 5 : ) represent the bulk of systems with large fractional radii. The turnover at larger R 1 + R 2 a values is a consequence of our input criteria to exclude morph &gt; 0.6 binaries, which are more likely to be tidally interacting.

Figure 8 illustrates the morphology (top), inclination (middle), and eccentricity (bottom) of EBs in our sample as a function of their period. Morphology parameter values are truncated at our input cut-off of 0.6, and decline with increasing period, such that P &amp; 10 d systems are typically well-detached binaries. This trend encapsulates the effect of tidal forces which scale inversely and steeply with orbital separation. Shorter period binaries ( P &lt; 5 d) have higher ellipsoidal amplitudes and are more likely to interact, and thus have larger morphology values in general. Roughly 18% of our total sample have morph GLYPH&lt;21&gt; 0.5.

As expected, binary inclination tapers to GLYPH&lt;24&gt; edge-on (90 ) at longer orbital periods (middle panel), due to observational GLYPH&lt;14&gt; bias in the geometry of detecting eclipses. The orbital separations of very short period binaries ( P &lt; 4 d) are small enough that

MNRAS 000 , 1-37 (2019)

## 8 D. Windemuth et al.

grazing eclipses still occur for orbits that are inclined significantly with respect to the observer ( i GLYPH&lt;24&gt; 70 GLYPH&lt;14&gt; ). Eclipsing binaries that harbour circumbinary planets (Doyle et al. 2011; Welsh et al. 2012; Orosz et al. 2012,b; Schwamb et al. 2013; Kostov et al. 2014; Welsh et al. 2014; Kostov et al. 2016) are highlighted as orange stars; they are well-aligned systems with P GLYPH&lt;25&gt; 10 GLYPH&lt;0&gt; 50 d.

The eccentricity-period plot (bottom panel) shows the prevalence of tidal circularization in binaries with orbital periods . 10 d. This result is consistent with findings from previous studies of Kepler binaries (Van Eylen et al. 2016; Price-Whelan &amp; Goodman 2018), although we do not recover a significant difference in the fraction of eccentric binaries with respect to hot and cold stellar components (Van Eylen et al. 2016).

## 3.3 Stellar properties

Next, we present the distributions of stellar parameters ( GLYPH&lt;28&gt; = log 10 (age [yr]), mass ratio Q , surface gravity log g , and temperature T ) of our EB sample in Figure 9. The age distribution (upper left panel) is bimodal, with a strong peak between GLYPH&lt;24&gt; 1-10 Gyr as expected, and another broad but much smaller distribution at younger ages, centered around GLYPH&lt;28&gt; GLYPH&lt;24&gt; 7, or 10 Myr. Were the stellar distribution drawn from a uniform distribution in formation timescale over 10 Gyr, which approximates the star formation of the Galactic disk from which this sample is drawn, then we would only expect GLYPH&lt;25&gt; 1% of stars with ages &lt; 100 Myr and GLYPH&lt;25&gt; 0.1% with ages less than 10 7 yr, rather than the GLYPH&lt;25&gt; 10% for GLYPH&lt;28&gt; &lt; 7 and GLYPH&lt;25&gt; 25% for GLYPH&lt;28&gt; &lt; 8 which comprise our sample. The over-representation of young ages corresponds to instances where single-star isochrones poorly predict the observed stellar properties, and we discuss the age results further in § 4.2.3

The log g distribution (upper right panel) shows that the majority of detached EBs in the Kepler sample are dwarfs ( log g GLYPH&lt;24&gt; 4.5), although about 10% of them have started to evolve off the main sequence as subgiants with log g &lt; 4 . We decompose the log g distributions of primary (dark grey) and secondary (light grey) components, and stack them vertically. The log g distributions for EB primaries and secondaries peak at slightly different values: 4.2 and 4.5, respectively. For reference, we overplot (in dark orange) the distribution of log g from the full Kepler Objects of Interest catalogue (KOI), i.e., single stars. The combined (primary + secondary) log g distribution is similar to the underlying KOI log g distribution in that both KOI and KEBLAT EB surface gravity distributions log g GLYPH&lt;24&gt; 4 5 : , consistent with main sequence dwarfs. However, there is a discrepancy between KOI and overall KEBLAT EB samples in the log g GLYPH&lt;25&gt; 3 8 : GLYPH&lt;0&gt; 4 2 : range, where slightly evolved sub-giants are overrepresented in our sample compared to KOIs. The log g distribution for our sample's primary stars also peaks at a lower log g value than the full KOI sample. These differences may not be significant, given typical uncertainties in KOI log g values are 0.4 (Brown et al. 2011), although it may be related to the higher eclipse probability for stars of larger radius.

The mass ratio distribution ( Q = M 2 GLYPH&lt;157&gt; M 1 ; lower left panel) shows a general positive trend. Typically, Q &gt; 1 systems indicate one or both of the components are not normal stars, i.e., non-MS or mass exchange has modified a system's nascent binary mass ratio. For plotting purposes here, however, we invert the mass ratios for Q &gt; 1 , and stack the resulting histogram on top of the Q GLYPH&lt;20&gt; 1 population, so that the values are in the range [0, 1]. The combined distribution contains a dearth of low-mass companions ( Q &lt; 0 3 : ) and strongly favours similar-mass binaries ( Q GLYPH&lt;24&gt; 0 95 : ), consistent with spectroscopic results of solar-type binaries by Raghavan et al. (2010) and Tokovinin et al. (2006). This increasing Q trend, if physical , should reflect formation values, as our morphology cut disfavours semi-detached binaries and excludes overcontact binaries, e.g., W Uma systems, which are able to exchange mass and evolve away from their nascent Q values (see, e.g., Rucinski 2001; Yakut &amp; Eggleton 2005; Gazeas &amp; Niarchos 2006). We further discuss the reliability of mass ratio inference in § 4.1.1.

Here, we briefly consider how observational biases (e.g., eclipse probability) and the input selection criteria (e.g., morph &lt; 0.6, SE depths) may affect this apparent mass ratio trend; quantifying the combined selection function of Kepler , the VKEBC, and our sample is beyond the scope of this paper. Since eclipse probability and morphology values decrease with increasing period (see Figure 8), we examine the distribution of mass ratios in different period bins. The lower right panel of Figure 9 shows Q distributions (for M 2 GLYPH&lt;157&gt; M 1 GLYPH&lt;20&gt; 1 only) for very short ( P &lt; 4 d; blue solid), short ( 4 GLYPH&lt;20&gt; P &lt; 10 d; orange dashed), and long ( P &gt; 10 d; green dotted) period binaries, where mass bins are widened relative to the cumulative Q distribution to account for low bin counting. Within Poisson counting noise, we find no significant deviations in Q distributions across different binary period bins. We expect eclipse probability to be relatively uniform across Q and not affect their distribution per orbital period bin, except at very low mass ratios which preferentially produce extremely shallow secondary eclipses.

## 3.4 H-R diagram for Kepler EBs

Using predictions from stellar isochrones based on best-fit SED+LC models, we plot luminosity in absolute Kepler magnitudes, temperature, and metallicity for our EB sample in Figure 10. The resulting H-R diagram shows a spread in the main sequence, with a smattering of evolved stars, but no prominent red giant branch. The lack of giants is likely due to a combination of bias in isochrone fitting and target selection. The target selection is affected both by 1) colour cuts in the Kepler catalogue and 2) the morphology cut we made for our sample of GLYPH&lt;24&gt; detached EBs. In the first case, Kepler was designed to detect Earth-size planets around Sun-like stars and selected for FGK dwarfs. In the second case, our morphology cut may bias the sample

because red giants are low density, and thus much more easily distorted by gravity and produce (over-)contact morphologies, which are preferentially excluded. Although there is some giant contamination, the majority of Kepler input catalogue targets are dwarfs. Berger et al. (2018) recently measured precise radii for GLYPH&lt;24&gt; 200,000 Kepler stars using Gaia DR2 parallaxes and properties from the DR25 Kepler Stellar Properties Catalog (Mathur et al. 2017), and found that only 23% and 12% of the stars are subgiants and giants, respectively. We find GLYPH&lt;24&gt; 18 % and GLYPH&lt;24&gt; 10 % of the primaries and secondaries in our sample have log g &lt; 4 (see upper right panel of Figure 9). We discuss how isochrone fitting may bias the age and log g results for KEBLAT EBs in Section 4.

## 3.5 ETV systems and triple candidates

In addition to extracting binary properties, we visually inspected the light curve fits of all EBs modeled and identified 84 systems with eclipse timing variations (ETVs), potentially due to perturbations by a tertiary companion or apsidal motion from tidal deformation of binary components. Out of these 84 systems, one is a confirmed circumbinary planet host (i.e. Kepler-16/KIC 12644769, see Figure 5) and 48 have been previously identified and characterized (Gies et al. 2012; Rappaport et al. 2013; Conroy et al. 2014; Orosz 2015; Borkovits et al. 2016), leaving 35 newly identified candidates. We list all ETV systems identified in this work in Table 5, where additional lettering denotes provenance of previous identification, and reported periods, mass ratios, and eccentricities correspond to maximum-likelihood solutions.

The ETVs typically manifest in the light curve residuals as additional scatter around eclipse ingress and egress, with a 'pinching' at the center of eclipse, see, e.g., Figures 3 and 5. Because ETVs cause eclipses to shift back and forth in time, the residuals take on a shape proportional to the derivative of the light curve. They are largest when the slope of the eclipse is steepest but decline to zero mid-eclipse where the slope is zero, which is responsible for the pinch. Stellar pulsations or starspot modulations may confound ETV signals or lead to false positives; however, these stellar variations tend to affect the entire range of the eclipse and take on a shape similar proportional to the shape of the eclipse. We note that because the systems reported here are inspected by eye, they are typically cases with strong ETV signals and do not comprise a homogeneous sample. Both Conroy et al. (2014) and Orosz (2015) refer to a large, comprehensive study of triples in well-detached Kepler EBs still in development by Orosz, in prep, which would be a source of detailed characterization not carried out here.

## 3.6 APOGEE-derived Mass Ratios

As a part of our investigation into the reliability of KEBLAT (non-RV) mass ratio estimates as compared to RV-derived Q , we measured and modeled RVs from APOGEE spectra for three SB2 systems in our sample, KICs 5284133, 5460835, and 6610219 (see § 2.7 for a description of the RV extraction and model). This is an independent analysis from KEBLAT, without LC or SED information or priors, and relies on software from Kiefer et al. (2018). Absolute masses of the individual components require the inclination constraint from the light curve and cannot be derived from the RV alone; however, the mass ratios depend only on the ratio of RV semi-amplitudes. We present mass ratio solutions for the three SB2s in Table 6, and show the best-fit radial velocity solutions in Figure 11. We discuss these results in the context our larger SED+LC modeling effort in § 4.1.1.

## 4 DISCUSSION

Here, we examine the robustness and limitations of our binary model and results. In particular, we determine the accuracy our EB properties by comparing them to various previous studies published in literature, including RV studies which yield direct measurements on mass (see § 4.1). We are able to reliably infer masses using only photometry for detached, main-sequence binaries. Our solutions are not broadly applicable for EBs with components modified by mass transfer or with red giant companions, as they are not well captured by theoretical isochrones. In § 4.2, we discuss these caveats and devise diagnostics to identify them using observables and model parameters.

## 4.1 Reliability of parameter estimation

## 4.1.1 Mass Measurements

We compile from several radial velocity studies a list of Kepler EBs with published mass solutions common to our sample in Table 7. These masses represent RV 'truths' with which we compare our photometry-only mass inferences. We adopt the mass ratio Q , a quantity independent of inclination, to compare our results; this allows us to use the APOGEE-derived mass ratios from our independent RV analysis.

Figure 12 compares KEBLAT mass ratios to RV solutions from the literature for an overlap sample of 55 EBs. The majority of data points lie close the one-to-one line, with a cluster of discrepant values around KEBLAT Q GLYPH&lt;24&gt; 0.75. We

## 10 D. Windemuth et al.

hypothesize that differences in the sub-population of binaries within the overlap sample may explain the varying levels of Q agreement. We use a mixture model to test our hypothesis and quantify the degree of fitness. Following the convention of Hogg et al. (2010), we include a linear trend (the 'foreground' model) and assume the outliers are drawn from a Gaussian distribution (the 'background' model). We sample the mixture model with emcee and compute the marginalized posterior probability that each data point belongs to the foreground model P FG . We use a threshold P FG = 10 GLYPH&lt;0&gt; 16 to split the sample into 'foreground' (42 EBs) and 'background' (13 EBs) populations for additional analysis.

We investigate the EBs from the 'background' or outlier distribution and find they cluster toward two distinct groups: (i) binaries with morphology values &gt; 0 5 : (e.g., ellipsoidal or Algol-type binaries) and (ii) binaries with one or more red giant components. These systems are not robustly fit by KEBLAT because of the inherent model assumptions; our method does not treat non-Keplerian effects such as tidal deformations that may be present in higher morphology systems, and it is reliant on theoretical isochrones, which cannot self-consistently describe systems that have undergone mass transfer. Moreover, the post-MS evolutionary phases are short lived, where small perturbations in mass or age can lead to large differences, or discontinuities in radius, making it difficult to capture via isochrone fitting. While KEBLAT successfully retrieved the mass ratios of some morph &gt; 0.5 and red giant binaries, it typically failed to reproduce physically self-consistent solutions. Instead, the solutions tended toward young, Q &gt; 1 massive stars with large radius ratios to satisfy the light curve eclipse constraints. One underlying symptom of these mass ratio outliers is that they typically have artificially young derived ages; we discuss this in further detail in § 4.2.3.

We perform linear regression on the 'foreground' EB population (e.g., with outliers removed) in our mixture model and use a systematic uncertainty term GLYPH&lt;27&gt; Q ; sys to quantify the accuracy of KEBLAT-derived masses with respect to RV solutions. We find that they correlate linearly with slope m = 0 95 : GLYPH&lt;6&gt; 0 06 : , intercept b = 0 01 : GLYPH&lt;6&gt; 0 05 : , and systematic error GLYPH&lt;27&gt; Q ; sys = 0 1 : GLYPH&lt;6&gt; 0 02 : . This agreement indicates that KEBLAT mass ratios are reliable to GLYPH&lt;24&gt; 0.1 for well-detached, main-sequence binaries. Additionally, we demonstrate the agreement for absolute masses between KEBLAT and RV analyses. Figure 13 shows primary (filled circle) and secondary (open circle) masses in the overlap sample, after removing EBs with red giant components and morph &gt; 0.5. The inferred masses from photometry-only are in accord with RV values to within 15% of the mass of the star (represented by the grey region). We note that not all of the RV studies included simultaneous light curve modeling in their analysis, which may introduce sin i corrections in the RV absolute masses at the 2% level; this bias is small relative to the scatter between RV and KEBLAT values, and so we do not account for it.

In addition to the RV studies, a small number of Kepler EBs in our sample overlapped with systems observed by the Trans-atlantic Exoplanet Survey (TrES) and modeled by Devor et al. (2008). Their method (MECI) was similar to ours, and used ground based light curves and archival 2MASS photometry in conjunction with isochrones to predict binary masses, but does not use parallax as a constraint upon the model. We applied a quality cut to their catalogue, and selected only for overlap binaries with reduced GLYPH&lt;31&gt; 2 &lt; 2 . Figure 14 shows the agreement in primary and secondary mass between KEBLAT and MECI analysis; in general, the masses agree to 20%, despite both analyses relying only on photometry and adopting different isochrones, data sets, and model details.

For mass inference, the comparison results above indicate that masses derived from photometry here are reliable for well-detached, main-sequence binaries, which comprise the majority of our sample. We recommend approaching systems that have morphology parameter &gt; 0 5 : , young age estimates ( GLYPH&lt;28&gt; . 7 5 : ), and large radius ratio ( R 2 GLYPH&lt;157&gt; R 1 &amp; 2 ; indicative of red giant secondaries) with caution.

## 4.1.2 Orbital elements and stellar radii &amp; temperatures

We compare our inferred eccentricities, inclinations, radii, and temperatures to solutions from the RV studies mentioned above. We supplement this comparison analysis with values from coarse ensemble parameter studies using light curve (Kjurkchieva et al. 2017) and SED (Armstrong et al. 2014) information. Because the literature studies used here have varying degrees of robustness based on their model assumptions and data quality, we only consider their qualitative agreement. In general, KEBLAT solutions agree with published values from small-sample RV studies to within the measurement uncertainties, and share bulk trends with existing larger-scale light curve or SED studies of Kepler EBs.

Figure 15 shows that eccentricity, inclination, and radii values between KEBLAT and literature studies generally lie close to 1:1 relationships, denoted by the dashed black lines. The inclination values show good agreement with RV analyses (see lower left panel), with the exception of Matson et al. (2017) values; we suspect this is because their inclinations are based on machine learning results from phenomenological LC modeling rather than a physical LC model. The eccentricity values (upper left panel) are also generally in very good agreement with published RV work. Our values are broadly consistent with those from Kjurkchieva et al. (2017), although there is large scatter as their approach suffers from using binned phase data and only fitting e and ! while holding other parameters fixed at approximate relation values (see their § 2). That is, they assume the eccentricity vectors do not correlate with radius, temperature, and other orbital parameters. However, binning or down-sampling the LC can change the eclipse profiles, especially at ingress and egress where the shapes are sensitive to

eccentricity vectors, inclination, and stellar radii. Therefore, we expect our full forward modeling on the entire light curve to yield more robust results.

The upper and lower right panels show that our inferred radii demonstrate very good fidelity to results from RV studies. In particular, the red giant systems (Gaulme et al. 2016) show much better agreement in radii than mass (see Figure 12); while the eclipse shape tightly constrains binary radii, the isochrone fitting prefers a more massive MS secondary to a low-mass post-MS component to reproduce its large secondary radius. This may be due to large systematic errors in the isochrones for post-MS evolutionary phases (del Burgo &amp; Allende Prieto 2018) and/or insufficient parameter sampling of complex isochrone morphologies (see § 4.2.3 for further discussion).

Figure 16 compares general Kepler EB properties derived in this study to those reported in Armstrong et al. (2014). The goal of their work was to generate binary temperatures for the entire VKEBC, using Castelli-Kurucz (Castelli, &amp; Kurucz 2003) model atmospheres to fit SEDs constructed from UBVJHK magnitudes. Their analysis differs significantly from ours; in particular, they did not utilize the light curves, except to use eclipse depths reported in the VKEBC to constrain T 2 GLYPH&lt;157&gt; T 1 . As shown in the upper left panel, the effective temperatures of primary components inferred by KEBLAT generally agrees with Armstrong et al. (2014) values as well, although there is a discrepant clump near KEBLAT log T 1 GLYPH&lt;24&gt; 4 1 : vs. Armstrong log T 1 GLYPH&lt;24&gt; 3 8 : . This latter clump represents the eclipsing Algol population in our sample, previously mentioned in § 4.1.1 and discussed further in § 4.2.3. Armstrong et al. (2014) used Gaussian priors centered around KIC values, with an adaptive upper limit of 13000K (see their Table 1), or log T GLYPH&lt;24&gt; 4.1, which explains the dearth of log T 1 &gt; 4 . The temperature ratios generally agree (lower left panel), although there is a large scatter. This is to be expected, since Armstrong et al. (2014) inferred the ratio of temperatures based on relative depth of eclipses as reported by the VKEBC, whereas here, we derive T 2 GLYPH&lt;157&gt; T 1 from a full forward model based on light curve and SED data. The upper right panel shows that primary radii values, normalized by estimated EB distance, agree well, while radius ratio as derived from KEBLAT and Armstrong et al. (2014) has significant scatter (see lower right panel). In particular, the inset overlay shows that for a small range of KEBLAT R 2 GLYPH&lt;157&gt; R 1 values, there is a large, nearly uniform spread (0.5-1.0) in Armstrong radius ratios, which the authors state are poorly constrained in their study (see, e.g., Figures 4 and 5 and section 5.2.2 in Armstrong et al. (2014)). Furthermore, we note that the data points with very large radius ratio values ( R 2 GLYPH&lt;157&gt; R 1 &gt; 2 ) are systems with confirmed red giant components (Gaulme et al. 2016).

## 4.2 Limitations &amp; Diagnostics

We have demonstrated that our models produce good fits to the data for the majority of our sample and accurate parameter inferences for a subset of EBs with RV overlap. Here, we examine our model assumptions, the inherent limitations to our method, and their effects on our results. We recommend using age ( GLYPH&lt;28&gt; . 7 5 : ) and morphology ( morph &gt; 0.5) as diagnostics for systems which may be poorly described by our model, and comment on future improvements.

## 4.2.1 Correlated noise in the light curve

In addition to stochastic 'white' noise, Kepler light curves exhibit correlated noise, due to both instrumental systematics, e.g., telescope drift, and astrophysical effects, e.g., stellar variability, contaminant light from nearby sources (see Gilliland et al. (2015) for a discussion of Kepler noise properties); sources of astrophysical 'noise' may be equivalently considered as 'signals.' A low-order, local polynomial fitting around each eclipse is typically sufficient for light curves of detached systems and minimal stellar variability. For EBs with strong quasi-periodic out of eclipse variation, such as those with morph &gt; 0 5 : which comprise the majority of poor fits in Figure 6, a low-order polynomial may not sufficiently capture the variability, giving rise to larger residuals.

For single star targets, stellar variability, typically in the form of starspot rotation (Giles et al. 2017), stellar oscillation and granulation (Bastien et al. 2013; North et al. 2017), flares (Davenport 2016), confound astrophysical signals. For EBs in our sample, the presence of a secondary star adds two levels of confusion: 1) stellar variability associated with the companion itself and 2) light curve effects due to dynamical interactions and binary evolution, such as ellipsoidal variations, gravity darkening, Doppler beaming, and reflection (Faigler et al. 2013). For case 2, the amplitudes of Doppler beaming and reflection/ellipsoidal effects for similar-mass, P . 5 d solar-type binaries can be of order 1 and 10 ppt, respectively (Zucker et al. 2007). Binaries in our sample with morphology parameter &amp; 0.5 are most affected by these effects, in particular ellipsoidal variations, and account for some of the population with larger residuals in Figure 6 (see orange vs. blue bins).

For the goal of large ensemble binary parameter inference from only photometry using MCMC, including a full physical treatment of tidal distortions and other photodynamical effects would be significantly more computationally expensive. Wellestablished numerical codes that include detailed physics already exist, such as ELC (Orosz, &amp; Hauschildt 2000) and the Wilson, &amp; Devinney (1971) based code PHOEBE (Prˇ sa et al. 2016), where stellar surfaces are modeled on a mesh or grid and allowed to deform. These numerical codes provide higher fidelity to the physics of close binaries at the expense of computational speed. A more computationally tractable method may be to use analytic approximations to explicitly describe ellipsoidal, beaming,

## 12 D. Windemuth et al.

and reflection effects, such as the BEER algorithm (Faigler, &amp; Mazeh 2011). However, this approach would still benefit from an appropriate treatment of additional correlated instrumental noise and stellar variability.

In principle, using Gaussian process (GP) kernels to model instrumental and astrophysical noise in conjunction with a simple physical binary model may lead to better light curve fitting fidelity. The kernel parameters may additionally give insight into the physical processes that underlie the observed quasi-periodic variations. For example, Angus et al. (2018) applied GPs to infer stellar rotation periods for a subset of Kepler objects of interests. In practice, however, we found that including GPs in our method to model Kepler EB light curves significantly increased computational time and required more user input on specific systems to assess kernel type and initial parameter choice; using GPs on K2 or TESS light curves, which typically has GLYPH&lt;24&gt; month long coverage rather than GLYPH&lt;24&gt; year, may be more computationally tractable. We leave implementing GPs on correlated noise as future work, and note that celerite , a newly developed algorithm that scales with the number of data points O' N ' , might make significant improvements to computation time and fitting flexibility.

## 4.2.2 Third Light Contamination

Here, we consider the effects of'third light contamination,'or photometric contamination associated with background/foreground source(s) coincident to or companion(s) gravitationally bound to the binary. As mentioned in § 2.3 and § 2.4, we fit for quarterly Kepler contamination values during model optimization; however, we do not explicitly and self-consistently model flux dilution in the SED by additional sources, due to the limited number of SED measurements, as well as increased model complexity associated with specifying a third isochrone and differential point spread functions (PSFs) with wavelength. Instead, the SED module fits for a single systematic error term, i.e., any additional error beyond observational uncertainties, including third light contamination. This systematic error term should sufficiently capture contributions to the SED from relatively faint contaminators.

Kepler photometry is more susceptible to flux contamination than the Johnson/SDSS/2MASS bands because of the instrument's large pixel scale and aperture radius 2 (e.g., Schwamb et al. 2013; Morton 2012). While WISE bands 1 and 2 share similarly large aperture radii ( GLYPH&lt;24&gt; 6 00 ; Kennedy, &amp; Wyatt 2012), the PSFs from optical and near-IR observations are small. As a result, many of the sources that apppear blended in Kepler and WISE photometry are expected to be resolved and contribute minimally to the bulk of the SED data in the optical and near-IR.

Unresolved contaminators (e.g., point sources at small angular separations &lt; 1 00 ) that are bright relative to the binary (e.g., giants or similar mass dwarfs) will introduce wavelength- dependent bias to the observed SED. If their contribution deviates significantly from Kepler quarterly estimates, they will also dilute the expected depth of eclipses. As we explain below, overall, we expect that only a small fraction (few percent) of our sample are deleteriously affected by bright, unresolved contaminants.

We estimate the extent to which these sources may affect our sample in two ways: 1) by using order of magnitude statistical arguments and 2) by identifying correlated characteristics in inferred model parameters. In the former case, we make the assumption that the majority of these unresolved sources are gravitationally bound, co-eval tertiaries. Tokovinin et al. (2006) found that the fraction of SB2s with tertiary companions in the period range P tertiary = 2 GLYPH&lt;0&gt; 10 5 yr depends strongly on the inner binary period, reaching GLYPH&lt;24&gt; 90% for P . 3 d and GLYPH&lt;24&gt; 35% for P &amp; 10 d. Studies of triples around largely short period ( P &lt; 10 d) Kepler EBs, with P tertiary . 4 yr, have converged on an occurrence rate of GLYPH&lt;24&gt; 15% (e.g., Gies et al. 2012; Rappaport et al. 2013; Conroy et al. 2014). Based on these values, we assume that GLYPH&lt;24&gt; 30% of our EBs host additional unresolved companions. Since the majority of tertiaries have detected masses smaller than the binary mass (Tokovinin 2008; Borkovits et al. 2015), we estimate that 1/3 of triples, or GLYPH&lt;24&gt; 10% of EBs in our sample may harbour companions with significant flux contribution.

In our analysis, we expect that systems with substantial wavelength dependent flux contamination will introduce bias to the total spectral shape (skew) and to the distance modulus (offset). As a result, they may exhibit large inferred systematic SED errors GLYPH&lt;27&gt; sys SED ; and distance errors GLYPH&lt;27&gt; d . We find that 5% of our sample have GLYPH&lt;27&gt; sys SED ; &gt; 0 1 : mag and GLYPH&lt;27&gt; d &gt; 200 pc. This value of a few percent appears consistent with the order of magnitude estimate from the statistical approach. However, as discussed in § 2.4, § 4.2.1, and § 4.2.3, other effects can also cause large GLYPH&lt;27&gt; d and GLYPH&lt;27&gt; sys SED ; values.

As an example, we note that KIC 4862625 is a binary with a pair of companions separated by 0.7 00 in a hierarchical quadruple configuration. The companion binary contributes GLYPH&lt;24&gt; 10% and GLYPH&lt;24&gt; 20% of the total flux in the Kepler and K s bands (Schwamb et al. 2013). Despite this, our stellar inferences ( M 1 = 1 33 : M GLYPH&lt;12&gt; , M 2 = 0 38 : M GLYPH&lt;12&gt; ) agree remarkably well with published values ( M 1 = 1 38 : GLYPH&lt;0&gt; 1 53 : M GLYPH&lt;12&gt; , M 2 = 0 387 : GLYPH&lt;0&gt; 0 408 : M GLYPH&lt;12&gt; ; Schwamb et al. 2013; Kostov et al. 2013).

An inherent limitation to our model is that it relies on theoretical stellar models, which may contain inaccuracies and generate fitting biases. In particular, the complex morphology of isochrones challenges uniform parameter space sampling (Jørgensen &amp; Lindegren 2005), and the predictive power of isochrones is weak for low mass dwarfs (the radius inflation problem; Stassun et al. 2012), post-MS evolutionary phases (e.g., red giants; Cassisi 2017), and stars that exchange(d) mass. These factors lead to discrepant mass values when compared to RV studies ( 4.1.1) and the population of over-represented young ( § GLYPH&lt;28&gt; &lt; 8 ) stars in our general sample ( 3.3). While young ages may be used as consistency checks to the fidelity of stellar isochrones, stringent § tests and calibration of stellar evolution models will require the bulk collection of precise radial-velocity measurements of these eclipsing binaries, such as expanding upon the initial survey of EBs with APOGEE (Zasowski et al. 2013; Badenes et al. 2018).

The morphology of stellar isochrones, in which small perturbations in mass, age, or metallicity can lead to large differences in radius, temperature, and luminosity, challenges uniform sampling of input parameters based on observables. This is particularly true during post-MS evolutionary phases, which apply to RG binaries in our sample. While masses do not always agree, KEBLAT radius values track RV solutions remarkably well for them (see Figure 15). This indicates that the isochrones are converging to young, higher mass MS secondaries with large radius values in order to reconcile R 2 GLYPH&lt;157&gt; R 1 GLYPH&lt;29&gt; 1 based on the eclipses. One means of improving isochrone interpolation for future implementations may be to adopt an 'Equivalent Evolutionary Phases' method (Dotter 2016), in which isochrones are transformed onto a uniform basis that samples evenly different evolutionary phases, rather than evenly in absolute timesteps. Another potential solution for better sampling is to fit the isochrones with a multi-dimensional Gaussian process and compute derivatives of the GP with respect to model parameters; these derivatives would enable more robust optimization and Hamiltonian Monte Carlo analysis to more efficiently sample parameter space.

Another reason why isochrone fitting might fail for red giants is that theoretical radius predictions may be too uncertain (e.g., poor constraints on mass loss, convective overshooting parameters) for the level of precision needed to fit flat-bottomed eclipses in the Kepler light curves. Astereoseismic studies of red giants, which have been revolutionized in the era of Kepler and CoRoT , can provide more reliable and physically motivated constraints on mass relative to our method, and better calibrate stellar models.

In addition to post-MS evolution uncertainty, theoretical models do not reproduce observations of M-dwarfs well. Several studies have shown that low-mass ( M GLYPH&lt;25&gt; 0 1 : GLYPH&lt;0&gt; 0 8 : M GLYPH&lt;12&gt; ) dwarfs have radii inflated relative to theoretical models, up to GLYPH&lt;24&gt; 15 % (Stassun et al. 2012; Kraus et al. 2011; Cruz et al. 2018). This radius inflation is coupled by a corresponding decrease in effective temperature, such that the mass-luminosity relationship (for MS stars) is preserved. The presence of radius discrepancies appears to be agnostic to stellar multiplicity (Kesseli et al. 2018). Multiple studies agree that radius inflation is associated with magnetic field activity, although they diverge in their conclusions about the relative importance of rotationinduced surface magnetic fields and their potential effects on stellar radii via convective inhibition, dark spot coverage, and tidal synchronization in binaries (Chabrier et al. 2007; Stassun et al. 2012; Feiden &amp; Chaboyer 2014; Kesseli et al. 2018).

The PARSEC isochrones we used in our SED+LC model included improved calibration for low mass stars with radius discrepancy at the GLYPH&lt;24&gt; 5 % level (Chen et al. 2014). Binaries with anomalously young M-dwarfs may be explained by the model attempting to compensate for a larger radius than theoretical model predicts; stars that have not yet started core H fusion, i.e., reached the main sequence, are still contracting, and thus have larger radii compared to their MS counterparts. For M dwarfs, the pre-MS stage may last GLYPH&lt;24&gt; 10 Myr to GLYPH&lt;24&gt; 1 Gyr, consistent with inferred youthful ages. Allowing a factor for radius inflation might help in modeling these stars. However, because the luminosity-mass relation for M dwarfs is well calibrated, a radius inflation factor necessarily requires a corresponding scaling relation to lower the effective temperature, which in turn affects the predicted SED. This makes treating the issue in a self-consistent way difficult.

Finally, our model assumes that stars fall on coeval isochrones, an assumption which can be strongly violated in binaries which have undergone mass-transfer. The majority of binaries with anomalously young ages in our general sample and with discrepant masses in our RV comparison sample appear to be Algol-type binaries. These are systems with primary masses between 2-4 M GLYPH&lt;12&gt; (Ibanoˇlu g et al. 2006), V-shaped eclipses, and typically exhibiting semi-detached morphologies (see e.g., Figure 4). In these eclipsing Algol systems, mass transfer occurs from the more evolved and initially more massive component which has filled its Roche lobe onto the less massive companion; as the companion accretes more mass, it is 'rejuvenated' becoming hotter, more massive, and delaying evolution away from main sequence. Stellar evolutionary tracks of single-stars cannot accurately capture masses and ages of stars in Algol-type binaries because they do not take into account modification due to binary evolution. Because the primary components of Algol-type binaries are typically 'main sequence' stars with high mass (O, B, A dwarfs), we expect the ages of this population to be skewed toward young populations. One possible ad-hoc solution to these rejuvenated systems is to allow the ages of both stars to vary independently, although such a treatment may not be physically self-consistent, and the additional degree of freedom may pose further difficulties in efficiently sampling the full parameter space. Alternatively, the same mechanism that enables mass transfer, i.e., when a star fills its Roche lobe,

## 14 D. Windemuth et al.

can be used to constrain mass ratios in EBs with semi-detached or over-contact light curve morphologies; this method has demonstrated success for totaling eclipsing SD or OC systems (Wyithe &amp; Wilson 2002).

## 5 CONCLUSIONS

We have modeled the light curves and SEDs of 728 GLYPH&lt;24&gt; detached Kepler eclipsing binaries, using Bayesian inference for full binary parameter estimation. Our forward model couples Padova stellar isochrones to a Mandel &amp; Agol (2002) transit model to solve the system's stellar and orbital parameters with only photometric information. The entire set of posteriors and maximum likelihood solutions are available online.

We are interested in obtaining reliable orbital and mass solutions, in particular as the basis of an automated and optimized search for circumbinary planets. We find that binary orbital parameters are typically very well constrained by the light curve data, and comparisons for a subset of our sample with previous studies indicate that the majority are reliable.

Moreover, we show that detached, main sequence binaries give accurate masses relative to RV 'truths.' Binary systems which have been affected by stellar evolution, or contain stars which are otherwise discrepant with isochrone predictions, however, may give poor agreement. The latter exhibit symptoms of higher morphology values ( morph &amp; 0.5) and artificially young stellar ages ( GLYPH&lt;28&gt; . 7 5 : ), which may be used as a filter for quality check. These symptoms are manifestations of the inherent limitations associated with our approach: namely, that we do not account for non-Keplerian effects like ellipsoidal variations, and that we rely on co-eval isochrone models; the inferred stellar parameters are only as accurate as input theoretical models. On the other hand, discrepant results that arise from these limitations may help inform which binary targets warrant more time-intensive follow-up observations, to better calibrate theoretical stellar models. For future work, our methodology may be improved by additional Gaia systemic radial velocities. Currently Gaia does not provide these for targets flagged as binaries; once Gaia releases its spectroscopic catalogue, systemic RV can be estimated and, in conjunction with proper motion, the full 3-dimensional motion can provide additional constraints to binary ages based on galactic kinematics. Additionally, as we have discussed in § 4.2, incorporating Gaussian processes into the time series and isochrone modeling may improve the robustness of our model to non-Keplerian light curve effects and post-MS stars.

We conclude by summarizing our key takeaways below:

- (i) We have created the largest, homogeneous catalogue of Kepler eclipsing binaries with full system parameters and posteriors from Bayesian inference. This catalogue is of interest to identify compelling targets for RV follow-up (e.g., M-dwarf binaries), and to binary population studies which require mass estimates and precise measurements of orbital elements (e.g., eccentricities), such as tidal theory and circumbinary planet searches.
- (ii) We additionally identify 35 new systems with eclipse timing variations, perhaps arising from apsidal motion due to binary tides and/or tertiary companions.
- (iii) We demonstrate that we can reliably retrieve mass using only photometry for well-detached, main sequence binaries, using comparisons to RV 'truths.' These systems make up the majority of our sample.
- (iv) For a subset of binaries with post-main sequence components or significant tidal deformations (e.g., semi-detached systems where mass exchange have occurred), our inferred masses are not universally reliable, as our model does not account for non-spherical distortions and are limited by the fidelity of input stellar isochrones. We recommend using GLYPH&lt;28&gt; . 7 5 : and morph &gt; 0.5 to diagnose potentially problematic systems. Fortunately, photometric masses may be derived using asteroseismology (Gaulme et al. 2016) for red giants and by exploiting Roche lobe-filling configurations for semi-detached and over-contact binaries (Wyithe &amp; Wilson 2002).
- (v) The technique used here can be applied to OGLE, TESS, and other large time-domain photometric surveys coming online, such as LSST, to characterize large numbers of eclipsing binaries. As many of these surveys have all-sky coverage, this poses an intriguing opportunity to probe binary galactic archeaology.

## 6 ACKNOWLEDGEMENTS

This work was supported by NASA grants NNX12AF20G, NNX13AF62G, NNX15AT44H (NESSF), and NSF grant 1615315. This work was facilitated through the use of advanced computational, storage, and networking infrastructure provided by the Hyak supercomputer system at the University of Washington, and made use of the gaia-kepler.fun crossmatch database created by Megan Bedell. We thank the anonymous referee for their helpful inputs and constructive criticism. DW would like to thank Dr. Leslie Hebb and Dr. James Davenport for useful discussions on binaries.

Table 1. Model Parameters

| Parameter                    | Meaning                                                        | Module y   | Constrained by                 | Priors                                                            |
|------------------------------|----------------------------------------------------------------|------------|--------------------------------|-------------------------------------------------------------------|
| P                            | period (d)                                                     | LC         | time of eclipses               | flat (0, 4000)                                                    |
| e sin !                      | transformation of eccentricity                                 | LC         | ratio of PE + SE durations     | flat (-0.95, 0.95)                                                |
| e cos !                      | & longitude of periastron                                      | LC         | time of PE + SE                | flat (-0.95, 0.95)                                                |
| t PE                         | time of primary eclipse (BJD-2454833)                          | LC         | LC minimum                     | flat (0, 4000)                                                    |
| b                            | impact parameter a cos i R                                     | LC         | shape of eclipse               | flat (-10, 10)                                                    |
| q 11 ; q 12                  | 1 transformed quadratic limb darkening                         | LC         | shape of eclipse               | flat [0, 1]                                                       |
| q 21 ; q 22                  | coeffs for primary & secondary                                 | LC         | shape of eclipse               | flat [0, 1]                                                       |
| F rat                        | Kepler flux ratio F 2 F 1                                      | LC only    | relative depth of eclipse      | flat [1e-8, 100]                                                  |
| R sum                        | sum of radii ( R GLYPH<12> )                                   | LC only    | duration of eclipse            | flat [0.1, 1e6]                                                   |
| R rat                        | ratio of radii R 2 R                                           | LC only    | eclipse shape + relative depth | flat [1e-6, 1000]                                                 |
| M sum                        | 1 sum of masses ( M GLYPH<12> )                                | SED, LC*   | shape of SED                   | flat [0.2, 24]                                                    |
| Q                            | ratio of masses M 2 M 1                                        | SED        | shape of SED                   | flat [0.0085, 2]                                                  |
| z                            | binary metallicity, defined as 1-X-Y ( z GLYPH<12> = 0 : 015 ) | SED        | shape of SED                   | flat [0.001, 0.06]                                                |
| GLYPH<28>                    | binary log 10 age ( log 10 (yr))                               | SED        | shape of SED                   | flat [6, 10.1]                                                    |
| d                            | distance to system (pc)                                        | SED        | SED offset                     | Gaussian z d Gaia ; GLYPH<27> d                                   |
| E ' B GLYPH<0> V '           | reddening, assuming R V = 3 : 1                                | SED        | shape of SED                   | Gaussian E ' B GLYPH<0> V ' Schlafly ; GLYPH<27> E ' B GLYPH<0> V |
| h 0                          | dust vertical scaleheight (pc)                                 | SED        | shape of SED                   | fixed at 119                                                      |
| GLYPH<27> d                  | uncertainty in Gaia distance prior ( ln )                      | SED        |                                | flat [-1, 7]                                                      |
| GLYPH<27> E ' B GLYPH<0> V ' | uncertainty in reddening prior ( ln )                          | SED        |                                | flat [-12, 2]                                                     |
| GLYPH<27> sys ; SED          | systematic SED error ( ln )                                    | SED        |                                | flat [-25, 0]                                                     |
| GLYPH<27> sys ; LC           | systematic LC error ( ln )                                     | LC         |                                | flat [-25, -4]                                                    |
| c j                          | crowding parameter for quarter j                               | LC         | depth of eclipses              | flat [0, 1]                                                       |

Note. y All parameters belong to the simultaneous SED+LC module, unless otherwise specified, e.g., 'LC only.'

*When only light curve data available, M sum is a proxy for a , to obtain the sizes of stars relative to their separation.

z When no Gaia parallax available, the distance prior is flat on [10, 15000] pc.

Table 2. Eclipsing Binary Orbital &amp; Distance Parameter Posteriors

| KIC                                          | P (d)                                                                                                                                                                                              | t PE (BJD-2454833)   | t PE (BJD-2454833)                                                                                                                                                                                        | e sin !                           | e sin !                                                                                   | e cos !                                | e cos !                                                                                             | i (rad)                             | i (rad)                                                                                     | d (pc)                                                                                                |    |
|----------------------------------------------|----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|----------------------|-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|-----------------------------------|-------------------------------------------------------------------------------------------|----------------------------------------|-----------------------------------------------------------------------------------------------------|-------------------------------------|---------------------------------------------------------------------------------------------|-------------------------------------------------------------------------------------------------------|----|
| 10031409 10198109 12356914 10619109 12644769 | 4.143879363 + 4 : 5 e GLYPH<0> 4 : 6 e 17.91874241 + 3 : 8 e GLYPH<0> 3 : 9 e 27.30822669 + 9 : 3 e GLYPH<0> 9 : 4 e 2.045166167 + 3 : 4 e GLYPH<0> 4 : 7 e 41.07759097 + 1 : 4 e GLYPH<0> 1 : 4 e |                      | 124.5491837 + 9 : 6 e GLYPH<0> 1 : 0 e 146.3930259 + 1 : 8 e GLYPH<0> 1 : 8 e 143.5001244 + 2 : 7 e GLYPH<0> 2 : 8 e 122.139347 + 1 : 7 e - 5 GLYPH<0> 1 : 6 e - 5 132.6583057 + 2 : 8 e GLYPH<0> 2 : 9 e | 0.0047 0.2787 0.4245 0.0676 .1641 | + 0 : 0014 : 0012 + 0 : 0025 : 0024 + 0 : 0154 : 0153 + 0 : 0043 : 0043 + 0 : 0005 : 0005 | 0.000007 0.2484 0.1213 0.00136 .017971 | + 0 : 000004 : 000004 + 0 : 0002 : 0002 + 0 : 0009 : 0010 + 0 : 00008 : 00013 + 0 : 000004 : 000004 | 1.4923 1.5825 1.5642 1.2260 1.57659 | + 0 : 0002 : 0002 + 0 : 0003 : 0003 + 0 : 0014 : 0013 + 0 : 0027 : 0027 + 0 : 00002 : 00001 | 941 + 23 GLYPH<0> 22 637 + 48 GLYPH<0> 22 1344 + 1 GLYPH<0> 14 1081 + 24 GLYPH<0> 22 68.6 + 0 : 5 : 5 |    |

Note. - Posteriors for orbital &amp; distance parameters for a subset of our sample which corresponds to Figures 1 to 5, to demonstrate the table's form and content. Full posteriors for all parameters available in online table. See Table 1 for details on parameter meanings.

Table 3. Eclipsing Binary Stellar Parameter Posteriors

| KIC                                          | z                                                                                                                                                                            | GLYPH<28> (log(yr))                                                                                                                                          | M 1 ' M GLYPH<12> '                                                                                                                                         | M 2 ' M GLYPH<12> '                                                                                                                                         | R 1 ' R GLYPH<12> '                                                                                                                                                     | R 2 ' R GLYPH<12> '                                                                                                                                                     | F 2 GLYPH<157> F 1                                                                                                                                                                        |
|----------------------------------------------|------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|--------------------------------------------------------------------------------------------------------------------------------------------------------------|-------------------------------------------------------------------------------------------------------------------------------------------------------------|-------------------------------------------------------------------------------------------------------------------------------------------------------------|-------------------------------------------------------------------------------------------------------------------------------------------------------------------------|-------------------------------------------------------------------------------------------------------------------------------------------------------------------------|-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| 10031409 10198109 12356914 10619109 12644769 | 0.011 + 0 : 002 GLYPH<0> 0 : 002 0.006 + 0 : 001 GLYPH<0> 0 : 002 0.007 + 0 : 002 GLYPH<0> 0 : 001 0.0023 + 0 : 0003 GLYPH<0> 0 : 0002 0.0100 + 0 : 00005 GLYPH<0> 0 : 00004 | 9.37 + 0 : 16 GLYPH<0> 0 : 20 9.57 + 0 : 28 GLYPH<0> 0 : 71 9.77 + 0 : 20 GLYPH<0> 0 : 33 6.25 + 0 : 01 GLYPH<0> 0 : 01 10.0796 + 0 : 0003 GLYPH<0> 0 : 0006 | 1.13 + 0 : 04 GLYPH<0> 0 : 05 1.01 + 0 : 15 GLYPH<0> 0 : 08 0.94 + 0 : 08 GLYPH<0> 0 : 06 2.48 + 0 : 04 GLYPH<0> 0 : 05 0.6002 + 0 : 0005 GLYPH<0> 0 : 0003 | 1.08 + 0 : 04 GLYPH<0> 0 : 05 0.36 + 0 : 03 GLYPH<0> 0 : 01 0.60 + 0 : 02 GLYPH<0> 0 : 02 1.87 + 0 : 02 GLYPH<0> 0 : 02 0.1913 + 0 : 0003 GLYPH<0> 0 : 0002 | 1.183 + 0 : 009 GLYPH<0> 0 : 009 1.078 + 0 : 045 GLYPH<0> 0 : 025 1.026 + 0 : 014 GLYPH<0> 0 : 018 2.619 + 0 : 016 GLYPH<0> 0 : 020 0.6092 + 0 : 0004 GLYPH<0> 0 : 0002 | 1.085 + 0 : 022 GLYPH<0> 0 : 028 0.347 + 0 : 014 GLYPH<0> 0 : 008 0.587 + 0 : 014 GLYPH<0> 0 : 015 3.038 + 0 : 018 GLYPH<0> 0 : 018 0.2161 + 0 : 0003 GLYPH<0> 0 : 0002 | 0.7689 + 0 : 0245 GLYPH<0> 0 : 0289 0.00473 + 0 : 00002 GLYPH<0> 0 : 00002 0.0602 + 0 : 0020 GLYPH<0> 0 : 0013 0.2093 + 0 : 0131 GLYPH<0> 0 : 0118 0.01582 + 0 : 00003 GLYPH<0> 0 : 00003 |

Note. - Posteriors for stellar parameters for a subset of our sample which corresponds to Figures 1 to 5, to demonstrate the table's form and content. Full posteriors for all parameters available in online table. See Table 1 for details on parameter meanings.

Table 4. Maximum Likelihood Parameter Solutions

| Parameter                       |   KIC 10031409 |   KIC 10198109 |   KIC 12356914 |   KIC 10619109 |   KIC 12644769 |
|---------------------------------|----------------|----------------|----------------|----------------|----------------|
| M sum                           |       2.2751   |       1.2485   |       1.5222   |       4.359    |       0.7919   |
| Q                               |       0.9536   |       0.3731   |       0.6604   |       0.752    |       0.3188   |
| z                               |       0.0138   |       0.0079   |       0.0087   |       0.0023   |       0.01     |
| GLYPH<28>                       |       9.2752   |       9.9106   |       9.8953   |       6.2404   |      10.08     |
| d                               |     963.48     |     608.77     |    1343.87     |    1073.79     |      68.63     |
| E(B-V)                          |       0.2117   |       0.2395   |       0.2214   |       0.5907   |       0.0664   |
| h 0                             |     119        |     119        |     119        |     119        |     119        |
| P                               |       4.14388  |      17.9187   |      27.3082   |       2.04517  |      41.0776   |
| t PE                            |     124.549    |     146.393    |     143.5      |     122.139    |     132.658    |
| e sin !                         |       0.003968 |       0.276679 |       0.443474 |       0.066164 |      -0.163677 |
| e cos !                         |       6e-06    |       0.24854  |       0.120055 |       0.001243 |      -0.017973 |
| b                               |       0.94     |      -0.3405   |       0.3313   |       1.4283   |      -0.441    |
| q 11                            |       0.1732   |       0.2798   |       0.4472   |       0.9994   |       0.1763   |
| q 12                            |       0.3858   |       0.4381   |       0.8051   |       0.0002   |       0.999    |
| q 21                            |       0.3624   |       0.603    |       0.9237   |       0.599    |       0.2659   |
| q 22                            |       0.1984   |       0.9802   |       0.8589   |       0.319    |       0.8639   |
| ln GLYPH<27> sys ; LC           |      -7.2      |      -8.1      |      -4.3      |      -6.5      |      -8        |
| ln GLYPH<27> sys ; SED          |      -3.5      |      -3.7      |      -3        |      -2.4      |      -3.4      |
| ln GLYPH<27> E ' B GLYPH<0> V ' |     -11.1      |      -2        |      -2        |      -0.7      |      -6.4      |
| ln GLYPH<27> d                  |       0        |       6.7      |      -6.1      |       0        |       2.4      |
| c 0                             |       0.9917   |       0        |       0        |       0.9885   |       0.9847   |
| c 1                             |       0.9909   |       0.9927   |       0.9339   |       0.9885   |       0.9784   |
| c 2                             |       0.9898   |       0.9909   |       0.9846   |       0.9901   |       0.9859   |
| c 3                             |       0.9864   |       0.991    |       0.993    |       0.9897   |       0.9902   |
| c 4                             |       0.9867   |       0.9902   |       0.9793   |       0.9886   |       0.9902   |
| c 5                             |       0.991    |       0.9927   |       0.9362   |       0.9884   |       0.9849   |
| c 6                             |       0.9927   |       0.9925   |       0.9884   |       0.9897   |       0.9874   |
| c 7                             |       0        |       0.991    |       0.9932   |       0        |       0.9929   |
| c 8                             |       0.9888   |       0.9902   |       0.9617   |       0.9885   |       0.9849   |
| c 9                             |       0.9887   |       0.9927   |       0.887    |       0.9884   |       0.9842   |
| c 10                            |       0.9925   |       0.9925   |       0.9718   |       0.9897   |       0.9928   |
| c 11                            |       0        |       0.9909   |       0.9667   |       0        |       0.9888   |
| c 12                            |       0.9887   |       0.9892   |       0.9872   |       0.9885   |       0.9922   |
| c 13                            |       0.9899   |       0.9927   |       0.9442   |       0.9885   |       0.9814   |
| c 14                            |       0.9851   |       0.9925   |       0.9958   |       0.9897   |       0.9955   |
| c 15                            |       0        |       0.9909   |       1        |       0        |       0.9859   |
| c 16                            |       0.9877   |       0.9892   |       0.9902   |       0.9885   |       0.9949   |
| c 17                            |       0.9904   |       0.9927   |       0.9388   |       0.9885   |       0.9818   |
| morph y                         |       0.35     |       0.07     |       0.03     |       0.55     |       0.03     |

Note. - Maximum likelihood solutions for a subset of our sample which corresponds to Figures 1 to 5, to demonstrate the table's form and content. Full version available online. See Table 1 for details on parameter meanings and units. Values of 0.0 denotes no data was available. y morph is the morphology value taken directly from the Villanova Kepler Eclipsing Binary Catalogue (Prˇ sa et al. 2011; Kirk et al. 2016).

Table 5. ETV candidates identified in our sample

| KIC             | P (d)                  | Q         | e               | Other Sources y   |
|-----------------|------------------------|-----------|-----------------|-------------------|
| 10095512        | 6.01720794             | 0.83      | 2.0e-03         | brc               |
| 10215422        | 24.84708632            | 0.64      | 2.9e-01         |                   |
| 10268903        | 1.10397904             | 0.82      | 4.5e-03         | b                 |
| 10292238        | 143.11911400           | 0.86      | 6.2e-01         |                   |
| 10296163        | 9.29674544             | 0.52      | 3.8e-01         | b                 |
| 10352603        | 32.77898645            | 1.01      | 4.7e-01         |                   |
| 10549576        | 9.08935880             | 0.52      | 6.9e-02         | b                 |
| 10583181        | 2.69635500             | 0.58      | 3.2e-07         | b                 |
| 10619109        | 2.04516617             | 0.75      | 6.6e-02         |                   |
| 10686876        | 2.61841459             | 0.46      | 1.0e-03         | b                 |
| 10736223        | 1.10509420             | 0.72      | 1.2e-02         |                   |
| 10909274        | 39.23813527            | 1.00      | 7.3e-01         |                   |
| 10979716        | 10.68409498            | 0.80      | 1.5e-01         | b                 |
| 11071207        | 8.04963317             | 1.01      | 3.1e-01         |                   |
| 11234677        | 1.58742711             | 0.45      | 1.6e-01         | b                 |
| 11499757        | 12.31440178            | 0.89      | 2.6e-01         |                   |
| 11558882        | 73.91782410            | 0.60      | 4.2e-01         | b                 |
| 11724091        | 1.55909021             | 0.50      | 1.9e-02         |                   |
| 11923819        | 33.15943126            | 1.02      | 3.3e-01         |                   |
| 12356914        | 27.30822721            | 0.66      | 4.6e-01         | b                 |
| 12459731        | 14.18110967            | 0.77      | 4.1e-02         | b                 |
| 12644769        | 41.07759104            | 0.32      | 1.6e-01         |                   |
| 1995732 2305372 | 77.36197359 1.40469172 | 1.09 0.77 | 1.3e-01 2.8e-02 | b                 |
| 2306740         | 10.30698967            | 0.86      | 3.1e-01         |                   |
| 2576692         | 87.87820528            | 1.00      | 2.1e-01         | b                 |
| 3241619         | 1.70334728             | 0.76      | 6.5e-03         |                   |
| 3247294         | 67.42012708            | 1.33      | 6.0e-01         | o                 |
| 3440230         | 2.88110019             | 0.72      | 1.0e-01         | b                 |
| 3757778         | 36.51436879            | 0.43      | 3.6e-01         |                   |
| 4077442         | 0.69284258             | 0.74      | 1.2e-01         |                   |
| 4544587         | 2.18911101             | 0.88      | 2.8e-01         |                   |
| 4753988         | 7.30445219             | 0.71      | 1.1e-02         | b                 |
| 4948863         | 8.64355067             | 0.68      | 4.8e-03         | b                 |
| 5039441         | 2.15138428             | 0.52      | 2.8e-03         | brc               |
| 5113053         | 3.18509078             | 0.87      | 6.0e-05         |                   |
| 5269407         | 0.95887119             | 0.53      | 2.3e-04         | b                 |
| 5288543         | 3.45707832             | 0.80      | 6.4e-03         |                   |
| 5513861         | 1.51021130             | 0.95      | 2.2e-06         | b                 |
| 5553624         | 25.76208222            | 0.83      | 5.6e-01         |                   |
| 5632781         | 11.02520265            | 1.01      | 2.7e-01         |                   |
| 5731312         | 7.94636806             | 0.77      | 4.7e-01         | b                 |
| 5955321         | 11.63787579            | 0.98      | 4.8e-01         |                   |
| 5962716         | 1.80459191             | 0.71      | 2.8e-01         | b                 |
| 6029130         | 12.59165677            | 0.94      | 1.5e-02         |                   |
| 6042116         | 5.40715640             | 0.64      | 9.7e-02         |                   |
| 6449358         | 5.77679432             | 0.43      | 4.0e-04         |                   |
| 6464285         | 0.84365137             | 0.61      | 1.6e-02         |                   |
| 6525196         | 3.42059775             | 0.94      | 7.1e-05         | brc               |
| 6543674         | 2.39103076             | 0.97      | 4.5e-02         | b                 |
| 6545018         | 3.99145640             | 0.83      | 2.7e-03         | brc               |
| 6610219         | 11.30099291            | 1.01      | 2.1e-01         |                   |
| 6877673         | 36.75871004            | 0.81      | 2.0e-01         | b                 |
| 7021177         | 18.64532032            | 0.97      | 5.9e-01         |                   |
| 7025540         | 2.14821893             | 0.96      | 1.2e-02         |                   |
| 7137798         | 2.25353766             | 0.46      | 1.1e-01         |                   |

Table 5 (cont'd)

|     KIC |    P (d) |    Q |       e | Other Sources y   |
|---------|----------|------|---------|-------------------|
| 7177553 | 17.9965  | 1.45 | 0.42    | b                 |
| 7257373 | 10.4669  | 0.99 | 0.00096 |                   |
| 7385478 |  1.65547 | 0.48 | 0.12    | b                 |
| 7630658 |  2.15116 | 0.94 | 0.00022 | b                 |
| 7670617 | 24.7037  | 0.49 | 0.18    | b                 |
| 7812175 | 17.7941  | 0.49 | 0.14    | b                 |
| 7821010 | 24.2382  | 0.96 | 0.68    | b                 |
| 8143170 | 28.7875  | 0.68 | 0.12    | b                 |
| 8210721 | 22.6732  | 0.36 | 0.13    | b                 |
| 8411947 |  1.79768 | 0.86 | 0.018   |                   |
| 8429450 |  2.70515 | 0.87 | 0.0097  | b                 |
| 8444552 |  1.1781  | 0.79 | 0.15    | b                 |
| 8553788 |  1.60616 | 0.7  | 0.16    | b                 |
| 8553907 | 42.0322  | 1.02 | 0.51    | o                 |
| 8701327 |  8.50884 | 0.91 | 0.0053  |                   |
| 9007918 |  1.38721 | 0.48 | 0.00028 | b                 |
| 9053086 |  1.27484 | 0.55 | 0.15    |                   |
| 9110346 |  1.79055 | 0.93 | 9.4e-05 | b                 |
| 9392702 |  3.90931 | 0.58 | 0.016   | b                 |
| 9451096 |  1.25039 | 0.64 | 0.024   | brc               |
| 9637299 |  1.88244 | 0.59 | 0.0031  |                   |
| 9711751 |  1.71153 | 0.79 | 5.3e-05 | b                 |
| 9714358 |  6.4742  | 0.2  | 0.014   | brc               |
| 9715925 |  6.30828 | 0.77 | 0.22    | b                 |
| 9777062 | 19.23    | 0.94 | 0.36    |                   |
| 9850387 |  2.7485  | 0.93 | 0.0028  | b                 |

Note. - Systems visually identified as exhibiting eclipse timing variations (ETVs). Note the period, mass ratio, and eccentricity included in this table correspond to maximum-likelihood solutions and not the median posterior values. al. (2016), Rappaport et al. (2013), Conroy et al. (2014), and Orosz y Triple candidates previously identified by Borkovits et (2015) contain flags 'b', 'r', 'c', 'o', respectively.

Table 6. Mass Ratios from APOGEE RVs

|     KIC |   P (d) |   GLYPH<1> P (d) |       Q |   GLYPH<1> Q |
|---------|---------|------------------|---------|--------------|
| 5284133 |  8.7778 |         0.0178   | 0.5358  |      0.011   |
| 5460835 | 21.47   |         2.64     | 0.8285  |      0.0346  |
| 6610219 | 11.3007 |         0.000301 | 0.96134 |      0.00248 |

Table 7. RV-derived EB Mass Values from Literature

| KIC               | P ' d '               | M 1 ' M GLYPH<12> '   | GLYPH<1> M 1 ' M GLYPH<12> '   | M 2 ' M GLYPH<12> '   | GLYPH<1> M 2 ' M GLYPH<12> '   | Provenance                                                   |
|-------------------|-----------------------|-----------------------|--------------------------------|-----------------------|--------------------------------|--------------------------------------------------------------|
| 2305372           | 1.40469145            | 1.2                   | 0.1                            | 0.62                  | 0.04                           | Matson et al. (2017)                                         |
| 3241619           | 1.70334707            | 1.24                  | 0.04                           | 0.86                  | 0.02                           | Matson et al. (2017)                                         |
| 3440230           | 2.88110031            | 1.6                   | 0.1                            | 0.4                   | 0.03                           | Matson et al. (2017)                                         |
| 4285087           | 4.4860314             | 1.137                 | 0.013                          | 1.103                 | 0.014                          | Clark Cunningham et al. (2019)                               |
| 4544587           | 2.1891143             | 1.69                  | 0.1                            | 1.42                  | 0.09                           | Matson et al. (2017)                                         |
| 4574310           | 1.30622004            | 1.38                  | 0.06                           | 0.31                  | 0.01                           | Matson et al. (2017)                                         |
| 4665989           | 2.248067589           | 1.77                  | 0.09                           | 1.32                  | 0.05                           | Matson et al. (2017)                                         |
| 4848423           | 3.0036461             | 1.22                  | 0.05                           | 1.08                  | 0.04                           | Matson et al. (2017)                                         |
| 4851217           | 2.47028283            | 1.43                  | 0.05                           | 1.55                  | 0.05                           | Matson et al. (2017)                                         |
| 4862625           | 20.000214             | 1.47                  | 0.08                           | 0.37                  | 0.035                          | Kostov et al. (2013)                                         |
| 5284133           | 8.7845758             | -                     | -                              | -                     | -                              | F. Kiefer, independent analysis                              |
| 5285607           | 3.89940111            | 1.557                 | 0.038                          | 1.346                 | 0.033                          | Clark Cunningham et al. (2019)                               |
| 5444392           | 1.51952889            | 1.17                  | 0.01                           | 1.19                  | 0.1                            | Matson et al. (2017)                                         |
| 5460835           | 21.5392662            | -                     | -                              | -                     | -                              | F. Kiefer, independent analysis                              |
| 5473556           | 11.258818             | 1.2207                | 0.0112                         | 0.9678                | 0.0039                         | Kostov et al. (2016)                                         |
| 5513861           | 1.51020825            | 1.5                   | 0.04                           | 1.32                  | 0.03                           | Matson et al. (2017)                                         |
| 5738698           | 4.80877396            | 1.52                  | 0.03                           | 1.44                  | 0.02                           | Matson et al. (2017)                                         |
| 5786154           | 2.00827091            | 1.06                  | 0.06                           | 1.02                  | 0.04                           | Gaulme et al. (2016)                                         |
| 6131659           | 17.5278276            | 0.942                 | 0.010                          | 0.703                 | 0.008                          | Clark Cunningham et al. (2019)                               |
| 6206751           | 1.24534281            | 1.5                   | 0.05                           | 0.198                 | 0.007                          | Matson et al. (2017)                                         |
| 6525196           | 3.42059774            | 1.0351                | 0.0055                         | 0.9712                | 0.0039                         | Heglyph[suppress]miniak l et al. (2017)                      |
| 6610219           | 21.5392662            | -                     | -                              | -                     | -                              | F. Kiefer, independent analysis                              |
| 6762829           | 18.79537              | 0.949                 | 0.06                           | 0.2492                | 0.06                           | Orosz et al. (2012b)                                         |
| 6778289           | 30.1301383            | 1.51                  | 0.022                          | 1.091                 | 0.018                          | Clark Cunningham et al. (2019)                               |
| 6781535           | 9.12208635            | 1.01                  | 0.03                           | 1.03                  | 0.03                           | Clark Cunningham et al. (2019)                               |
| 6864859           | 40.8778419            | 1.354                 | 0.029                          | 1.411                 | 0.028                          | Clark Cunningham et al. (2019)                               |
| 7037405           | 207.1083              | 1.25                  | 0.04                           | 1.14                  | 0.02                           | Gaulme et al. (2016)                                         |
| 7605600           | 3.32619385            | 0.53                  | 0.02                           | 0.17                  | 0.01                           | Han et al. (2019)                                            |
| 7821010           | 24.238235             | 1.289                 | 0.015                          | 1.231                 | 0.015                          | Heglyph[suppress]miniak l et al. (2017)                      |
| 7943602           | 14.69199              | 1.0                   | 0.1                            | 0.78                  | 0.05                           | Gaulme et al. (2016)                                         |
| 8262223           | 1.61301476            | 1.96                  | 0.006                          | 0.21                  | 0.001                          | Guo et al. (2017)                                            |
| 8552540           | 1.06193426            | 1.32                  | 0.03                           | 1.04                  | 0.02                           | Matson et al. (2017)                                         |
| 8560861           | 31.9732937            | 1.93                  | 0.12                           | 1.06                  | 0.08                           | Borkovits et al. (2014)                                      |
| 8572936           | 27.7958103            | 1.0479                | 0.0033                         | 1.0208                | 0.0022                         | Welsh et al. (2012)                                          |
| 8823397           | 1.5065037             | 2.1                   | 0.2                            | 0.21                  | 0.02                           | Matson et al. (2017)                                         |
| 9159301           | 3.04477215            | 1.61                  | 0.08                           | 0.4                   | 0.02                           | Matson et al. (2017)                                         |
| 9246715           | 171.2768599           | 2.149                 | 0.007                          | 2.171                 | 0.007                          | Gaulme et al. (2016)                                         |
| 9602595           | 3.5565129             | 3.0                   | 0.1                            | 0.60                  | 0.03                           | Matson et al. (2017)                                         |
| 9632895           | 27.322037             | 0.934                 | 0.01                           | 0.1938                | 0.002                          | Welsh et al. (2014)                                          |
| 9641031           | 2.17815425            | 1.2041                | 0.0076                         | 0.9498                | 0.0046                         | Heglyph[suppress]miniak l et al. (2017)                      |
| 9837578           | 20.733666             |                       |                                |                       | 0.0045                         | Welsh et al. (2012)                                          |
| 9970396           | 235.2985              | 0.8877 1.14           | 0.0053 0.03                    | 0.8094 1.02           | 0.02                           | Gaulme et al. (2016)                                         |
| 10001167          | 120.3903              | 0.81                  | 0.05                           | 0.79                  | 0.03                           | Gaulme et al. (2016)                                         |
| 10020423          | 7.44837695            | 1.043                 | 0.055                          | 0.362                 | 0.013                          | Orosz et al. (2012)                                          |
| 10031808          | 8.5896432             |                       |                                |                       | 0.013                          | Heglyph[suppress]miniak l et al. (2017)                      |
| 10156064          | 4.855936446           | 1.741 2.1             | 0.009 0.1                      | 1.798 1.44            | 0.08                           | Matson et al. (2017)                                         |
| 10191056          | 2.427494881           | 1.59                  | 0.32                           | 1.427                 | 0.036                          | Heglyph[suppress]miniak l et al. (2017)                      |
| 10581918          | 1.80186366            |                       |                                | 0.169                 | 0.009                          | Matson et al. (2017)                                         |
| 10619109          | 2.04516616            | 1.3 1.5               | 0.06 0.4                       | 0.31                  | 0.07                           | Matson et al. (2017)                                         |
| 10736223          | 1.105094186           | 1.6                   | 0.1                            | 0.35                  | 0.03                           | Matson et al. (2017)                                         |
| 10935310          | 4.128795224           | 0.92                  | 0.05                           | 0.5                   | 0.03                           | Han et al. (2017)                                            |
| 10987439          | 10.67459809           | 0.9862                |                                |                       | 0.0045                         | Heglyph[suppress]miniak l et al. (2017)                      |
| 11922782 12351927 | 3.512934275 10.116146 | 1.067                 | 0.0034 0.010                   | 1.4215 0.836          | 0.006                          | Heglyph[suppress]miniak l et al. (2017) Kostov et al. (2014) |
| 12644769          | 41.07922              | 0.82 0.6897           | 0.015 0.0035                   | 0.5423                | 0.008                          | Doyle et al. (2011)                                          |
|                   |                       |                       |                                | 0.20255               | 0.00066                        |                                                              |

Note. - Note for binary studies which contain multiple mass estimates (e.g., from asteroseismology and RV), we adopted the RV-only derived values for consistency. GLYPH&lt;1&gt; M 1 and GLYPH&lt;1&gt; M 2 are uncertainties associated with published mass values.

Figure 1. Maximum-likelihood (ML) solution for KIC 10031419, a G-dwarf binary in a GLYPH&lt;24&gt; 4 d, circular orbit. The upper panels show the normalized, detrended, and phase-folded primary and secondary eclipse flux data (black) and the light curve fit corresponding to the ML model (red), while the bottom left panel shows the SED fit. Data points with Kepler quality flags &gt; 8 are masked for visual clarity here (although they are not removed during fitting, see § 2.2.1). The ML parameters are reported for reference. Both the light curve and SED residuals belong to the mean and mode bin from Figure 6, and indicate a good fit to the data; the scatter during eclipse is consistent with starspot modulations.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00139v1-figure-44a235ed-60a1-4ba7-a7d9-5e4b13e96177.jpg)

## KIC 10198109 simultaneous

ml 1.18e+00

Icerr = 3.13e-04

m2 3.93e-01

isoerr = 4.78e-02

z0 = 1.93e-03

msum 1.57e+00

age = 9.07e+00

mrat = 3.34e-01

dist = 7.20e+02

rsum 1.49e+00

ebv 6.22e-01

rrat 3.22e-01

ho = 1.19e+02

rl = 1.13e+00

period 1.79e+01

r2 = 3.63e-01

tpe 1.46e+02

inc = 1.58e+00

esinw 2.77e-01

frat = 4.73e-03

ecosw 2.48e-01

b = 3.43e-01

8.77e-01

94 = 6.39e-01

Figure 2. Same as Figure 1, except for KIC 10198109, an eccentric, 18 d binary with Q GLYPH&lt;24&gt; 0 3 : . The in-eclipse scatter is small, . 1 ppt, and again likely due to starspots.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00139v1-figure-e3a57f6c-300f-42e0-a567-6b9f29511081.jpg)

ml =8.79e-01

Icerr = 1.34e-02

m2 5.72e-01

isoerr = 4.14e-02

z0 = 5.52e-03

msum 1.45e+00

age = 9.89e+00

mrat = 6.51e-01

dist = 1.28e+03

rsum 1.56e+00

ebv =

2.68e-01

rrat 5.63e-01

ho = 1.19e+02

rl = 1.00e+00

period 2.73e+01

r2 = 5.63e-01

tpe 1.44e+02

inc = 1.57e+00

esinw 4.10e-01

frat = 5.96e-02

ecosw 1.22e-01

b = 2.47e-01

5,88e-01

2.27e-01

5.46e-01

Figure 3. Same as Figure 1, but for KIC 12356914, a pair of low mass stars in an eccentric, 27 d orbit around each other. The strong eclipse timing variations present in primary and secondary eclipse pull GLYPH&lt;27&gt; sys LC ; toward very large values, and indicate the presence of a third companion in a wider orbit around the central binary, consistent with the more rigorous analysis of Borkovits et al. (2016).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00139v1-figure-9f4a964e-e683-44a0-b7dc-c7a0be265774.jpg)

## KIC 10619109 simultaneous

Figure 4. Same as Figure 1, but for KIC 10619109, a 2 d, nearly circular orbit binary. Note the relatively large light curve scatter out of eclipse and 0.2 mag scatter in the SED residuals. This is an example of a marginally good fit, but with anomalously young ages, relatively high mass components, and morph parameter &gt; 0.5, indicative of an Algol-type system (see § 4.2.3 for further discussion). The outlier data points are an artifact of poor polynomial fitting, due to missing data near/during a particular eclipse.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00139v1-figure-1e99a0e9-9d8f-4239-a42a-2eddebcb7966.jpg)

Figure 5. Same as Figure 1, but for KIC 12644769, a system in which two M dwarfs orbit each other every 41 days in a slightly eccentric fashion. The scatter in-eclipse is relatively small at &lt; 1 ppt, but nevertheless exhibit interesting phenomena. Residuals during primary eclipse are consistent with starspot modulations (on the primary star), while residuals at secondary ingress and egress indicate small ETVs by a small tertiary component. Indeed, KIC 12644769, aka Kepler-16, is a known circumbinary planet host (Doyle et al. 2011).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00139v1-figure-27e56fbb-b422-4e49-862b-9b53e8453712.jpg)

Figure 6. The median absolute data-model residual distribution (black curve) for our EB fits in log space; we decompose the distributions between morph &lt; 0.5 (blue) and morph &gt; 0.5 (orange) systems and stack them vertically. The majority of systems modeled exhibit residual levels comparable to light curve 'jitter,' quantified by the median absolute difference (MAD) in out-of-eclipse flux (red dashed curve). However, this proxy for LC noise does not capture correlated noise well, so light curves which exhibit e.g., large Doppler or ellipsoidal amplitudes, starspot variation, quasi-periodic stellar variability, or third light dilution that varies within each quarter, will have underestimated noise values. This contributes to the difference between in-eclipse model residuals and inherent light curve jitter distributions; systems with large eclipse timing variations, starspot or ellipsoidal variations compose the tail of large residuals ( &amp; 0 01 : ). Indeed, morph &gt; 0.5, e.g., short period EBs likely to exhibit ellipsoidal variations and/or have undergone interactions, have larger residuals relative to morph &lt; 0.5 systems.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00139v1-figure-03f7b0ab-f1b2-4ba0-85f7-f7e2f6d4675f.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00139v1-figure-779d7250-7f2a-4744-9914-fd123fad292b.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00139v1-figure-b5b98acd-1ed9-42c5-b942-4ee3ba8ed86b.jpg)

Figure 7. Distribution of orbital parameters: period, inclination, eccentricity, and orbital separation relative to stellar size, for EBs in our sample. The period distribution peaks around 3 days, and follows a gamma distribution in shape. The inclination distribution peaks near edge-on and declines exponentially. Note that our fits allow inclinations from 0 GLYPH&lt;14&gt; to 180 , but here we wish to illustrate the GLYPH&lt;14&gt; distribution of edge-on vs. grazing systems, so we folded the data about 90 degrees. The eccentricity distribution is bimodal, peaking around log e GLYPH&lt;24&gt; GLYPH&lt;0&gt; 1 6 : ; GLYPH&lt;0&gt; 0 5 : corresponding to e GLYPH&lt;25&gt; 0.03, 0.3. We also plot the distribution of orbital semi-major axis relative to stellar radius, which to first order is a proxy of light curve morphology; about 3/4 of our sample are have ' R 1 + R 2 'GLYPH&lt;157&gt; a &lt; 0 2 : .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00139v1-figure-beb17c8e-db6b-4410-a932-04bb282316d3.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00139v1-figure-5704d9d0-2016-4e97-8d65-b43e1c575d55.jpg)

Figure 8. Morphology (top), inclination (middle), and eccentricity (bottom) of EBs in our sample as a function of their orbital period.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00139v1-figure-5baac668-8915-4524-aeb4-2df1c2b1d4ec.jpg)

Figure 9. Distribution of stellar parameters: age, log g , and mass ratio Q , for EBs in our sample. The age distribution (upper left) is bimodal, such that the majority of EBs contain GLYPH&lt;24&gt; billion-year-old stars, with a small but notable excess of young ( &lt; 100 Myr) stars (see text for discussion). The surface gravity distribution (upper right) peaks around log g GLYPH&lt;24&gt; 4 3 : , indicating the prevalence of main sequence binaries in our sample. The mass ratio distribution (lower left) shows an increasing slope toward similar-mass binaries, consistent with observations of solar-type binaries in the field (Raghavan et al. 2010). Note that for systems where Q &gt; 1 (light grey), we invert the mass ratio such that Q = M 1 GLYPH&lt;157&gt; M 2 to keep values in range [0, 1]. The mass ratio distribution for Q GLYPH&lt;20&gt; 1 (lower right) is relatively uniform across P &lt; 4 , P = 4 GLYPH&lt;0&gt; 10 , and P &gt; 10 d binaries.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00139v1-figure-776a0ba7-8a9e-4f84-bc50-bfacd587f8af.jpg)

Figure 10. EBs from our sample, with temperature and absolute Kepler magnitude of each stellar component, plotted on the HR diagram; the color bar denotes inferred metallicity z of the system. While there is a small population of sub-giants starting to turn off the main sequence, there is an absence of red giant branch. The lack of giants may be due to a combination of isochrone fitting bias and Kepler target selection (see text for discussion).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00139v1-figure-de3139b6-cdf0-4bf2-845e-c04e5d1cc46a.jpg)

## 30 D. Windemuth et al.

Figure 11. Radial velocity solutions in an independent analysis to three SB2 systems in our sample, using the method of Kiefer et al. (2018).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00139v1-figure-3511e318-ac11-4f7d-98ec-25fd657de085.jpg)

Figure 12. Comparison of RV-derived mass ratio values as reported in literature vs. KEBLAT values using the SED+LC fitting method. The dashed black line denotes 1:1 relationship, with light grey regions representing GLYPH&lt;6&gt; 0.1 uncertainty in Q . In general, there is broad agreement with literature values. Binaries which are not well described by stellar isochrones (represented by lighter alpha values), e.g., red giants and binaries that exchange(d) mass, give discrepant mass ratio values.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00139v1-figure-6d37afef-6705-4b54-8314-37b9c0b75061.jpg)

Figure 13. Comparison of absolute mass values between KEBLAT and RV-derived values from literature, after removing EBs with red giant components and morph &gt; 0.5 (see § 4.1.1). The closed and open circles represent primary and secondary components of the binary. The dashed black line denotes 1:1 relationship, with light grey regions representing 15% of the mass uncertainties. The photometric masses show good agreement with RV values.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00139v1-figure-4fd5bcb4-c4ed-4493-a7a0-6e4e9ddc6f83.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00139v1-figure-041c1f5f-fe0a-45b7-aa35-ef58d75cf7eb.jpg)

/circledot

Figure 14. Comparison of mass values from our analysis and that of Devor et al. (2008) for a small number of overlap binaries. Both studies used light curve and SED data with stellar isochrones to derive masses, although predicated on different data and details of each method are different. The mass values show relatively good agreement.

0

Figure 15. Comparison of eccentricity, inclination, and radii values from KEBLAT and various literature studies, which show overall good agreement. The values are collected from the same RV studies from mass comparison, plotted in Fig. 12, and we supplement additional eccentricity estimates from (Kjurkchieva et al. 2017). While the masses of red giant components (Gaulme et al. 2016) were poorly inferred, the absolute radii show good fidelity to literature values. There is also greater scatter among inclination values from Matson et al. (2017) (blue), however these were fixed to be Slawson et al. (2011) neural network inference values based on phenomenological light curve modeling rather than physical model, which may explain some of the discrepancies.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00139v1-figure-83c7638b-e59d-4e3e-83f4-7f29b15555f1.jpg)

Figure 16. Comparison of parameters as derived in this paper and values presented in Armstrong et al. (2014) for temperature of the primary star (top left), primary radius normalized by distance to system (top right), temperature ratio (bottom left), and radius ratio (bottom right). While T 1 and R 1 GLYPH&lt;157&gt; d show bulk agreement, their temperature and radius ratios estimates are relatively crude and show significant scatter with respect to values derived in this work (see text for discussion).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00139v1-figure-9c87d21f-5272-4fe2-9e07-06406a587176.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00139v1-figure-dea98b07-2667-4758-ab30-c9a513c6f70e.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00139v1-figure-21889010-9678-4e2d-986a-28a08fdf6ac7.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00139v1-figure-f5c95e63-a246-4da3-bcd1-f31e01a99a95.jpg)

## 36 D. Windemuth et al.

## REFERENCES

| Abt, H. A. 1979, AJ, 84, 1591                                                                                                                                                                     |
|---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Ahn, C. P., Alexandroff, R., Allende Prieto, C., et al. 2012, ApJS, 203, 21 Angus, R., Morton, T., Aigrain, S., et al. 2018, MNRAS, 474, 2094.                                                    |
| Armstrong, D. J., Osborn, H. P., Brown, D. J. A., et al. 2014, MNRAS, 444, 1873                                                                                                                   |
| Badenes, C., Mazzola, C., Thompson, T. A., et al. 2018, ApJ, 854, 147                                                                                                                             |
| Bailer-Jones, C. A. L., Rybizki, J., Fouesneau, M., Mantelet, G., & Andrae, R. 2018, arXiv:1804.10121                                                                                             |
| Bastien, F. A., Stassun, K. G., Basri, G., & Pepper, J. 2013, Nature, 500, 427                                                                                                                    |
| Berger, T. A., Huber, D., Gaidos, E., & van Saders, J. L. 2018, arXiv:1805.00231                                                                                                                  |
| Birkby, J., Nefs, B., Hodgkin, S., et al. 2012, MNRAS, 426, 1507.                                                                                                                                 |
| Borkovits, T., Derekas, A., Fuller, J., et al. 2014, MNRAS, 443, 3068                                                                                                                             |
| Borkovits, T., Rappaport, S., Hajdu, T., & Sztakovics, J. 2015, MNRAS, 448, 946                                                                                                                   |
| Borkovits, T., Hajdu, T., Sztakovics, J., et al. 2016, MNRAS, 455, 4136                                                                                                                           |
| Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427,                                                                                                                                    |
| 127 Brown, T. M., Latham, D. W., Everett, M. E., & Esquerdo, G. A. 2011, AJ, 142, 112                                                                                                             |
| Cassisi, S. 2017, European Physical Journal Web of Conferences, 160, 04002                                                                                                                        |
| Castelli, F., & Kurucz, R. L. 2003, Modelling of Stellar Atmospheres, A20.                                                                                                                        |
| Chabrier, G., Gallardo, J., & Baraffe, I. 2007, A&A, 472, L17                                                                                                                                     |
| Chen, Y., Girardi, L., Bressan, A., et al. 2014, MNRAS, 444, 2525                                                                                                                                 |
| Claret, A., & Torres, G. 2018, ApJ, 859, 100.                                                                                                                                                     |
| Clark Cunningham, J. M., Rawls, M. L., Windemuth, D., et al. 2019, arXiv e-prints, arXiv:1907.00095                                                                                               |
| Conroy, K. E., Prˇa, s A., Stassun, K. G., et al. 2014, AJ, 147, 45 Cruz, P., Diaz, M., Birkby, J., et al. 2018, MNRAS, 476, 5253                                                                 |
| Cutri, R. M., & et al. 2014, VizieR Online Data Catalog,                                                                                                                                          |
| 2328,                                                                                                                                                                                             |
| Davenport, J. R. A. 2016, ApJ, 829, 23                                                                                                                                                            |
| del Burgo, C. and Allende Prieto, C. (2018) MNRAS,                                                                                                                                                |
| 479: 1953 J., 2005 ApJ, 628 411                                                                                                                                                                   |
| Devor, Devor, J., Charbonneau, D., O'Donovan, F.                                                                                                                                                  |
| T., Mandushev, G., & Torres, G. 2008, AJ,                                                                                                                                                         |
| 135, 850 Dotter, A. 2016, ApJS, 222, 8                                                                                                                                                            |
| Doyle, L. R., Carter, J. A., Fabrycky, D. C., et al. 2011, Science, 333, 1602 Duchˆne, e G., & Kraus, A. 2013, ARA&A, 51, 269                                                                     |
| Duquennoy, A., & Mayor, M. 1991, A&A, 248, 485 Eldridge, J. J., Stanway, E. R., Xiao, L., et al. 2017, Publ. Astron. Soc.                                                                         |
| Australia, 34, e058 Everett, M. E., Howell, S. B., & Kinemuchi, K. 2012, PASP, 124, 316                                                                                                           |
| Faigler, S., & Mazeh, T. 2011, MNRAS, 415, 3921 Faigler, S., Tal-Or, L., Mazeh, T., Latham, D. W., & Buchhave, L. A. 2013, ApJ, 815,                                                              |
| 771, 26                                                                                                                                                                                           |
| Faigler S., Kull I., Mazeh T., Kiefer F., Latham D. W., Bloemen S.,                                                                                                                               |
| 2015, ApJ, 26 Feiden, G. A., & Chaboyer, B. 2014, ApJ, 789, 53                                                                                                                                    |
| Fitzpatrick, E. L. 1999, PASP, 111, 63                                                                                                                                                            |
| Fleming, D. P., Barnes, R., Graham, D. E., Luger, R., &                                                                                                                                           |
| Quinn, T. 86                                                                                                                                                                                      |
| R. 2018, ApJ, 858, Fleming, D. P., Barnes, R., Davenport, J. R. A., et al. 2019, arXiv e-prints , arXiv:1903.05686. Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 |
| Gaulme, P., McKeever, J., Jackiewicz, J., et al. 2016, ApJ, 832, 121 Gazeas, K. D., & Niarchos, P. G. 2006, MNRAS, 370, L29                                                                       |
| Gies, D. R., Williams, S. J., Matson, R. A., et al. 2012, AJ, 143, 137                                                                                                                            |
| Giles, H. A. C., Collier Cameron, A., & Haywood, R. D. 2017, MNRAS,                                                                                                                               |
| 472, 1618 Gilliland, R. L., Chaplin, W. J., Jenkins, J. M., Ramsey, L. W., & Smith, J. C.                                                                                                         |
| 2015, AJ, 150, 133 Guo, Z., Gies, D. R., Matson, R. A., & Garc´a ı Hern´ndez, a A. 2016, ApJ, 826, 69 Guo, Z., Gies, D. R., Matson, R. A., et al. 2017, ApJ, 837, 114                             |
| Halbwachs J.-L., et al., 2016, MNRAS, 455, 3303 Han, E., Muirhead, P. S., Swift, J. J., et al. 2017, AJ, 154, 100 Han, E., Muirhead, P. S., & Swift, J. J. 2019, arXiv e-prints,                  |
| Hauschildt P. H., Allard F., Baron E., 1999, ApJ, 512, 377                                                                                                                                        |
| arXiv:1907.07180 Heglyph[suppress]miniak, l K. G., Ukita, N., Kambe, E., et al. 2017, MNRAS, 468, 1726                                                                                            |
| Hodgson, R. M., Bailey, D. G., Naylor, M. J., Ng, A. L. M., McNeil, S. J.,                                                                                                                        |
| Image Vision Comput., Hogg, D. W., Bovy, J., & Lang, D. 2010, arXiv:1008.4686                                                                                                                     |
| 1985, Ibanoˇlu, g C., Soydugan, F., Soydugan, E., & Dervi¸oˇlu, s g A. 2006, MNRAS, 373,                                                                                                          |
| 435 Izzard, R. G., & Halabi, G. M. 2018, arXiv:1808.06883                                                                                                                                         |
| Jørgensen, B. R., & Lindegren, L. 2005, A&A, 436, 127 Kennedy, G. M., & Wyatt, M. C. 2012, MNRAS, 426, 91                                                                                         |
| Kesseli, A. Y., Muirhead, P. S., Mann, A. W., & Mace, G. 225                                                                                                                                      |
| 2018, AJ, 155, Kiefer, F., Halbwachs, J.-L., Arenou, F., et al. 2016, MNRAS, 458, 3272 Kiefer F., et al., 2018, MNRAS, 474, 731                                                                   |
| Kipping, D. M. 2010, MNRAS, 408, 1758. Kipping, D. M. 2013, MNRAS, 435, 2152.                                                                                                                     |
| B., Conroy, K., Prˇa, s A., et al. 2016, AJ,                                                                                                                                                      |
| 151,                                                                                                                                                                                              |
| Kirk, 68                                                                                                                                                                                          |

Kjurkchieva, D., Vasileva, D., &amp; Dimitrov, D. 2016, AJ, 152, 189

- Kjurkchieva, D., Vasileva, D., &amp; Atanasova, T. 2017, AJ, 154, 105
- Kostov, V. B., McCullough, P. R., Hinse, T. C., et al. 2013, ApJ, 770, 52
- Kostov, V. B., McCullough, P. R., Carter, J. A., et al. 2014, ApJ, 784, 14
- Kostov, V. B., Orosz, J. A., Welsh, W. F., et al. 2016, ApJ, 827, 86
- Kraus, A. L., Tucker, R. A., Thompson, M. I., Craine, E. R., &amp; Hillenbrand, L. A. 2011, ApJ, 728, 48
- Kruse, E., &amp; Agol, E. 2014, Science, 344, 275
- Lurie, J. C., Vyhmeister, K., Hawley, S. L., et al. 2017, AJ, 154, 250
- Majewski, S. R., Schiavon, R. P., Frinchaboy, P. M., et al. 2017, AJ, 154, 94.
- Mandel, K., &amp; Agol, E. 2002, ApJ, 580, L171
- Martin, D. V., Mazeh, T., &amp; Fabrycky, D. C. 2015, MNRAS, 453, 3554
- Mathur, S., Huber, D., Batalha, N. M., et al. 2017, ApJS, 229, 30
- Matijeviˇ, G., Prˇa, A., Orosz, J. A., et al. 2012, AJ, 143, 123 c s
- Matson, R. A., Gies, D. R., Guo, Z., &amp; Orosz, J. A. 2016, AJ, 151, 139
- Matson, R. A., Gies, D. R., Guo, Z., &amp; Williams, S. J. 2017, AJ, 154, 216
- Mor´ e, Jorge J. 'The Levenberg-Marquardt algorithm: implementation and theory.' Numerical analysis. Springer, Berlin, Heidelberg, 1978. 105-116.
- Morton, T. D. 2012, ApJ, 761, 6
- Newville, M., Stensitzki, T., Allen, D. B., &amp; Ingargiola, A. 2014, lmfit, v.0.9.2, Zenodo, doi:10.5281/zenodo.11813˝ u
- North, T. S. H., Chaplin, W. J., Gilliland, R. L., et al. 2017, MNRAS, 465, 1308
- Orosz, J. A., &amp; Hauschildt, P. H. 2000, A&amp;A, 364, 265
- Orosz, J. A., Welsh, W. F., Carter, J. A., et al. 2012, Science, 337, 1511
- Orosz, J. A., Welsh, W. F., Carter, J. A., et al. 2012b, ApJ, 758, 87
- Orosz, J. A. 2015, Living Together: Planets, Host Stars and Binaries, 496, 55
- Prˇ sa, A., Batalha, N., Slawson, R. W., et al. 2011, AJ, 141, 83
- Prˇ sa, A., Conroy, K. E., Horvat, M., et al. 2016, ApJS, 227, 29
- Price-Whelan, A. M., &amp; Goodman, J. 2018, arXiv:1804.06841
- Raghavan, D., McAlister, H. A., Henry, T. J., et al. 2010, ApJS, 190, 1
- Rappaport, S., Deck, K., Levine, A., et al. 2013, ApJ, 768, 33
- Rucinski, S. M. 2001, AJ, 122, 1007
- Schlafly, E. F., &amp; Finkbeiner, D. P. 2011, ApJ, 737, 103
- Schwamb, M. E. et al. 2013, ApJ, 768, 127
- Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163
- Slawson, R. W., Prˇ sa, A., Welsh, W. F., et al. 2011, AJ, 142, 160
- Southworth, J., Maxted, P. F. L., &amp; Smalley, B. 2005, A&amp;A, 429, 645.
- Stassun, K. G., Hebb, L., Lopez-Morales, M., &amp; Prsa, A. 2009, arXiv:0902.2548
- Stassun, K. G., Kratter, K. M., Scholz, A., &amp; Dupuy, T. J. 2012, ApJ, 756, 47
- Stassun, K. G., &amp; Torres, G. 2016, AJ, 152, 180
- Tegmark, M., Strauss, M. A., Blanton, M. R., et al. 2004, Phys. Rev. D, 69, 103501
- Tokovinin, A., Thomas, S., Sterzik, M., &amp; Udry, S. 2006, A&amp;A, 450, 681
- Tokovinin, A. 2008, MNRAS, 389, 925
- Van Eylen, V., Winn, J. N., &amp; Albrecht, S. 2016, ApJ, 824, 15
- Welsh, W. F., Orosz, J. A., Carter, J. A., et al. 2012, Nature, 481, 475
- Welsh, W. F., Orosz, J. A., Carter, J. A., &amp; Fabrycky, D. C. 2014, Proceedings of the International Astronomical Union, vol. 293, pp. 125-32
- Wilson, R. E., &amp; Devinney, E. J. 1971, ApJ, 166, 605
- Wyithe, J. S. B., &amp; Wilson, R. E. 2002, ApJ, 571, 293
- Yakut, K., &amp; Eggleton, P. P. 2005, ApJ, 629, 1055
- Zasowski, G., Johnson, J. A., Frinchaboy, P. M., et al. 2013, AJ, 146, 81
- Zucker S., Mazeh T., 1994, ApJ, 420, 806
- Zucker S., Mazeh T., Santos N. C., Udry S., Mayor M., 2003, A&amp;A, 404, 775
- Zucker, S., Mazeh, T., &amp; Alexander, T. 2007, ApJ, 670, 1326

This paper has been typeset from a T E X/LT E X file prepared by the author. A
Diana Windemuth Eric Agol Aleezah Ali Flavien Kiefer
Modeling Kepler Eclipsing Binaries: Homogeneous Inference of Orbital & Stellar Properties
We report on the properties of eclipsing binaries from the Kepler mission
with a newly developed photometric modeling code, which uses the light curve,
spectral energy distribution of each binary, and stellar evolution models to
infer stellar masses without the need for radial velocity measurements. We
present solutions and posteriors to orbital and stellar parameters for 728
systems, forming the largest homogeneous catalogue of full Kepler binary
parameter estimates to date. Using comparisons to published radial velocity
measurements, we demonstrate that the inferred properties (e.g., masses) are
reliable for well-detached main-sequence binaries, which make up the majority
of our sample. The fidelity of our inferred parameters degrades for a subset of
systems not well described by input isochrones, such as short-period binaries
that have undergone interactions, or binaries with post-main sequence
components. Additionally, we identify 35 new systems which show evidence of
eclipse timing variations, perhaps from apsidal motion due to binary tides or
tertiary companions. We plan to subsequently use these models to search for and
constrain the presence of circumbinary planets in Kepler eclipsing binary
systems.