With hydrodynamical simulations we examine the evolution of a protoplanetary disc around $\alpha$ Centauri B including the effect of the eccentric orbit binary companion $\alpha$ Centauri A. The initially circular orbit disc undergoes two types of eccentricity growth. First, the eccentricity oscillates on the orbital period of the binary, $P_{\rm orb}$, due to the eccentricity of the binary orbit. Secondly, for a sufficiently small disc aspect ratio, the disc undergoes global forced eccentricity oscillations on a time-scale of around $20\,P_{\rm orb}$. These oscillations damp out through viscous dissipation leaving a quasi-steady eccentricity profile for the disc that oscillates only on the binary orbital period. The time-averaged global eccentricity is in the range 0.05-0.1, with no precession in the steady state. The periastrons of the gas particles are aligned to one another. The higher the disc viscosity, the higher the disc eccentricity. With $N$-body simulations we examine the evolution of a disc of planetesimals that forms with the orbital properties of the quasi-steady protoplanetary disc. We find that the average magnitude of the eccentricity of particles increases and their periastrons become misaligned to each other once they decouple from the gas disc. The low planetesimal collision velocity required for planet formation suggests that for planet formation to have occurred in a disc of planetesimals formed from a protoplanetary disc around $\alpha$ Centauri B, said disc's viscosity must be have been small and planet formation must have occurred at orbital radii smaller than about $2.5\,\rm au$. Planet formation may be easier with the presence of gas.

An intriguing question in the context of dynamics arises: Could a moon possess a moon itself? Such a configuration does not exist in the Solar System, although this may be possible in theory. Kollmeier et al. (2019) determined the critical size of a satellite necessary to host a long-lived sub-satellite, or submoon. However, the orbital constraints for these submoons to exist are still undetermined. Domingos et al. (2006) indicated that moons are stable out to a fraction of the host planet Hill radius $R_{H,p}$, which in turn depends on the eccentricity of its host's orbit. Motivated by this, we simulate a system of exomoons and submoons for $10^5$ planetary orbits, while considering many initial orbital phases to obtain the critical semimajor axis in terms of $R_{H,p}$ or the hosts satellite's Hill radius $R_{H,sat}$, respectively. We find that, assuming circular coplanar orbits, the stability limit for exomoons is 0.40 $R_{H,p}$ and for a submoon is 0.33 $R_{H,sat}$. Additionally, we discuss the observational feasibility of detecting these sub-satellites through photometric, radial velocity, or direct imaging observations using the Neptunes-sized exomoon candidates Kepler 1625b-I (Teachey et al. 2018) and identify how stability can shape the identification of future candidates.

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.

Previous studies have shown that planets that rotate retrograde (backwards with respect to their orbital motion) generally experience less severe obliquity variations than those that rotate prograde (the same direction as their orbital motion). Here we examine retrograde-rotating planets on eccentric orbits and find a previously unknown secular spin-orbit resonance that can drive significant obliquity variations. This resonance occurs when the frequency of the planet's rotation axis precession becomes commensurate with an orbital eigenfrequency of the planetary system. The planet's eccentricity enables a participating orbital frequency through an interaction in which the apsidal precession of the planet's orbit causes a cyclic nutation of the planet's orbital angular momentum vector. The resulting orbital frequency follows the relationship $f = 2 \dot{\varpi} - \dot{\Omega}$, where $\dot{\varpi}$ and $\dot{\Omega}$ are the rates of the planet's changing longitude of periapsis and ascending node, respectively. We test this mechanism by simulating cases of a simple Earth-Jupiter system, and confirm the predicted resonance. Over the course of 100 Myr, the test Earths with rotation axis precession rates near the predicted resonant frequency experienced pronounced obliquity variations of order $10^\circ$-$30^\circ$. These variations can be significant, and suggest that while retrograde rotation is a stabilizing influence most of the time, retrograde rotators can experience large obliquity variations if they are on eccentric orbits and enter this spin-orbit resonance.

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.

