The Transiting Exoplanet Survey Satellite (TESS) will be conducting a nearly all-sky photometric survey over two years, with a core mission goal to discover small transiting exoplanets orbiting nearby bright stars. It will obtain 30-minute cadence observations of all objects in the TESS fields of view, along with 2-minute cadence observations of 200,000 to 400,000 selected stars. The choice of which stars to observe at the 2-min cadence is driven by the need to detect small transiting planets, which leads to the selection of primarily bright, cool dwarfs. We describe the catalogs assembled and the algorithms used to populate the TESS Input Catalog (TIC). We also describe a ranking system for prioritizing stars according to the smallest transiting planet detectable, and assemble a Candidate Target List (CTL) using that ranking. We discuss additional factors that affect the ability to photometrically detect and dynamically confirm small planets, and we note additional stellar populations of interest that may be added to the final target list. The TIC is available on the STScI MAST server, and an enhanced CTL is available through the Filtergraph data visualization portal system at the URL https://filtergraph.vanderbilt.edu/tess_ctl .

EarthFinder is a Probe Mission concept selected for study by NASA for input to the 2020 astronomy decadal survey. This study is currently active and a final white paper report is due to NASA at the end of calendar 2018. We are tasked with evaluating the scientific rationale for obtaining precise radial velocity (PRV) measurements in space, which is a two-part inquiry: What can be gained from going to space? What can't be done form the ground? These two questions flow down to these specific tasks for our study - Identify the velocity limit, if any, introduced from micro- and macro-telluric absorption in the Earth's atmosphere; Evaluate the unique advantages that a space-based platform provides to emable the identification and mitigation of stellar acitivity for multi-planet signal recovery.

We present a possible NASA key project using the iSHELL near-infrared high-resolution echelle spectrograph on the NASA Infrared Telescope Facility for precise radial velocity follow-up of candidate transiting exoplanets identified by the NASA TESS mission. We briefly review key motivations and challenges with near-infrared radial velocities. We present the current status of our preliminary radial velocity analysis from the first year on sky with iSHELL.

The Kilodegree Extremely Little Telescope (KELT) project has been conducting a photometric survey for transiting planets orbiting bright stars for over ten years. The KELT images have a pixel scale of ~23 arsec per pixel (similar to TESS) and a large point spread function, and the KELT reduction pipeline uses a weighted photometric aperture with radius 3 arcmin. At this angular scale, multiple stars are typically blended in the photometric apertures. In order to identify false positives and confirm transiting exoplanets, we have assembled a follow-up network (KELT-FUN) to conduct imaging with higher spatial resolution, cadence, and photometric precision than the KELT telescopes, as well as spectroscopic observations of the candidate host stars. The KELT-FUN team has followed-up over 1,600 planet candidates since 2011, resulting in more than 20 planet discoveries. We present an all-sky catalog of 1,081 bright stars (6<V<13) that show transit-like features in the KELT light curves. These stars were originally identified as planet candidates, but were subsequently determined to be astrophysical false positives (FPs) after photometric and/or spectroscopic follow-up observations. The remaining ~500 retired planet candidates have been classified as false alarms (instrumental or systematic noise). The KELT-FUN team continues to pursue KELT and other planet candidates and will eventually follow up certain classes of TESS candidates. The KELT FP catalog will help minimize the duplication of follow-up observations by current and future transit surveys such as TESS.

The NASA Exoplanet Exploration Program's SAG15 group has solicited, collected, and organized community input on high-level science questions that could be addressed with future direct imaging exoplanet missions and the type and quality of data answering these questions will require. Input was solicited through a variety of forums and the report draft was shared with the exoplanet community continuously during the period of the report development (Nov 2015 -- May 2017). The report benefitted from the input of over 50 exoplanet scientists and from multiple open-forum discussions at exoplanet and astrobiology meetings. The SAG15 team has identified three group of high-level questions, those that focus on the properties of planetary systems (Questions A1--A2), those that focus on the properties of individual planets (Questions B1--B4), and questions that relate to planetary processes (Questions C1--C4). The questions in categories A, B, and C require different target samples and often different observational approaches. For each questions we summarize the current body of knowledge, the available and future observational approaches that can directly or indirectly contribute to answering the question, and provide examples and general considerations for the target sample required.

