We report the discovery of a Neptune-size (R_p = 3.87 +/- 0.06 R_Earth) transiting circumbinary planet, Kepler-1661 b, found in the Kepler photometry. The planet has a period of ~175 days and its orbit precesses with a period of only 35 years. The precession causes the alignment of the orbital planes to vary, and the planet is in a transiting configuration only ~7% of the time as seen from Earth. As with several other Kepler circumbinary planets, Kepler-1661 b orbits close to the stability radius, and is near the (hot) edge of habitable zone. The planet orbits a single-lined, grazing eclipsing binary, containing a 0.84 M_Sun and 0.26 M_Sun pair of stars in a mildly eccentric (e=0.11), 28.2-day orbit. The system is fairly young, with an estimated age of ~1-3 Gyrs, and exhibits significant starspot modulations. The grazing-eclipse configuration means the system is very sensitive to changes in the binary inclination, which manifests itself as a change in the eclipse depth. The starspots contaminate the eclipse photometry, but not in the usual way of inducing spurious eclipse timing variations. Rather, the starspots alter the normalization of the light curve, and hence the eclipse depths. This can lead to spurious eclipse depth variations, which are then incorrectly ascribed to binary orbital precession.

There is a well-known stellar parameter discrepancy for late K and M dwarfs, in that the observed radii and temperatures are often respectively larger and cooler than predicted by theory by several percent. In an on-going effort to elucidate this issue, we examine the double-lined Kepler eclipsing binary star system KIC 8736245. We supplement the near-continuous 4-year Kepler light curve with ground-based multicolor photometry from Mount Laguna Observatory and spectroscopy from the Hobby-Eberly Telescope. The binary has an edge-on, circular 5.07 day orbit with stellar masses equal to $0.987\pm0.009$ and $0.782\pm0.009\,\text{M}_\odot$ and radii of $1.311 \pm 0.006$ and $0.804 \pm 0.004\,\text{R}_\odot$, respectively, and an estimated age of 7-9 Gyr. We find that the stellar radii are consistent with theoretical models within the uncertainties, whereas the temperature of the secondary star is $\sim$6% cooler than predicted. An important aspect of this work is that the uncertainties derived from a single epoch (individual night of observations) underestimates the overall system parameter uncertainties due to the effect of the 1-4% fluctuations caused by stellar activity. Our error estimates come from the spread in parameters measured at 8 epochs. From the periodicities in the light curve and from the eclipse times, we measure candidate spin periods to be approximately 4.98 and 5.87 days for the primary and secondary star. Surprisingly, these imply super- and sub-synchronous rotation compared to the orbital period. Thus KIC 8736245 serves as an interesting case study for the exchange of angular momentum and general stellar astrophysics as stars in binaries evolve off the main sequence.

Recently there has been a renewed interest in the power-2 limb darkening law for modeling exoplanet transits. This law provides a better match to the intensities generated by spherical stellar atmosphere models than other 2-parameter laws. To help facilitate a wider use of the power-2 law we correct a minor error and, expanding on previous work, suggest a parametrization that can improve the sampling required by some numerical methods such as MCMC.

Of the nine confirmed transiting circumbinary planet systems, only Kepler-47 is known to contain more than one planet. Kepler-47 b (the "inner planet") has an orbital period of 49.5 days and a radius of about $3\,R_{\oplus}$. Kepler-47 c (the "outer planet") has an orbital period of 303.2 days and a radius of about $4.7\,R_{\oplus}$. Here we report the discovery of a third planet, Kepler-47 d (the "middle planet"), which has an orbital period of 187.4 days and a radius of about $7\,R_{\oplus}$. The presence of the middle planet allows us to place much better constraints on the masses of all three planets, where the $1\sigma$ ranges are less than $26\,M_{\oplus}$, between $7-43\,M_{\oplus}$, and between $2-5\,M_{\oplus}$ for the inner, middle, and outer planets, respectively. The middle and outer planets have low bulk densities, with $\rho_{\rm middle} < 0.68$ g cm$^{-3}$ and $\rho_{\rm outer} < 0.26$ g cm$^{-3}$ at the $1\sigma$ level. The two outer planets are "tightly packed," assuming the nominal masses, meaning no other planet could stably orbit between them. All of the orbits have low eccentricities and are nearly coplanar, disfavoring violent scattering scenarios and suggesting gentle migration in the protoplanetary disk.

We present here an efficient method for computing the visible flux for each body during a multi-body eclipsing event for all commonly used limb darkening laws. Our approach follows the idea put forth by Pal (2012) to apply Green's Theorem on the limb darkening integral, thus transforming the two-dimensional flux integral over the visible disk into a one-dimensional integral over the visible boundary. We implement this idea through an iterative process which combines a fast method for describing the visible boundary of each body with a fast numerical integration scheme to compute the integrals. For the two-body case, our method compares well in speed with both that of Mandel & Agol (2002) and that of Gimenez (2006a). The strength of the method is that it works for any number of spherical bodies, with a computational accuracy that is adjustable through the use of a tolerance parameter. Most significantly, the method offers two main advantages over previously used techniques: (i) it can employ a multitude of limb darkening laws, including all of the commonly used ones; (ii) it can compute the Rossiter-McLaughlin effect for rigid body rotation with an arbitrary orientation of the rotation axis, using any of these limb darkening laws. In addition, we can compute the Rossiter-McLaughlin effect for stars exhibiting differential rotation, using the quadratic limb darkening law. We provide the mathematical background for the method and explain in detail how to implement the technique with the help of several examples and codes which we make available.