Planets that orbit only one of the stars in stellar binary systems (i.e., circumstellar) are dynamically constrained to a limited range of orbital parameters and thus understanding conditions on their stability is of great importance in exoplanet searches. We perform $\sim$700 million N-body simulations to identify how stability regions depend on properties of the binary, as well as, the starting planetary inclination and mean longitude relative to the binary orbit. Moreover, we provide grid interpolation maps and lookup tables for the community to use our results. Through Monte-Carlo methods we determine that planets with a semimajor axis $a_p$ $\lesssim$8\% of the binary semimajor axis $a_{bin}$ will likely be stable given the known distribution of binary star parameters. This estimate varies in the Lidov-Kozai regime or for retrograde orbits to 4\% or 10\% of $a_{bin}$, respectively. Our method to quickly determine the circumstellar stability limit is important for interpreting observations of binaries using direct imaging with JWST, photometry with TESS, or even astrometry with Gaia.

Changes in planetary obliquity, or axial tilt, influence the climates on Earth-like planets. In the solar system, the Earth's obliquity is stabilized due to interactions with our moon and the resulting {small amplitude variations ($\sim$2.4\degree)} are beneficial for advanced life. Most Sun-like stars have at least one stellar companion and the habitability of circumstellar exoplanets is shaped by their stellar companion. We show that a stellar companion can dramatically change whether {Earth-like obliquity stability is} possible through planetary orbital precession relative to the binary orbit or resonant pumping of the obliquity through spin-orbit interactions. We present a new formalism for the planetary spin precession that accounts for orbital misalignments between the planet and binary. Using numerical modeling in $\alpha$ Centauri AB we show: a stark contrast between the planetary obliquity variations depending on the host star, planetary neighbors limit the possible spin states for {Earth-like obliquity stability}, and the presence of a moon can destabilize the obliquity, defying our Earth-based expectations. An Earth-like rotator orbiting the primary star would experience {small} obliquity variations for 87\%, 74\%, or 54\% of Solar type binaries, depending on the mass of the primary (0.8, 1.0, or 1.2 M$_\odot$, respectively). Thus, Earth-like planets likely experience much larger obliquity variations, with more extreme climates, unless they are in specific states, such as orbiting nearly planar with the binary and rotating retrograde (backwards) like Venus.

Alpha Centauri A is the closest solar-type star to the Sun and offers an excellent opportunity to detect the thermal emission of a mature planet heated by its host star. The MIRI coronagraph on JWST can search the 1-3 AU (1"-2") region around alpha Cen A which is predicted to be stable within the alpha Cen AB system. We demonstrate that with reasonable performance of the telescope and instrument, a 20 hr program combining on-target and reference star observations at 15.5 um could detect thermal emission from planets as small as ~5 RE. Multiple visits every 3-6 months would increase the geometrical completeness, provide astrometric confirmation of detected sources, and push the radius limit down to ~3 RE. An exozodiacal cloud only a few times brighter than our own should also be detectable, although a sufficiently bright cloud might obscure any planet present in the system. While current precision radial velocity (PRV) observations set a limit of 50-100 ME at 1-3 AU for planets orbiting alpha Cen A, there is a broad range of exoplanet radii up to 10 RE consistent with these mass limits. A carefully planned observing sequence along with state-of-the-art post-processing analysis could reject the light from alpha Cen A at the level of ~10^-5 at 1"-2" and minimize the influence of alpha Cen B located 7-8" away in the 2022-2023 timeframe. These space-based observations would complement on-going imaging experiments at shorter wavelengths as well as PRV and astrometric experiments to detect planets dynamically. Planetary demographics suggest that the likelihood of directly imaging a planet whose mass and orbit are consistent with present PRV limits is small, ~5%, and possibly lower if the presence of a binary companion further reduces occurrence rates. However, at a distance of just 1.34 pc, alpha Cen A is our closest sibling star and certainly merits close scrutiny.