Future generations of precise radial velocity (RV) surveys aim to achieve sensitivity sufficient to detect Earth mass planets orbiting in their stars' habitable zones. A major obstacle to this goal is astrophysical radial velocity noise caused by active areas moving across the stellar limb as a star rotates. In this paper, we quantify how stellar activity impacts exoplanet detection with radial velocities as a function of orbital and stellar rotational periods. We perform data-driven simulations of how stellar rotation affects planet detectability and compile and present relations for the typical timescale and amplitude of stellar radial velocity noise as a function of stellar mass. We show that the characteristic timescales of quasi-periodic radial velocity jitter from stellar rotational modulations coincides with the orbital period of habitable zone exoplanets around early M-dwarfs. These coincident periods underscore the importance of monitoring the targets of RV habitable zone planet surveys through simultaneous photometric measurements for determining rotation periods and activity signals, and mitigating activity signals using spectroscopic indicators and/or RV measurements at different wavelengths.

Given that low-mass stars have intrinsically low luminosities at optical wavelengths and a propensity for stellar activity, it is advantageous for radial velocity (RV) surveys of these objects to use near-infrared (NIR) wavelengths. In this work we describe and test a novel RV extraction pipeline dedicated to retrieving RVs from low mass stars using NIR spectra taken by the CSHELL spectrograph at the NASA Infrared Telescope Facility, where a methane isotopologue gas cell is used for wavelength calibration. The pipeline minimizes the residuals between the observations and a spectral model composed of templates for the target star, the gas cell, and atmospheric telluric absorption; models of the line spread function, continuum curvature, and sinusoidal fringing; and a parameterization of the wavelength solution. The stellar template is derived iteratively from the science observations themselves without a need for separate observations dedicated to retrieving it. Despite limitations from CSHELL's narrow wavelength range and instrumental systematics, we are able to (1) obtain an RV precision of 35 m/s for the RV standard star GJ 15 A over a time baseline of 817 days, reaching the photon noise limit for our attained SNR, (2) achieve ~3 m/s RV precision for the M giant SV Peg over a baseline of several days and confirm its long-term RV trend due to stellar pulsations, as well as obtain nightly noise floors of ~2 - 6 m/s, and (3) show that our data are consistent with the known masses, periods, and orbital eccentricities of the two most massive planets orbiting GJ 876. Future applications of our pipeline to RV surveys using the next generation of NIR spectrographs, such as iSHELL, will enable the potential detection of Super-Earths and Mini-Neptunes in the habitable zones of M dwarfs.

We present the results of a precise near-infrared (NIR) radial velocity (RV) survey of 32 low-mass stars with spectral types K2-M4 using CSHELL at the NASA IRTF in the $K$-band with an isotopologue methane gas cell to achieve wavelength calibration and a novel iterative RV extraction method. We surveyed 14 members of young ($\approx$ 25-150 Myr) moving groups, the young field star $\varepsilon$ Eridani as well as 18 nearby ($<$ 25 pc) low-mass stars and achieved typical single-measurement precisions of 8-15 m s$^{-1}$ with a long-term stability of 15-50 m s$^{-1}$. We obtain the best NIR RV constraints to date on 27 targets in our sample, 19 of which were never followed by high-precision RV surveys. Our results indicate that very active stars can display long-term RV variations as low as $\sim$ 25-50 m s$^{-1}$ at $\approx$ 2.3125 $\mu$m, thus constraining the effect of jitter at these wavelengths. We provide the first multi-wavelength confirmation of GJ 876 bc and independently retrieve orbital parameters consistent with previous studies. We recovered RV variability for HD 160934 AB and GJ 725 AB that are consistent with their known binary orbits, and nine other targets are candidate RV variables with a statistical significance of 3-5$\sigma$. Our method combined with the new iSHELL spectrograph will yield long-term RV precisions of $\lesssim$ 5 m s$^{-1}$ in the NIR, which will allow the detection of Super-Earths near the habitable zone of mid-M dwarfs.