KIC 9832227 is a contact binary whose 11 hr orbital period is rapidly changing. Based on the apparent exponential decay of its period, the two stars were predicted to merge in early 2022 resulting in a rare red nova outburst. Fortunately KIC 832227 was observed in 2003 as part of the NASA Ames pre-Kepler Vulcan Project to search for transiting exoplanets. We find that the Vulcan timing measurement does not agree with the previous exponential decay model. This led us to re-evaluate the other early epoch non-Kepler data sets, the Northern Sky Variability Survey (NSVS) and Wide Angle Search for Planets (WASP) survey. We find that the WASP times are in good agreement with the previous prediction, but the NSVS eclipse time differs by nearly an hour. The very large disagreement of the Vulcan and NSVS eclipse times with an exponentially decaying model forces us to reject the merger hypothesis. Although period variations are common in contact binaries, the physical cause of the period changes in KIC 9832227 remains unexplained; a third star scenario is unlikely. This study shows the data collected by the Vulcan photometer to be extremely valuable for extending the baseline for measurements of variable stars in the Kepler field.

The body of work presented here revolves around the investigation of the existence and nature of extra-solar planetary systems. The fitting of stellar radial velocity time series data is attempted by constructing a model to quantify the orbital properties of a star-planetary system. This is achieved with the Planetary Orbit Fitting Process (POFP). Though specific to the investigated problem, the POFP is founded on two separate, more general ideas. One is a Solver producing the gravitational dynamics of a Three-Body system by integrating its Newtonian equations of motion. The other is an independent optimisation scheme. Both have been devised using MATLAB. Applying the optimisation to the Solver results in a realistic Three-Body dynamics that best describes the radial velocity data under the model-specific orbital-observational constraints. Combining these aspects also allows for the study of dynamical instability derived from interaction, which is reaffirmed as a necessary criterion for evaluating the fit. The validity of POFP solutions with respect to the observations and other models is discussed in this context. The underlying generality and fundamental principles demonstrate a larger frame of operation where problems in Physics and Mathematics can be solved with a multitude of techniques.

In this letter we report the possible existence of a second planet in the transiting extrasolar planet system HD 17156 and its interactive dynamics with the previously known planet. The analysis is achieved through the \POFP\ optimization software which is based on a full integration of the system's multiple-body Newtonian equations of motion. The two-planet solution yields a significantly improved fit to the previously published radial velocities. The two planets are strongly interacting and exchange angular momentum in a 5:1 mean motion resonance, yet remain stable as they mutually excite orbital eccentricities and periastron advances.

We report the discovery of a new Kepler transiting circumbinary planet (CBP). This latest addition to the still-small family of CBPs defies the current trend of known short-period planets orbiting near the stability limit of binary stars. Unlike the previous discoveries, the planet revolving around the eclipsing binary system Kepler-1647 has a very long orbital period (~1100 days) and was at conjunction only twice during the Kepler mission lifetime. Due to the singular configuration of the system, Kepler-1647b is not only the longest-period transiting CBP at the time of writing, but also one of the longest-period transiting planets. With a radius of 1.06+/-0.01 RJup it is also the largest CBP to date. The planet produced three transits in the light-curve of Kepler-1647 (one of them during an eclipse, creating a syzygy) and measurably perturbed the times of the stellar eclipses, allowing us to measure its mass to be 1.52+/-0.65 MJup. The planet revolves around an 11-day period eclipsing binary consisting of two Solar-mass stars on a slightly inclined, mildly eccentric (e_bin = 0.16), spin-synchronized orbit. Despite having an orbital period three times longer than Earth's, Kepler-1647b is in the conservative habitable zone of the binary star throughout its orbit.

We present the discovery of KIC 9632895b, a 6.2 Earth-radius planet in a low-eccentricity, 240.5-day orbit about an eclipsing binary. The binary itself consists of a 0.93 and 0.194 solar mass pair of stars with an orbital period of 27.3 days. The plane of the planet's orbit is rapidly precessing, and its inclination only becomes sufficiently aligned with the primary star in the latter portion of the Kepler data. Thus three transits are present in the latter half of the light curve, but none of the three conjunctions that occurred during the first half of the light curve produced transits. The precession period is ~103 years, and during that cycle, transits are visible only ~8% of the time. This has the important implication that for every system like KIC 9632895 that we detect, there are ~12 circumbinary systems that exist but are not currently exhibiting transits. The planet's mass is too small to noticeably perturb the binary, consequently its mass is not measurable with these data; but our photodynamical model places a 1-sigma upper limit of 16 Earth masses. With a period 8.8 times that of the binary, the planet is well outside the dynamical instability zone. It does, however, lie within the habitable zone of the binary, and making it the third of ten Kepler circumbinary planets to do so.

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