Variations in the axial tilt, or obliquity, of terrestrial planets can affect their climates and therefore their habitability. Kepler-62f is a 1.4 R$_\oplus$ planet orbiting within the habitable zone of its K2 dwarf host star (Borucki et al. 2013). We perform N-body simulations that monitor the evolution of obliquity of Kepler-62f for 10 million year timescales to explore the effects on model assumptions, such as the masses of the Kepler-62 planets and the possibility of outer bodies. Significant obliquity variation occurs when the rotational precession frequency overlaps with one or more of the secular orbital frequencies, but most variations are limited to $\lesssim$10$^\circ$. Moderate variations ($\sim$10$^\circ - 20^\circ$) can occur over a broader range of initial obliquities when the relative nodal longitude ($\Delta\Omega$) overlaps with the frequency and phase of a given secular mode. However, we find that adding outer gas giants on long period orbits ($\gtrsim$ 1000 days) can produce large ($\sim$60$^\circ$) variations in obliquity if Kepler-62f has a very rapid (4 hr) rotation period. The possibility of giant planets on long period orbits impacts the climate and habitability of Kepler-62f through variations in the latitudinal surface flux, where the timescale for large variation can occur on million year timescales.

Kepler-1647b is the most recently discovered planet that transits two stars, i.e., a circumbinary planet (CBP). Due to its large orbital separation, Kepler-1647b stands out from the rest of the Kepler CBPs, which mostly reside on much tighter orbits near the stability limit. The large separation of Kepler-1647b challenges inward disk migration as a dominant formation pathway, suggested by the other Kepler CBPs. In this paper, we consider the possibility of an undetected planet near the stability limit by examining observational consequences of such a planet. We calculate the transit probability of the putative planet, transit timing variations (TTVs) of the known planet, and eclipsing timing variations (ETVs) of the host binary caused by the putative planet. We find the presence of a $\gtrsim$30M$_{\oplus}$ inner planet to be highly unlikely near the stability limit. In addition, we provide future TTV observation windows, which will further constrain possible undetected planets with lower masses.

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.

Modern studies of the early solar system routinely invoke the possibility of an orbital instability among the giant planets triggered by gravitational interactions between the planets and a massive exterior disk of planetesimals. Previous works have suggested that this instability can be substantially delayed ($\sim$100s Myr) after the formation of the giant planets. Bodies in the disk are typically treated in a semi-active manner, wherein their gravitational force on the planets is included, but interactions between the planetesimals are ignored. We perform $N$-body numerical simulations using GENGA, which makes use of GPUs to allow for the inclusion of all gravitational interactions between bodies. Although our simulated Kuiper belt particles are more massive than the probable masses of real primordial Kuiper belt objects, our simulations indicate that the self-stirring of the primordial Kuiper belt is very important to the dynamics of the giant planet instability. We find that interactions between planetesimals dynamically heat the disk and typically prevent the outer solar system instability from being delayed by more than a few tens of million years after giant planet formation. Longer delays occur in a small fraction of systems that have at least 3.5 AU gaps between the planets and planetesimal disk. Our final planetary configurations match the solar system at a rate consistent with other previous works in most regards. Pre-instability heating of the disk typically yields final Jovian eccentricities comparable to the modern solar system value, which has been a difficult constraint to match in past works.

The stability limit for circumbinary planets (CBPs) is not well defined and can depend on initial parameters defining either the planetary orbit or the inner binary orbit. We expand on the work of Holman & Wiegert (1999, AJ 117, 621) to develop numerical tools for quick, easy, and accurate determination of the stability limit. The results of our simulations, as well as our numerical tools, are available to the community through $\texttt{Zenodo}$ and $\texttt{GitHub}$, respectively. We employ a grid interpolation method based on $\sim$150 million full N-body simulations of initially circular, coplanar systems and compare to the 9 known Kepler CBP systems. Using a formalism from planet packing studies, we find that 55% of the Kepler CBP systems allow for an additional equal-mass planet to potentially exist on an interior orbit relative to the observed planet. Therefore, we do $\textit{not}$ find strong evidence for a pile-up in the Kepler CBP systems and more detections are needed to adequately characterize the formation mechanisms for the CBP population. Observations from the Transiting Exoplanet Survey Satellite are expected to substantially increase the number of detections using the unique geometry of CBP systems, where multiple transits can occur during a single conjunction.