We present the results of two 2.3 micron near-infrared radial velocity surveys to detect exoplanets around 36 nearby and young M dwarfs. We use the CSHELL spectrograph (R ~46,000) at the NASA InfraRed Telescope Facility, combined with an isotopic methane absorption gas cell for common optical path relative wavelength calibration. We have developed a sophisticated RV forward modeling code that accounts for fringing and other instrumental artifacts present in the spectra. With a spectral grasp of only 5 nm, we are able to reach long-term radial velocity dispersions of ~20-30 m/s on our survey targets.

Theoretical models and spectroscopic observations of newborn stars suggest that protoplantary disks have an inner "wall" at a distance set by the disk interaction with the star. Around T Tauri stars, the size of this disk hole is expected to be on a 0.1-AU scale that is unresolved by current adaptive optics imaging, though some model-dependent constraints have been obtained by near-infrared interferometry. Here we report the first measurement of the inner disk wall around a solar-mass young stellar object, YLW 16B in the {\rho} Ophiuchi star-forming region, by detecting the light travel time of the variable radiation from the stellar surface to the disk. Consistent time lags were detected on two nights, when the time series in H (1.6 {\mu}m) and K (2.2 {\mu}m) bands were synchronized while the 4.5 {\mu}m emission lagged by 74.5 +/- 3.2 seconds. Considering the nearly edge-on geometry of the disk, the inner rim should be 0.084 AU from the protostar on average, with an error of order 0.01 AU. This size is likely larger than the range of magnetospheric truncations, and consistent with an optically and geometrically thick disk front at the dust sublimation radius at ~1500 K. The widths of the cross-correlation functions between the data in different wavebands place possible new constraints on the geometry of the disk.

Detecting transient light curves (e.g., transiting planets) requires high precision data, and thus it is important to effectively filter systematic trends affecting ground based wide field surveys. We apply an implementation of the Trend Filtering Algorithm (TFA) (Kovacs et al. 2005) to the 2MASS calibration catalog and select Palomar Transient Factory (PTF) photometric time series data. TFA is successful at reducing the overall dispersion of light curves, however it may over filter intrinsic variables and increase "instantaneous" dispersion when a template set is not judiciously chosen. In an attempt to rectify these issues we modify the original literature TFA by including measurement uncertainties in its computation, including ancillary data correlated with noise, and algorithmically selecting a template set using clustering algorithms as suggested by various authors. This approach may be particularly useful for appropriately accounting for variable photometric precision surveys and/or combined data-sets. In summary, our contributions are to provide a MATLAB software implementation of TFA and a number of modifications tested on synthetics and real data, summarize the performance of TFA and various modifications on real ground based data sets (2MASS and PTF), and assess the efficacy of TFA and modifications using synthetic light curve tests consisting of transiting and sinusoidal variables. While the transiting variables test indicates that these modifications confer no advantage to transit detection, the sinusoidal variables test indicates potential improvements in detection accuracy.

The MINiature Exoplanet Radial Velocity Array (MINERVA) is a US-based observational facility dedicated to the discovery and characterization of exoplanets around a nearby sample of bright stars. MINERVA employs a robotic array of four 0.7 m telescopes outfitted for both high-resolution spectroscopy and photometry, and is designed for completely autonomous operation. The primary science program is a dedicated radial velocity survey and the secondary science objective is to obtain high precision transit light curves. The modular design of the facility and the flexibility of our hardware allows for both science programs to be pursued simultaneously, while the robotic control software provides a robust and efficient means to carry out nightly observations. In this article, we describe the design of MINERVA including major hardware components, software, and science goals. The telescopes and photometry cameras are characterized at our test facility on the Caltech campus in Pasadena, CA, and their on-sky performance is validated. New observations from our test facility demonstrate sub-mmag photometric precision of one of our radial velocity survey targets, and we present new transit observations and fits of WASP-52b -- a known hot-Jupiter with an inflated radius and misaligned orbit. The facility is now in the process of being relocated to its final destination at the Fred Lawrence Whipple Observatory in southern Arizona, and science operations will begin in 2015.

