We report the detection of the first circumbinary planet found by TESS. The target, a known eclipsing binary, was observed in sectors 1 through 12 at 30-minute cadence and in sectors 4 through 12 at two-minute cadence. It consists of two stars with masses of 1.1 MSun and 0.3 MSun on a slightly eccentric (0.16), 14.6-day orbit, producing prominent primary eclipses and shallow secondary eclipses. The planet has a radius of ~6.9 REarth and was observed to make three transits across the primary star of roughly equal depths (~0.2%) but different durations -- a common signature of transiting circumbinary planets. Its orbit is nearly circular (e ~ 0.09) with an orbital period of 95.2 days. The orbital planes of the binary and the planet are aligned to within ~1 degree. To obtain a complete solution for the system, we combined the TESS photometry with existing ground-based radial-velocity observations in a numerical photometric-dynamical model. The system demonstrates the discovery potential of TESS for circumbinary planets, and provides further understanding of the formation and evolution of planets orbiting close binary stars.

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.

The Kepler-11 planetary system contains six transiting planets ranging in size from 1.8 to 4.2 times the radius of Earth. Five of these planets orbit in a tightly-packed configuration with periods between 10 and 47 days. We perform a dynamical analysis of the system based upon transit timing variations observed in more than three years of \ik photometric data. Stellar parameters are derived using a combination of spectral classification and constraints on the star's density derived from transit profiles together with planetary eccentricity vectors provided by our dynamical study. Combining masses of the planets relative to the star from our dynamical study and radii of the planets relative to the star from transit depths together with deduced stellar properties yields measurements of the radii of all six planets, masses of the five inner planets, and an upper bound to the mass of the outermost planet, whose orbital period is 118 days. We find mass-radius combinations for all six planets that imply that substantial fractions of their volumes are occupied by constituents that are less dense than rock. The Kepler-11 system contains the lowest mass exoplanets for which both mass and radius have been measured.

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

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.

When an extrasolar planet passes in front of its star (transits), its radius can be measured from the decrease in starlight and its orbital period from the time between transits. Multiple planets transiting the same star reveal more: period ratios determine stability and dynamics, mutual gravitational interactions reflect planet masses and orbital shapes, and the fraction of transiting planets observed as multiples has implications for the planarity of planetary systems. But few stars have more than one known transiting planet, and none has more than three. Here we report Kepler spacecraft observations of a single Sun-like star that reveal six transiting planets, five with orbital periods between 10 and 47 days plus a sixth one with a longer period. The five inner planets are among the smallest whose masses and sizes have both been measured, and these measurements imply substantial envelopes of light gases. The degree of coplanarity and proximity of the planetary orbits imply energy dissipation near the end of planet formation.

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.