We continue to investigate the binary system Kepler-16, consisting of a K-type main-sequence star, a red dwarf, and a circumbinary Saturnian planet. As part of our study, we describe the system's habitable zone based on different climate models. We also report on stability investigations for possible Earth-mass Trojans while expanding a previous study by B. L. Quarles and collaborators given in 2012. For the climate models we carefully consider the relevance of the system's parameters. Furthermore, we pursue new stability simulations for the Earth-mass objects starting along the orbit of Kepler-16b. The eccentricity distribution as obtained prefers values close to circular, whereas the inclination distribution remains flat. The stable solutions are distributed near the co-orbital Lagrangian points, thus enhancing the plausibility that Earth-mass Trojans might be able to exist in the Kepler-16(AB) system.

We perform long-term simulations, up to ten billion years, of closely-spaced configurations of 2 -- 6 planets, each as massive as the Earth, traveling on nested orbits about either stellar component in $\alpha$ Centauri AB. The innermost planet initially orbits at either the inner edge of its star's empirical habitable zone (HZ) or the inner edge of its star's conservative HZ. Although individual planets on low inclination, low eccentricity, orbits can survive throughout the habitable zones of both stars, perturbations from the companion star require that the minimum spacing of planets in multi-planet systems within the habitable zones of each star must be significantly larger than the spacing of similar multi-planet systems orbiting single stars in order to be long-lived. The binary companion induces a forced eccentricity upon the orbits of planets in orbit around either star. Planets on appropriately-phased circumstellar orbits with initial eccentricities equal to their forced eccentricities can survive on more closely spaced orbits than those with initially circular orbits, although the required spacing remains higher than for planets orbiting single stars. A total of up to nine planets on nested prograde orbits can survive for the current age of the system within the empirical HZs of the two stars, with five of these orbiting $\alpha$ Centauri B and four orbiting $\alpha$ Centauri A.