We investigate the possibility of constraining the sin i degeneracy of alpha Cen B b -- with orbital period P=3.24 d; a = 0.042 AU; msini = 1.1 M_earth -- to estimate the true mass of the newly reported terrestrial exoplanet in the nearest stellar system to our Sun. We present detailed numerical simulations of the dynamical stability of the exoplanet in the alpha Cen AB binary system for a range of initial inclinations, eccentricities, and semi-major axes. The system represents a benchmark case for the interplay of the Kozai mechanism, general relativistic and tidal forces. From our simulations, there is only a small boundary in initial inclinations and initial semi-major axes that result in the migration via the Kozai mechanism of alpha Cen B b to its present location. Inside this boundary, the planet orbit is stable for up to 1 Gyr against the Kozai mechanism, and outside this boundary the planet collides with alpha Cen B or is ejected. In our three simulations where the planet migrates in towards the star via the Kozai mechanism, the final inclination is 46-53 degrees relative to the AB orbital plane, lower than the initial inclination of 75 degrees in each case. We discuss inclination constraints from the formation of alpha Cen B b in situ at its present location, migration in a proto-planetary disk, or migration in resonance with additional planets. We conclude that alpha Cen B b probably has a mass of less than 2.7 M_earth, implying a likely terrestrial composition warranting future confirmation.

The final assembly of terrestrial planets occurs via massive collisions, which can launch copious clouds of dust that are warmed by the star and glow in the infrared. We report the real-time detection of a debris-producing impact in the terrestrial planet zone around a 35-million year-old solar analog star. We observed a substantial brightening of the debris disk at 3-5 {\mu}m, followed by a decay over a year, with quasi-periodic modulations of the disk flux. The behavior is consistent with the occurrence of a violent impact that produced vapor out of which a thick cloud of silicate spherules condensed that were ground into dust by collisions. These results demonstrate how the time domain can become a new dimension for the study of terrestrial planet formation.

[Abridged] The Study Analysis Group 8 of the NASA Exoplanet Analysis Group was convened to assess the current capabilities and the future potential of the precise radial velocity (PRV) method to advance the NASA goal to "search for planetary bodies and Earth-like planets in orbit around other stars.: (U.S. National Space Policy, June 28, 2010). PRVs complement other exoplanet detection methods, for example offering a direct path to obtaining the bulk density and thus the structure and composition of transiting exoplanets. Our analysis builds upon previous community input, including the ExoPlanet Community Report chapter on radial velocities in 2008, the 2010 Decadal Survey of Astronomy, the Penn State Precise Radial Velocities Workshop response to the Decadal Survey in 2010, and the NSF Portfolio Review in 2012. The radial-velocity detection of exoplanets is strongly endorsed by both the Astro 2010 Decadal Survey "New Worlds, New Horizons" and the NSF Portfolio Review, and the community has recommended robust investment in PRVs. The demands on telescope time for the above mission support, especially for systems of small planets, will exceed the number of nights available using instruments now in operation by a factor of at least several for TESS alone. Pushing down towards true Earth twins will require more photons (i.e. larger telescopes), more stable spectrographs than are currently available, better calibration, and better correction for stellar jitter. We outline four hypothetical situations for PRV work necessary to meet NASA mission exoplanet science objectives.