KIC 6131659 is a long-period (17.5 days) eclipsing binary discovered by the Kepler mission. We analyzed six quarters of Kepler data along with supporting ground-based photometric and spectroscopic data to obtain accurate values for the mass and radius of both stars, namely M_1=0.922 +/- 0.007 M_sun, R_1=0.8800 +/- 0.0028 R_sun, and M_2=0.685 +/- 0.005 M_sun, R_2=0.6395 +/- 0.0061 R_sun. There is a well-known issue with low mass (M <<0.8 M_sun) stars (in cases where the mass and radius measurement uncertainties are smaller than two or three percent) where the measured radii are almost always 5 to 15 percent larger than expected from evolutionary models, i.e. the measured radii are all above the model isochrones in a mass-radius plane. In contrast, the two stars in KIC 6131659 were found to sit on the same theoretical isochrone in the mass-radius plane. Until recently, all of the well-studied eclipsing binaries with low-mass stars had periods less than about three days. The stars in such systems may have been inflated by high levels of stellar activity induced by tidal effects in these close binaries. KIC 6131659 shows essentially no evidence of enhanced stellar activity, and our measurements support the hypothesis that the unusual mass-radius relationship observed in most low-mass stars is influenced by strong magnetic activity created by the rapid rotation of the stars in tidally-locked, short-period systems. Finally, using short cadence data, we show that KIC 6131657 has one of the smallest measured non-zero eccentricities of a binary with two main sequence stars, where e cos omega=(4.57 +/-0.02)*10^-5.

We report the detection of Kepler-47, a system consisting of two planets orbiting around an eclipsing pair of stars. The inner and outer planets have radii 3.0 and 4.6 times that of the Earth, respectively. The binary star consists of a Sun-like star and a companion roughly one-third its size, orbiting each other every 7.45 days. With an orbital period of 49.5 days, eighteen transits of the inner planet have been observed, allowing a detailed characterization of its orbit and those of the stars. The outer planet's orbital period is 303.2 days, and although the planet is not Earth-like, it resides within the classical "habitable zone", where liquid water could exist on an Earth-like planet. With its two known planets, Kepler-47 establishes that close binary stars can host complete planetary systems.

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

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

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

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

GU Boo is one of only a relatively small number of well studied double-lined eclipsing binaries that contain low-mass stars. L'opez-Morales & Ribas (2005, hereafter LR05) present a comprehensive analysis of multi-color light and radial velocity curves for this system. The GU Boo light curves presented in LR05 had substantial asymmetries, which were attributed to large spots. In spite of the asymmetry LR05 derived masses and radii accurate to about 2%. We obtained additional photometry of GU Boo using both a CCD and a single-channel photometer and modeled the light curves with the ELC software to determine if the large spots in the light curves give rise to systematic errors at the few percent level. We also modeled the original light curves from LR05 using models with and without spots. We derived a radius of the primary of 0.6329 +- 0.0026 Solar-Radii, 0.6413 +- 0.0049 SR, and 0.6373 +- 0.0029 SR from the CCD, photoelectric, and LR05 data, respectively. Each of these measurements agrees with the value reported by LR05 (R1=0.623 +- 0.016 SR) at the level of about 2%. In addition, the spread in these values is about 1-2% from the mean. For the secondary, we derive radii of 0.6074 +- 0.0035 SR, 0.5944 +- 0.0069 SR, and 0.5976 +- 0.0059 SR from the three respective data sets. The LR05 value is R2=0.620 +- 0.020 SR, which is about 2-3% larger than each of the three values we found. The spread in these values is about 2% from the mean. The systematic difference between our three determinations of the secondary radius and that of LR05 might be attributed to differences in the modelling process and codes used. Our own fits suggest that, for GU Boo at least, using accurate spot modelling of a single set of multicolor light curves results in radii determinations accurate at the 2% level.

In this letter we present several new three and four-planet solutions based on the most current available radial velocity data for HD160691 ($\mu$Ara). These solutions are optimised using the Planetary Orbit Fitting Process (POFP) which is programmed and executed in MATLAB. POFP is based on a full integration of the system's multiple-body Newtonian equations of motion and on a multi level optimisation utilizing a variety of algorithms. The POFP solutions are presented in the context of the Keplerian-based solutions already appearing in the literature which we have reproduced here for comparison. The different solutions and their properties are compared over all data sets separately and combined. The new solutions do not seem to exhibit instabilities and are both co-planar and three-dimensional. We also provide a comparative prediction of the published and new solutions showing their diversion in the near future. In the short term, this projection will allow to choose between the variety of solutions as further observations are made.