We extend our study of the extent of the regions within the $\alpha$ Centauri AB star system where small planets are able to orbit for billion-year timescales (Quarles & Lissauer 2016, AJ 151, 111) to investigate the effects of minimizing the forced eccentricity of initial trajectories. We find that initially prograde, circumstellar orbits require a piecewise quadratic function to accurately approximate forced eccentricity as a function of semimajor axis, but retrograde orbits can be modeled using a linear function. Circumbinary orbits in the $\alpha$ Centauri AB system are less affected by the forced eccentricity. {Planets on circumstellar orbits that begin with eccentricity vectors near their forced values are generally stable, up to $\sim$$10^9$ yr, out to a larger semimajor axis than are planets beginning on circular orbits. The amount by which the region of stability expands is much larger for retrograde orbits than it is for prograde orbits. The location of the stability boundary for two planet systems on prograde, circular orbits is much more sensitive to the initial eccentricity state than it is for analogous single planet systems.

TRAPPIST-1 is a nearby ultra-cool dwarf that is host to a remarkable planetary system consisting of seven transiting planets. The orbital properties and radii of the planets have been well-constrained, and recently the masses of the inner six planets have been measured with additional ground and space-based photometric observations. Large uncertainties in these mass measurements have prevented a robust analysis of the planetary compositions. Here we perform many thousands of N-body dynamical simulations with planet properties perturbed from the observed values and identify those that are stable for millions of years. This allows us to identify self-consistent orbital solutions that can be used in future studies. From our range of dynamical masses, we find that most of the planets are consistent with an Earth-like composition, where TRAPPIST-1f is likely to have a volatile-rich envelope.

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 evaluate the extent of the regions within the $\alpha$ Centauri AB star system where small planets are able to orbit for billion-year timescales, and we calculate the positions on the sky plane where planets on stable orbits about either stellar component may appear. We confirm the qualitative results of Wiegert and Holman (AJ 113, 1445, 1997) regarding the approximate size of the regions of stable orbits, which are larger for retrograde orbits relative to the binary than for prograde orbits. Additionally, we find that mean motion resonances with the binary orbit leave an imprint on the limits of orbital stability, and the effects of the Lidov-Kozai mechanism are also readily apparent.

We present the seventh Kepler planet candidate catalog, which is the first to be based on the entire, uniformly processed, 48 month Kepler dataset. This is the first fully automated catalog, employing robotic vetting procedures to uniformly evaluate every periodic signal detected by the Q1-Q17 Data Release 24 (DR24) Kepler pipeline. While we prioritize uniform vetting over the absolute correctness of individual objects, we find that our robotic vetting is overall comparable to, and in most cases is superior to, the human vetting procedures employed by past catalogs. This catalog is the first to utilize artificial transit injection to evaluate the performance of our vetting procedures and quantify potential biases, which are essential for accurate computation of planetary occurrence rates. With respect to the cumulative Kepler Object of Interest (KOI) catalog, we designate 1,478 new KOIs, of which 402 are dispositioned as planet candidates (PCs). Also, 237 KOIs dispositioned as false positives (FPs) in previous Kepler catalogs have their disposition changed to PC and 118 PCs have their disposition changed to FP. This brings the total number of known KOIs to 8,826 and PCs to 4,696. We compare the Q1-Q17 DR24 KOI catalog to previous KOI catalogs, as well as ancillary Kepler catalogs, finding good agreement between them. We highlight new PCs that are both potentially rocky and potentially in the habitable zone of their host stars, many of which orbit solar-type stars. This work represents significant progress in accurately determining the fraction of Earth-size planets in the habitable zone of Sun-like stars. The full catalog is publicly available at the NASA Exoplanet Archive.

Venus currently rotates slowly, with its spin controlled by solid-body and atmospheric thermal tides. However, conditions may have been far different 4 billion years ago, when the Sun was fainter and most of the carbon within Venus could have been in solid form, implying a low-mass atmosphere. We investigate how the obliquity would have varied for a hypothetical rapidly rotating Early Venus. The obliquity variation structure of an ensemble of hypothetical Early Venuses is simpler than that Earth would have if it lacked its large Moon (Lissauer et al., 2012), having just one primary chaotic regime at high prograde obliquities. We note an unexpected long-term variability of up to $\pm7^\circ$ for retrograde Venuses. Low-obliquity Venuses show very low total obliquity variability over billion-year timescales -- comparable to that of the real Moon-influenced Earth.

The primary Kepler Mission provided nearly continuous monitoring of ~200,000 objects with unprecedented photometric precision. We present the final catalog of eclipsing binary systems within the 105 square degree Kepler field of view. This release incorporates the full extent of the data from the primary mission (Q0-Q17 Data Release). As a result, new systems have been added, additional false positives have been removed, ephemerides and principal parameters have been recomputed, classifications have been revised to rely on analytical models, and eclipse timing variations have been computed for each system. We identify several classes of systems including those that exhibit tertiary eclipse events, systems that show clear evidence of additional bodies, heartbeat systems, systems with changing eclipse depths, and systems exhibiting only one eclipse event over the duration of the mission. We have updated the period and galactic latitude distribution diagrams and included a catalog completeness evaluation. The total number of identified eclipsing and ellipsoidal binary systems in the Kepler field of view has increased to 2878, 1.3% of all observed Kepler targets. An online version of this catalog with downloadable content and visualization tools is maintained at http://keplerEBs.villanova.edu.

The three-body problem, which describes three masses interacting through Newtonian gravity without any restrictions imposed on the initial positions and velocities of these masses, has attracted the attention of many scientists for more than 300 years. In this paper, we present a review of the three-body problem in the context of both historical and modern developments. We describe the general and restricted (circular and elliptic) three-body problems, different analytical and numerical methods of finding solutions, methods for performing stability analysis, search for periodic orbits and resonances, and application of the results to some interesting astronomical and space dynamical settings. We also provide a brief presentation of the general and restricted relativistic three-body problems, and discuss their astronomical applications.

\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.

The Kepler mission discovered 2842 exoplanet candidates with 2 years of data. We provide updates to the Kepler planet candidate sample based upon 3 years (Q1-Q12) of data. Through a series of tests to exclude false-positives, primarily caused by eclipsing binary stars and instrumental systematics, 855 additional planetary candidates have been discovered, bringing the total number known to 3697. We provide revised transit parameters and accompanying posterior distributions based on a Markov Chain Monte Carlo algorithm for the cumulative catalogue of Kepler Objects of Interest. There are now 130 candidates in the cumulative catalogue that receive less than twice the flux the Earth receives and more than 1100 have a radius less than 1.5 Rearth. There are now a dozen candidates meeting both criteria, roughly doubling the number of candidate Earth analogs. A majority of planetary candidates have a high probability of being bonafide planets, however, there are populations of likely false-positives. We discuss and suggest additional cuts that can be easily applied to the catalogue to produce a set of planetary candidates with good fidelity. The full catalogue is publicly available at the NASA Exoplanet Archive.

We present integrations of a model Solar System with five terrestrial planets (beginning ~30-50 Myr after the formation of primitive Solar System bodies) in order to determine the preferred regions of parameter space leading to a giant impact that resulted in the formation of the Moon. Our results indicate which choices of semimajor axes and eccentricities for Theia (the proto-Moon) at this epoch can produce a late Giant Impact, assuming that Mercury, Venus, and Mars are near the current orbits. We find that the likely semimajor axis of Theia, at the epoch when our simulations begin, depends on the assumed mass ratio of Earth-Moon progenitors (8/1, 4/1, or 1/1). The low eccentricities of the terrestrial planets are most commonly produced when the progenitors have similar semimajor axes at the epoch when our integrations commence. Additionally, we show that mean motion resonances among the terrestrial planets and perturbations from the giant planets can affect the dynamical evolution of the system leading to a late Giant Impact.

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.

Over 2500 eclipsing binaries were identified and characterized from the ultra-precise photometric data provided by the Kepler space telescope. Kepler is now beginning its second mission, K2, which is proving to again provide ultra-precise photometry for a large sample of eclipsing binary stars. In the 1951 light curves covering 12 days in the K2 engineering data-set, we have identified and determined the ephemerides for 31 eclipsing binaries that demonstrate the capabilities for eclipsing binary science in the upcoming campaigns in K2. Of those, 20 are new discoveries. We describe both manual and automated approaches to harvesting the complete set of eclipsing binaries in the K2 data, provide identifications and details for the full set of eclipsing binaries present in the engineering data-set, and discuss the prospects for application of eclipsing binary searches in the K2 mission.

Depending on the planetary orbit around the host star(s), a planet could orbit either one or both stars in a binary system as S-type or P-type, respectively. We have analysed the dynamics of the S-type planetary system in HD 196885 AB with an emphasis on a planet with a higher orbital inclination relative to the binary plane. The mean exponential growth factor of nearby orbits (MEGNO) maps are used as an indicator to determine regions of periodicity and chaos for the various choices of the planet's semimajor axis, eccentricity and inclination with respect to the previously determined observational uncertainties. We have quantitatively mapped out the chaotic and quasi-periodic regions of the system's phase space which indicate a likely regime of the planet's inclination. In addition, we inspect the resonant angle to determine whether alternation between libration and circulation occurs as a consequence of Kozai oscillations, a probable mechanism that can drive the planetary orbit to a very large inclination. Also, we demonstrate the possible higher mass limit of the planet and improve upon the current dynamical model based on our analysis.

The quest for Earth-like planets represents a major focus of current exoplanet research. While planets that are Earth-sized and smaller have been detected, these planets reside in orbits that are too close to their host star to allow liquid water on their surface. We present the detection of Kepler-186f, a 1.11+\-0.14 Earth radius planet that is the outermost of five planets - all roughly Earth-sized - that transit a 0.47+\-0.05 Rsun star. The intensity and spectrum of the star's radiation places Kepler-186f in the stellar habitable zone, implying that if Kepler-186f has an Earth-like atmosphere and H2O at its surface, then some of this H2O is likely to be in liquid form.

The orbits of two individual planets in two known binary star systems, \gamma Cephei and HD 196885 are numerically integrated using various numerical techniques to assess the chaotic or quasi-periodic nature of the dynamical system considered. The Hill stability (HS) function which measures the orbital perturbation of a planet around the primary star due to the secondary star is calculated for each system. The maximum Lyapunov exponent (MLE) time series are generated to measure the divergence/convergence rate of stable manifolds, which are used to differentiate between chaotic and non-chaotic orbits. Then, we calculate dynamical Mean Exponential Growth factor of Nearby Orbits (MEGNO) maps from solving the variational equations along with the equations of motion. These maps allow us to accurately differentiate between stable and unstable dynamical systems. The results obtained from the analysis of HS, MLE, and MEGNO maps are analysed for their dynamical variations and resemblance. The HS test for the planets shows stability and quasi-periodicity for at least ten million years. The MLE and the MEGNO maps have also indicated the local quasi-periodicity and global stability in a relatively short integration period. The orbital stability of the systems is tested using each indicator for various values of planet inclinations (i_{pl} \le 25^\circ) and binary eccentricities. The reliability of HS criterion is also discussed based on its stability results compared with the MLE and MEGNO maps.

The aim of our study is to investigate the possibility of habitable moons orbiting the giant planet HD 23079b, a Jupiter-mass planet, which follows a low-eccentricity orbit in the outer region of HD 23079's habitable zone. We show that HD 23079b is able to host habitable moons in prograde and retrograde orbits, as expected, noting that the outer stability limit for retrograde orbits is increased by nearly 90% compared to that of prograde orbits, a result consistent with previous generalized studies. For the targeted parameter space it was found that the outer stability limit for habitable moons varies between 0.05236 and 0.06955 AU (prograde orbits) and between 0.1023 and 0.1190 AU (retrograde orbits) depending on the orbital parameters of the Jupiter-type planet if a minimum mass is assumed. These intervals correspond to 0.306 and 0.345 (prograde orbits) and 0.583 and 0.611 (retrograde orbits) of the planet's Hill radius. Larger stability limits are obtained if an increased value for the planetary mass m_p is considered; they are consistent with the theoretically deduced relationship of m_p^{1/3}. Finally, we compare our results to the statistical formulae of Domingos et al. (2006) [MNRAS 373, 1227], indicating both concurrence and limitations.

Our aim is to identify and classify mean-motion resonances (MMRs) for the coplanar circular restricted three-body problem (CR3BP) for mass ratios between 0.10 and 0.50. Our methods include the maximum Lyapunov exponent, which is used as an indicator for the location of the resonances, the Fast Fourier Transform (FFT) used for determining what kind of resonances are present, and the inspection of the orbital elements to classify the periodicity. We show that the 2:1 resonance occurs the most frequently. Among other resonances, the 3:1 resonance is the second most common, and furthermore both 3:2 and 5:3 resonances occur more often than the 4:1 resonance. Moreover, the resonances in the coplanar CR3BP are classified based on the behaviour of the orbits. We show that orbital stability is ensured for high values of resonance (i.e., high ratios) where only a single resonance is present. The resonances attained are consistent with the previously established resonances for the Solar System, i.e., specifically, in regards to the asteroid belt. Previous work employed digital filtering and Lyapunov characteristic exponents to determine stochasticity of the eccentricity, which is found to be consistent with our usage of Lyapunov exponents as an alternate approach based on varying the mass ratio instead of the eccentricity. Our results are expected to be of principal interest to future studies, including augmentations to observed or proposed resonances, of extra-solar planets in binary stellar systems.

We demonstrate that habitable Earth-mass planets and moons can exist in the Kepler-16 system, known to host a Saturn-mass planet around a stellar binary, by investigating their orbital stability in the standard and extended habitable zone (HZ). We find that Earth-mass planets in satellite-like (S-type)orbits are possible within the standard HZ in direct vicinity of Kepler-16b, thus constituting habitable exomoons. However, Earth-mass planets cannot exist in planetary-like (P-type) orbits around the two stellar components within the standard HZ. Yet, P-type Earth-mass planets can exist superior to the Saturnian planet in the extended HZ pertaining to considerably enhanced back-warming in the planetary atmosphere if facilitated. We briefly discuss the potential detectability of such habitable Earth-mass moons and planets positioned in satellite and planetary orbits, respectively. The range of inferior and superior P-type orbits in the HZ is between 0.657 to 0.71 AU and 0.95 to 1.02 AU, respectively.

We provide a detailed theoretical study aimed at the observational finding about the nu Octantis binary system that indicates the possible existence of a Jupiter-type planet in this system. If a prograde planetary orbit is assumed, it has earlier been argued that the planet, if existing, should be located outside the zone of orbital stability. However, a previous study by Eberle & Cuntz (2010) [ApJ 721, L168] concludes that the planet is most likely stable if assumed to be in a retrograde orbit with respect to the secondary system component. In the present work, we significantly augment this study by taking into account the observationally deduced uncertainty ranges of the orbital parameters for the stellar components and the suggested planet. Furthermore, our study employs additional mathematical methods, which include monitoring the Jacobi constant, the zero velocity function, and the maximum Lyapunov exponent. We again find that the suggested planet is indeed possible if assumed to be in a retrograde orbit, but it is virtually impossible if assumed in a prograde orbit. Its existence is found to be consistent with the deduced system parameters of the binary components and of the suggested planet, including the associated uncertainty bars given by observations.

We establish a criterion for the stability of planetary orbits in stellar binary systems by using Lyapunov exponents and power spectra for the special case of the circular restricted 3-body problem (CR3BP). The centerpiece of our method is the concept of Lyapunov exponents, which are incorporated into the analysis of orbital stability by integrating the Jacobian of the CR3BP and orthogonalizing the tangent vectors via a well-established algorithm originally developed by Wolf et al. The criterion for orbital stability based on the Lyapunov exponents is independently verified by using power spectra. The obtained results are compared to results presented in the two previous papers of this series. It is shown that the maximum Lyapunov exponent can be used as an indicator for chaotic behaviour of planetary orbits, which is consistent with previous applications of this method, particularly studies for the Solar System. The chaotic behaviour corresponds to either orbital stability or instability, and it depends solely on the mass ratio of the binary components and the initial distance ratio of the planet relative to the stellar separation distance. Our theoretical results allow us to link the study of planetary orbital stability to chaos theory noting that there is a large array of literature on the properties and significance of Lyapunov exponents. Although our results are given for the special case of the CR3BP, we expect that it may be possible to augment the proposed Lyapunov exponent criterion to studies of planets in generalized stellar binary systems, which is strongly motivated by existing observational results as well as results expected from ongoing and future planet search missions.

We investigate the possibility of habitable Trojan planets in the HD 23079 star-planet system. This system consists of a solar-type star and a Jupiter-type planet, which orbits the star near the outer edge of the stellar habitable zone in an orbit of low eccentricity. We find that in agreement with previous studies Earth-mass habitable Trojan planets are possible in this system, although the success of staying within the zone of habitability is significantly affected by the orbital parameters of the giant planet and by the initial condition of the theoretical Earth-mass planet. In one of our simulations, the Earth-mass planet is captured by the giant planet and thus becomes a habitable moon.