\We present the sixth catalog of Kepler candidate planets based on nearly 4 years of high precision photometry. This catalog builds on the legacy of previous catalogs released by the Kepler project and includes 1493 new Kepler Objects of Interest (KOIs) of which 554 are planet candidates, and 131 of these candidates have best fit radii <1.5 R_earth. This brings the total number of KOIs and planet candidates to 7305 and 4173 respectively. We suspect that many of these new candidates at the low signal-to-noise limit may be false alarms created by instrumental noise, and discuss our efforts to identify such objects. We re-evaluate all previously published KOIs with orbital periods of >50 days to provide a consistently vetted sample that can be used to improve planet occurrence rate calculations. We discuss the performance of our planet detection algorithms, and the consistency of our vetting products. The full catalog is publicly available at the NASA Exoplanet Archive.

We identify nine young stellar objects (YSOs) in the NGC 2264 star-forming region with optical {\em CoRoT} light curves exhibiting short-duration, shallow, periodic flux dips. All of these stars have infrared (IR) excesses that are consistent with their having inner disk walls near the Keplerian co-rotation radius. The repeating photometric dips have FWHM generally less than one day, depths almost always less than 15%, and periods (3<P<11 days) consistent with dust near the Keplerian co-rotation period. The flux dips vary considerably in their depth from epoch to epoch, but usually persist for several weeks and, in two cases, were present in data collected on successive years. For several of these stars, we also measure the photospheric rotation period and find that the rotation and dip periods are the same, as predicted by standard "disk-locking" models. We attribute these flux dips to clumps of material in or near the inner disk wall, passing through our line of sight to the stellar photosphere. In some cases, these dips are also present in simultaneous {\em Spitzer} IRAC light curves at 3.6 and 4.5 microns. We characterize the properties of these dips, and compare the stars with light curves exhibiting this behavior to other classes of YSO in NGC 2264. A number of physical mechanisms could locally increase the dust scale height near the inner disk wall, and we discuss several of those mechanisms; the most plausible mechanisms are either a disk warp due to interaction with the stellar magnetic field or dust entrained in funnel-flow accretion columns arising near the inner disk wall.

We present results from the "Weather on Other Worlds" Spitzer Exploration Science program to investigate photometric variability in L and T dwarfs, usually attributed to patchy clouds. We surveyed 44 L3-T8 dwarfs, spanning a range of $J-K_s$ colors and surface gravities. We find that 14/23 (61%; 95% confidence interval: 41%-78%) of our single L3-L9.5 dwarfs are variable with peak-to-peak amplitudes between 0.2% and 1.5%, and 5/16 (31%; 95% confidence interval: 14%-56%) of our single T0-T8 dwarfs are variable with amplitudes between 0.8% and 4.6%. After correcting for sensitivity, we find that 80% (95% confidence interval: 53%-100%) of L dwarfs vary by >0.2%, and 36% (95% confidence interval: 19%-52%) of T dwarfs vary by >0.4%. Given viewing geometry considerations, we conclude that photospheric heterogeneities causing >0.2% 3-5-micron flux variations are present on virtually all L dwarfs, and probably on most T dwarfs. A third of L dwarf variables show irregular light curves, indicating that L dwarfs may have multiple spots that evolve over a single rotation. Also, approximately a third of the periodicities are on time scales >10 h, suggesting that slowly-rotating brown dwarfs may be common. We observe an increase in the maximum amplitudes over the entire spectral type range, revealing a potential for greater temperature contrasts in T dwarfs than in L dwarfs. We find a tentative association (92% confidence) between low surface gravity and high-amplitude variability among L3-L5.5 dwarfs. Although we can not confirm whether lower gravity is also correlated with a higher incidence of variables, the result is promising for the characterization of directly imaged young extrasolar planets through variability.

We present the Coordinated Synoptic Investigation of NGC 2264, a continuous 30-day multi-wavelength photometric monitoring campaign on more than 1000 young cluster members using 16 telescopes. The unprecedented combination of multi-wavelength, high-precision, high-cadence, and long-duration data opens a new window into the time domain behavior of young stellar objects. Here we provide an overview of the observations, focusing on results from Spitzer and CoRoT. The highlight of this work is detailed analysis of 162 classical T Tauri stars for which we can probe optical and mid-infrared flux variations to 1% amplitudes and sub-hour timescales. We present a morphological variability census and then use metrics of periodicity, stochasticity, and symmetry to statistically separate the light curves into seven distinct classes, which we suggest represent different physical processes and geometric effects. We provide distributions of the characteristic timescales and amplitudes, and assess the fractional representation within each class. The largest category (>20%) are optical "dippers" having discrete fading events lasting ~1-5 days. The degree of correlation between the optical and infrared light curves is positive but weak; notably, the independently assigned optical and infrared morphology classes tend to be different for the same object. Assessment of flux variation behavior with respect to (circum)stellar properties reveals correlations of variability parameters with H$\alpha$ emission and with effective temperature. Overall, our results point to multiple origins of young star variability, including circumstellar obscuration events, hot spots on the star and/or disk, accretion bursts, and rapid structural changes in the inner disk.

Based on more than four weeks of continuous high cadence photometric monitoring of several hundred members of the young cluster NGC 2264 with two space telescopes, NASA's Spitzer and the CNES CoRoT (Convection, Rotation, and planetary Transits), we provide high quality, multi-wavelength light curves for young stellar objects (YSOs) whose optical variability is dominated by short duration flux bursts, which we infer are due to enhanced mass accretion rates. These light curves show many brief -- several hour to one day -- brightenings at optical and near-infrared (IR) wavelengths with amplitudes generally in the range 5-50% of the quiescent value. Typically, a dozen or more of these bursts occur in a thirty day period. We demonstrate that stars exhibiting this type of variability have large ultraviolet (UV) excesses and dominate the portion of the u-g vs. g-r color-color diagram with the largest UV excesses. These stars also have large Halpha equivalent widths, and either centrally peaked, lumpy Halpha emission profiles or profiles with blue-shifted absorption dips associated with disk or stellar winds. Light curves of this type have been predicted for stars whose accretion is dominated by Rayleigh-Taylor instabilities at the boundary between their magnetosphere and inner circumstellar disk, or where magneto-rotational instabilities modulate the accretion rate from the inner disk. Amongst the stars with the largest UV excesses or largest Halpha equivalent widths, light curves with this type of variability greatly outnumber light curves with relatively smooth sinusoidal variations associated with long-lived hot spots. We provide quantitative statistics for the average duration and strength of the accretion bursts and for the fraction of the accretion luminosity associated with these bursts.

The Kepler mission discovery of candidate transiting exoplanets (KOIs) enables a plethora of ensemble analysis of the architecture and properties of exoplanetary systems. We compare the observed transit durations of KOIs to a synthetic distribution generated from the known eccentricities of radial velocity (RV) discovered exoplanets. We find that the Kepler and RV distributions differ at a statistically significant level. We identify three related systematic trends that are likely due to errors in stellar radii, which in turn affect the inferred exoplanet radii and the distribution thereof, and prevent a valid analysis of the underlying ensemble eccentricity distribution. First, 15% of KOIs have transit durations >20% longer than the transit duration expected for an edge-on circular orbit, including 92 KOIs with transit durations >50% longer, when only a handful of such systems are expected. Second, the median transit duration is too long by up to ~25%. Random errors of <50% in the stellar radius are not adequate to account for these two trends, and they are present for all spectral types in the Kepler sample. We identify that incorrect estimates of stellar metallicity and extinction could account for these anomalies, rather than astrophysical effects such as eccentric exoplanets improbably transiting near apastron. Third, we find that the median transit duration is correlated with stellar radius, when no such trend is expected. All three effects are still present, although less pronounced, when considering only multiple transiting KOI systems which are thought to have a low false positive rate. Improved stellar parameters for KOIs are necessary for the validity of future ensemble tests of exoplanetary systems found by Kepler.

With the loss of two reaction wheels, the period of Kepler's ultra-high precision photometric performance is at an end. Yet Kepler retains unique capabilities impossible to replicate from the ground or with existing or future space missions. This White Paper calls for the use of Kepler to conduct a survey in the ecliptic plane to search for planet transits around stars at high galactic latitudes and to study star forming regions to investigate physics of very young stars not studied by Kepler in its prime mission. Even with reduced photometric precision, Kepler's 1 m aperture will enable it to survey faint M stars to find ice giants and Super Earths in Habitable Zone orbits.

We report observations of a possible young transiting planet orbiting a previously known weak-lined T-Tauri star in the 7-10 Myr old Orion-OB1a/25-Ori region. The candidate was found as part of the Palomar Transient Factory (PTF) Orion project. It has a photometric transit period of 0.448413 +- 0.000040 days, and appears in both 2009 and 2010 PTF data. Follow-up low-precision radial velocity (RV) observations and adaptive optics imaging suggest that the star is not an eclipsing binary, and that it is unlikely that a background source is blended with the target and mimicking the observed transit. RV observations with the Hobby-Eberly and Keck telescopes yield an RV that has the same period as the photometric event, but is offset in phase from the transit center by approximately -0.22 periods. The amplitude (half range) of the RV variations is 2.4 km/s and is comparable with the expected RV amplitude that stellar spots could induce. The RV curve is likely dominated by stellar spot modulation and provides an upper limit to the projected companion mass of M_p sin i_orb < 4.8 +- 1.2 M_Jup; when combined with the orbital inclination, i orb, of the candidate planet from modeling of the transit light curve, we find an upper limit on the mass of the planetary candidate of M_p < 5.5 +- 1.4 M_Jup. This limit implies that the planet is orbiting close to, if not inside, its Roche limiting orbital radius, so that it may be undergoing active mass loss and evaporation.

Periodic variability in young stellar objects (YSOs) offers indirect evidence for an active dynamical mechanism. Starspots, accretion, stellar companions, and disk veiling can contribute to the photometric variability of YSOs. As part of an ongoing study of the Rho Oph star forming region, we report the discovery of 92.6 day periodic variations for the Class I YSO YLW 16A, observed over a period of three years. A SED model was fit to available photometric data for the object. We propose a triple-system with an inner binary with a period of 93 days eclipsed by a warped circum-binary disk. The nature of the secondary is unconstrained and could be stellar or sub-stellar. We report the discovery of a tertiary companion at a projected separation of ~40 AU that could account for the circum-binary disk warp. This light curve and model is similar to the model we proposed for WL 4 in previous work. Understanding these systems may lead to insights about the nature of stellar evolution and planetary formation, and provide valuable benchmarks for future theoretical modeling and near- and mid-infrared synoptic surveys of YSOs.

The Kepler Eclipsing Binary Catalog (KEBC)describes 2165 eclipsing binaries identified in the 115 deg^2 Kepler Field based on observations from Kepler quarters Q0, Q1, and Q2. The periods in the KEBC are given in units of days out to six decimal places but no period errors are provided. We present the PEC (Period Error Calculator) algorithm which can be used to estimate the period errors of strictly periodic variables observed by the Kepler Mission. The PEC algorithm is based on propagation of error theory and assumes that observation of every light curve peak/minimum in a long time-series observation can be unambiguously identified. The PEC algorithm can be efficiently programmed using just a few lines of C computer language code. The PEC algorithm was used to develop a simple model which provides period error estimates for eclipsing binaries in the KEBC with periods less than 62.5 days. KEBC systems with periods >=62.5 days have KEBC period errors of about 0.0144 days. Periods and period errors of 7 eclipsing binary systems in the KEBC were measured using the NASA Exoplanet Archive Periodogram Service and compared to period errors estimated using the PEC algorithm.