A dozen short-period detached binaries are known to host transiting circumbinary planets. In all circumbinary systems so far, the planetary and binary orbits are aligned within a couple of degrees. However, the obliquity of the primary star, which is an important tracer of their formation, evolution, and tidal history, has only been measured in one circumbinary system until now. EBLM J0608-59/TOI-1338 is a low-mass eclipsing binary system with a recently discovered circumbinary planet identified by TESS. Here, we perform high-resolution spectroscopy during primary eclipse to measure the projected stellar obliquity of the primary component. The obliquity is low, and thus the primary star is aligned with the binary and planetary orbits with a projected spin-orbit angle $\beta = 2.8 \pm 17.1$ deg. The rotation period of $18.1 \pm 1.6$ days implied by our measurement of $v\sin{i_\star}$ suggests that the primary has not yet pseudo-synchronized with the binary orbit, but is consistent with gyrochronology and weak tidal interaction with the binary companion. Our result, combined with the known coplanarity of the binary and planet orbits, is suggestive of formation from a single disc. Finally, we considered whether the spectrum of the faint secondary star could affect our measurements. We show through simulations that the effect is negligible for our system, but can lead to strong biases in $v\sin{i_\star}$ and $\beta$ for higher flux ratios. We encourage future studies in eclipse spectroscopy test the assumption of a dark secondary for flux ratios $\gtrsim 1$ ppt.

We present the discovery and validation of a three-planet system orbiting the nearby (31.1 pc) M2 dwarf star TOI-700 (TIC 150428135). TOI-700 lies in the TESS continuous viewing zone in the Southern Ecliptic Hemisphere; observations spanning 11 sectors reveal three planets with radii ranging from 1 R$_\oplus$ to 2.6 R$_\oplus$ and orbital periods ranging from 9.98 to 37.43 days. Ground-based follow-up combined with diagnostic vetting and validation tests enable us to rule out common astrophysical false-positive scenarios and validate the system of planets. The outermost planet, TOI-700 d, has a radius of $1.19\pm0.11$ R$_\oplus$ and resides in the conservative habitable zone of its host star, where it receives a flux from its star that is approximately 86% of the Earth's insolation. In contrast to some other low-mass stars that host Earth-sized planets in their habitable zones, TOI-700 exhibits low levels of stellar activity, presenting a valuable opportunity to study potentially-rocky planets over a wide range of conditions affecting atmospheric escape. While atmospheric characterization of TOI-700 d with the James Webb Space Telescope (JWST) will be challenging, the larger sub-Neptune, TOI-700 c (R = 2.63 R$_\oplus$), will be an excellent target for JWST and beyond. TESS is scheduled to return to the Southern Hemisphere and observe TOI-700 for an additional 11 sectors in its extended mission, which should provide further constraints on the known planet parameters and searches for additional planets and transit timing variations in the system.

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.

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.

Planets that revolve around a binary pair of stars are known as circumbinary planets. The orbital motion of the stars around their center of mass causes a periodic variation in the total instellation incident upon a circumbinary planet. This study uses both an analytic and numerical energy balance model to calculate the extent to which this effect can drive changes in surface temperature on circumbinary terrestrial planets. We show that the amplitude of the temperature variation is largely constrained by the effective heat capacity, which corresponds to the ocean-to-land ratio on the planet. Planets with large ocean fractions should experience only modest warming and cooling of only a few degrees, which suggests that habitability cannot be precluded for such circumbinary planets. Planets with large land fractions that experience extreme periodic forcing could be prone to changes in temperature of tens of degrees or more, which could drive more extreme climate changes that inhibit continuously habitable conditions.

We use a one-dimensional (1-D) cloud-free climate model to estimate habitable zone (HZ) boundaries for terrestrial planets of masses 0.1 M$_{E}$ and 5 M$_{E}$ around circumbinary stars of various spectral type combinations. Specifically, we consider binary systems with host spectral types F-F, F-G, F-K, F-M, G-G, G-K, G-M, K-K, K-M and M-M. Scaling the background N2 atmospheric pressure with the radius of the planet, we find that the inner edge of the HZ moves inwards towards the star for 5ME compared to 0.1ME planets for all spectral types. This is because the water-vapor column depth is smaller for larger planets and higher temperatures are needed before water vapor completely dominates the outgoing longwave radiation. The outer edge of the HZ changes little due to competing effects of the albedo and greenhouse effect. While these results are broadly consistent with the trend of single star HZ results for different mass planets, there are significant differences between single star and binary star systems for the inner edge of the HZ. Interesting combinations of stellar pairs from our 1-D model results can be used to explore for in-depth climate studies with 3-D climate models. We identify a common HZ stellar flux domain for all circumbinary spectral types

The Kepler, K2, and Transiting Exoplanet Survey Satellite (TESS) missions have provided a wealth of confirmed exoplanets, benefiting from a huge effort from the planet-hunting and follow-up community. With careful systematics mitigation, these missions provide precise photometric time series, which enable detection of transiting exoplanet signals. However, exoplanet hunting can be confounded by several factors, including instrumental noise, search biases, and host star variability. In this Letter, we discuss strategies to overcome these challenges using newly emerging techniques and tools. We demonstrate the power of new, fast open-source community tools (e.g., lightkurve, starry, celerite, exoplanet), and discuss four high signal-to-noise ratio (S/N) exoplanets that showcase specific challenges present in planet detection: K2-43c, K2-168c, K2-198c, and K2-198d. These planets have been undetected in several large K2 planet searches, despite having transit signals with S/N > 10. Two of the planets discussed here are new discoveries. In this work we confirm all four as true planets. Alongside these planet systems, we discuss three key challenges in finding small transiting exoplanets. The aim of this Letter is to help new researchers understand where planet detection efficiency gains can be made, and to encourage the continued use of K2 archive data. The considerations presented in this Letter are equally applicable to Kepler, K2, and TESS, and the tools discussed here are available for the community to apply to improve exoplanet discovery and fitting.

Planets occur most frequently around cool dwarfs, but only a handful of specific examples are known to orbit the latest-type M stars. Using TESS photometry, we report the discovery of two planets transiting the low-mass star called LP 791-18 (identified by TESS as TOI 736). This star has spectral type M6V, effective temperature 2960 K, and radius 0.17 R_Sun, making it the third-coolest star known to host planets. The two planets straddle the radius gap seen for smaller exoplanets; they include a 1.1 R_Earth planet on a 0.95 day orbit and a 2.3 R_Earth planet on a 5 day orbit. Because the host star is small the loss of light during these planets' transits is fairly large (0.4% and 1.7%). This has allowed us to detect both planets' transits from ground-based photometry, refining their radii and orbital ephemerides. In the future, radial velocity observations and transmission spectroscopy can both probe these planets' bulk interior and atmospheric compositions, and additional photometric monitoring would be sensitive to even smaller transiting planets.

We report the Transiting Exoplanet Survey Satellite (TESS) discovery of three terrestrial-sized planets transiting L 98-59 (TOI-175, TIC 307210830) -- a bright M dwarf at a distance of 10.6 pc. Using the Gaia-measured distance and broad-band photometry we find that the host star is an M3 dwarf. Combined with the TESS transits from three sectors, the corresponding stellar parameters yield planet radii ranging from 0.8REarth to 1.6REarth. All three planets have short orbital periods, ranging from 2.25 to 7.45 days with the outer pair just wide of a 2:1 period resonance. Diagnostic tests produced by the TESS Data Validation Report and the vetting package DAVE rule out common false positive sources. These analyses, along with dedicated follow-up and the multiplicity of the system, lend confidence that the observed signals are caused by planets transiting L 98-59 and are not associated with other sources in the field. The L 98-59 system is interesting for a number of reasons: the host star is bright (V = 11.7 mag, K = 7.1 mag) and the planets are prime targets for further follow-up observations including precision radial-velocity mass measurements and future transit spectroscopy with the James Webb Space Telescope; the near resonant configuration makes the system a laboratory to study planetary system dynamical evolution; and three planets of relatively similar size in the same system present an opportunity to study terrestrial planets where other variables (age, metallicity, etc.) can be held constant. L 98-59 will be observed in 4 more TESS sectors, which will provide a wealth of information on the three currently known planets and have the potential to reveal additional planets in the system.

We have adapted the algorithmic tools developed during the Kepler mission to vet the quality of transit-like signals for use on the K2 mission data. Using the four sets of publicly-available lightcurves on MAST, we produced a uniformly-vetted catalog of 772 transiting planet candidates from K2 as listed at the NASA Exoplanet archive in the K2 Table of Candidates. Our analysis marks 676 of these as planet candidates and 96 as false positives. All confirmed planets pass our vetting tests. 60 of our false positives are new identifications -- effectively doubling the overall number of astrophysical signals mimicking planetary transits in K2 data. Most of the targets listed as false positives in our catalog either show prominent secondary eclipses, transit depths suggesting a stellar companion instead of a planet, or significant photocenter shifts during transit. We packaged our tools into the open-source, automated vetting pipeline DAVE (Discovery and Vetting of Exoplanets) designed to streamline follow-up efforts by reducing the time and resources wasted observing targets that are likely false positives. DAVE will also be a valuable tool for analyzing planet candidates from NASA's TESS mission, where several guest-investigator programs will provide independent lightcurve sets -- and likely many more from the community. We are currently testing DAVE on recently-released TESS planet candidates and will present our results in a follow-up paper.

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.

Inspired by the recent Kepler discoveries of circumbinary planets orbiting nine close binary stars, we explore the fate of the former as the latter evolve off the main sequence. We combine binary star evolution models with dynamical simulations to study the orbital evolution of these planets as their hosts undergo common-envelope stages, losing in the process a tremendous amount of mass on dynamical timescales. Five of the systems experience at least one Roche-lobe overflow and common-envelope stages (Kepler-1647 experiences three), and the binary stars either shrink to very short orbits or coalesce; two systems trigger a double-degenerate supernova explosion. Kepler's circumbinary planets predominantly remain gravitationally bound at the end of the common-envelope phase, migrate to larger orbits, and may gain significant eccentricity; their orbital expansion can be more than an order of magnitude and can occur over the course of a single planetary orbit. The orbits these planets can reach are qualitatively consistent with those of the currently known post-common-envelope, eclipse-time variations circumbinary candidates. Our results also show that circumbinary planets can experience both modes of orbital expansion (adiabatic and non-adiabatic) if their host binaries undergo more than one common-envelope stage; multiplanet circumbinary systems like Kepler-47 can experience both modes during the same common-envelope stage. Additionally, unlike Mercury orbiting the Sun, a circumbinary planet with the same semi-major axis can survive the common envelope evolution of a close binary star with a total mass of 1 MSun.

In this brief communication we provide the rationale for, and the outcome of the International Astronomical Union (IAU) resolution vote at the XXIX-th General Assembly in Honolulu, Hawaii, in 2015, on recommended nominal conversion constants for selected solar and planetary properties. The problem addressed by the resolution is a lack of established conversion constants between solar and planetary values and SI units: a missing standard has caused a proliferation of solar values (e.g., solar radius, solar irradiance, solar luminosity, solar effective temperature and solar mass parameter) in the literature, with cited solar values typically based on best estimates at the time of paper writing. As precision of observations increases, a set of consistent values becomes increasingly important. To address this, an IAU Working Group on Nominal Units for Stellar and Planetary Astronomy formed in 2011, uniting experts from the solar, stellar, planetary, exoplanetary and fundamental astronomy, as well as from general standards fields to converge on optimal values for nominal conversion constants. The effort resulted in the IAU 2015 Resolution B3, passed at the IAU General Assembly by a large majority. The resolution recommends the use of nominal solar and planetary values, which are by definition exact and are expressed in SI units. These nominal values should be understood as conversion factors only, not as the true solar/planetary properties or current best estimates. Authors and journal editors are urged to join in using the standard values set forth by this resolution in future work and publications to help minimize further confusion.

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.

The Kepler Mission has provided unprecedented, nearly continuous photometric data of $\sim$200,000 objects in the $\sim$105 deg$^{2}$ field of view from the beginning of science operations in May of 2009 until the loss of the second reaction wheel in May of 2013. The Kepler Eclipsing Binary Catalog contains information including but not limited to ephemerides, stellar parameters and analytical approximation fits for every known eclipsing binary system in the Kepler Field of View. Using Target Pixel level data collected from Kepler in conjunction with the Kepler Eclipsing Binary Catalog, we identify false positives among eclipsing binaries, i.e. targets that are not eclipsing binaries themselves, but are instead contaminated by eclipsing binary sources nearby on the sky and show eclipsing binary signatures in their light curves. We present methods for identifying these false positives and for extracting new light curves for the true source of the observed binary signal. For each source, we extract three separate light curves for each quarter of available data by optimizing the signal-to-noise ratio, the relative percent eclipse depth and the flux eclipse depth. We present 289 new eclipsing binaries in the Kepler Field of View that were not targets for observation, and these have been added to the Catalog. An online version of this Catalog with downloadable content and visualization tools is maintained at http://keplerEBs.villanova.edu.

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 XXIXth IAU General Assembly in Honolulu adopted IAU 2015 Resolution B2 on recommended zero points for the absolute and apparent bolometric magnitude scales. The resolution was proposed by the IAU Inter-Division A-G Working Group on Nominal Units for Stellar and Planetary Astronomy after consulting with a broad spectrum of researchers from the astronomical community. Resolution B2 resolves the long-standing absence of an internationally-adopted zero point for the absolute and apparent bolometric magnitude scales. Resolution B2 defines the zero point of the absolute bolometric magnitude scale such that a radiation source with $M_{\rm Bol}$ = 0 has luminosity L$_{\circ}$ = 3.0128e28 W. The zero point of the apparent bolometric magnitude scale ($m_{\rm Bol}$ = 0) corresponds to irradiance $f_{\circ}$ = 2.518021002e-8 W/m$^2$. The zero points were chosen so that the nominal solar luminosity (3.828e26 W) adopted by IAU 2015 Resolution B3 corresponds approximately to $M_{\rm Bol}$(Sun) = 4.74, the value most commonly adopted in recent literature. The nominal total solar irradiance (1361 W/m$^2$) adopted in IAU 2015 Resolution B3 corresponds approximately to apparent bolometric magnitude $m_{\rm bol}$(Sun) = -26.832. Implicit in the IAU 2015 Resolution B2 definition of the apparent bolometric magnitude scale is an exact definition for the parsec (648000/$\pi$ au) based on the IAU 2012 Resolution B2 definition of the astronomical unit.

Astronomers commonly quote the properties of celestial objects in units of parameters for the Sun, Jupiter, or the Earth. The resolution presented here was proposed by the IAU Inter-Division Working Group on Nominal Units for Stellar and Planetary Astronomy and passed by the XXIXth IAU General Assembly in Honolulu. IAU 2015 Resolution B3 adopts a set of nominal solar, terrestrial, and jovian conversion constants for stellar and (exo)planetary astronomy which are defined to be exact SI values. While the nominal constants are based on current best estimates (CBEs; which have uncertainties, are not secularly constant, and are updated regularly using new observations), they should be interpreted as standard values and not as CBEs. IAU 2015 Resolution B3 adopts five solar conversion constants (nominal solar radius, nominal total solar irradiance, nominal solar luminosity, nominal solar effective temperature, and nominal solar mass parameter) and six planetary conversion constants (nominal terrestrial equatorial radius, nominal terrestrial polar radius, nominal jovian equatorial radius, nominal jovian polar radius, nominal terrestrial mass parameter, and nominal jovian mass parameter).

We report rotational periods for 16 young brown dwarfs in the nearby Upper Scorpius association, based on 72 days of high-cadence, high-precision photometry from the Kepler space telescope's K2 mission. The periods range from a few hours to two days (plus one outlier at 5 days), with a median just above one day, confirming that brown dwarfs, except at the very youngest ages, are fast rotators. Interestingly, four of the slowest rotators in our sample exhibit mid-infrared excess emission from disks; at least two also show signs of disk eclipses and accretion in the lightcurves. Comparing these new periods with those for two other young clusters and simple angular momentum evolution tracks, we find little or no rotational braking in brown dwarfs between 1-10 Myr, in contrast to low-mass stars. Our findings show that disk braking, while still at work, is inefficient in the substellar regime, thus provide an important constraint on the mass dependence of the braking mechanism.

We have studied the possibility that a third circumbinary planet in the Kepler-47 planetary system be the source of the single unexplained transiting event reported during the discovery of these planets. We applied the MEGNO technique to identify regions in the phase space where a third planet can maintain quasi-periodic orbits, and assessed the long-term stability of the three-planet system by integrating the entire 5 bodies (binary + planets) for 10 Myr. We identified several stable regions between the two known planets as well as a region beyond the orbit of Kepler-47c where the orbit of the third planet could be stable. To constrain the orbit of this planet, we used the measured duration of the unexplained transit event ($\sim 4.15$ hours) and compared that with the transit duration of the third planet in an ensemble of stable orbits. To remove the degeneracy among the orbits with similar transit durations, we considered the planet to be in a circular orbit and calculated its period analytically. The latter places an upper limit of 424 days on the orbital period of the third planet. Our analysis suggests that if the unexplained transit event detected during the discovery of the Kepler-47 circumbinary system is due to a planetary object, this planet will be in a low eccentricity orbit with a semimajor axis smaller than 1.24 AU. Further constraining of the mass and orbital elements of this planet requires a re-analysis of the entire currently available data including those obtained post-announcement of the discovery of this system. We present details of our methodology and discuss the implication of the results.

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 report the discovery of a transiting, Rp = 4.347+/-0.099REarth, circumbinary planet (CBP) orbiting the Kepler K+M Eclipsing Binary (EB) system KIC 12351927 (Kepler-413) every ~66 days on an eccentric orbit with ap = 0.355+/-0.002AU, ep = 0.118+/-0.002. The two stars, with MA = 0.820+/-0.015MSun, RA = 0.776+/-0.009RSun and MB = 0.542+/-0.008MSun, RB = 0.484+/-0.024RSun respectively revolve around each other every 10.11615+/-0.00001 days on a nearly circular (eEB = 0.037+/-0.002) orbit. The orbital plane of the EB is slightly inclined to the line of sight (iEB = 87.33+/-0.06 degrees) while that of the planet is inclined by ~2.5 degrees to the binary plane at the reference epoch. Orbital precession with a period of ~11 years causes the inclination of the latter to the sky plane to continuously change. As a result, the planet often fails to transit the primary star at inferior conjunction, causing stretches of hundreds of days with no transits (corresponding to multiple planetary orbital periods). We predict that the next transit will not occur until 2020. The orbital configuration of the system places the planet slightly closer to its host stars than the inner edge of the extended habitable zone. Additionally, the orbital configuration of the system is such that the CBP may experience Cassini-States dynamics under the influence of the EB, in which the planet's obliquity precesses with a rate comparable to its orbital precession. Depending on the angular precession frequency of the CBP, it could potentially undergo obliquity fluctuations of dozens of degrees (and complex seasonal cycles) on precession timescales.

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.

We report the discovery of a transiting, gas giant circumbinary planet orbiting the eclipsing binary KIC 4862625 and describe our independent discovery of the two transiting planets orbiting Kepler-47 (Orosz et al. 2012). We describe a simple and semi-automated procedure for identifying individual transits in light curves and present our follow-up measurements of the two circumbinary systems. For the KIC 4862625 system, the 0.52+/-0.018 RJup radius planet revolves every ~138 days and occults the 1.47+/-0.08 MSun, 1.7 +/-0.06 RSun F8 IV primary star producing aperiodic transits of variable durations commensurate with the configuration of the eclipsing binary star. Our best-fit model indicates the orbit has a semi-major axis of 0.64 AU and is slightly eccentric, e=0.1. For the Kepler-47 system, we confirm the results of Orosz et al. (2012). Modulations in the radial velocity of KIC 4862625A are measured both spectroscopically and photometrically, i.e. via Doppler boosting, and produce similar results.

Those that do not sow care little about such mundane things as equinoxes or planting seasons, or even crop rotation for that matter. Wherever and whenever the reavers reave, the mood is always foul and the nights are never warm or pleasant. For the rest of the good folks of Westeros, however, a decent grasp of the long-term weather forecast is a necessity. Many a maester have tried to play the Game of Weather Patterns and foretell when to plant those last turnip seeds, hoping for a few more years of balmy respite. Tried and failed. For other than the somewhat vague (if not outright meaningless) omens of "Winter is Coming", their meteorological efforts have been worse than useless. To right that appalling wrong, here we attempt to explain the apparently erratic seasonal changes in the world of G.R.R.M. A natural explanation for such phenomena is the unique behavior of a circumbinary planet. Thus, by speculating that the planet under scrutiny is orbiting a pair of stars, we utilize the power of numerical three-body dynamics to predict that, unfortunately, it is not possible to predict either the length, or the severity of any coming winter. We conclude that, alas, the Maesters were right -- one can only throw their hands in the air in frustration and, defeated by non-analytic solutions, mumble "Coming winter? May be long and nasty (~850 days, T<268K) or may be short and sweet (~600 days, T~273K). Who knows..."

With the increasing number of directly imaged giant exoplanets the current atmosphere models are often not capable of fully explaining the spectra and luminosity of the sources. A particularly challenging component of the atmosphere models is the formation and properties of condensate cloud layers, which fundamentally impact the energetics, opacity, and evolution of the planets. Here we present a suite of techniques that can be used to estimate the level of rotational modulations these planets may show. We propose that the time--resolved observations of such periodic photometric and spectroscopic variations of extrasolar planets due to their rotation can be used as a powerful tool to probe the heterogeneity of their optical surfaces. We address and discuss the following questions: a) what planet properties can be deduced from the light curve and/or spectra, and in particular can we determine rotation periods, spot--coverage, spot colors, spot spectra; b) what is the optimal configuration of instrument/wavelength/temporal sampling required for these measurements; and, c) can principal component analysis be used to invert the light curve and deduce the surface map of the planet. Our simulations describe the expected spectral differences between homogeneous (clear or cloudy) and patchy atmospheres, outline the significance of the dominant absorption features of water, methane, and CO and provide a method to distinguish these two types of atmospheres. Simulated photometry from current and future instruments is used to estimate the level of detectable photometric variations. We conclude that future instruments will be able to recover not only the rotation periods, cloud cover, cloud colors and spectra but even cloud evolution. We also show that a longitudinal map of the planet's atmosphere can be deduced from its disk--integrated light curves.

The paper studies the correction to the distance modulus induced by inhomogeneities and averaged over all directions from a given observer. The inhomogeneities are modeled as mass-compensated voids in random or regular lattices within Swiss-cheese universes. Void radii below 300 Mpc are considered, which are supported by current redshift surveys and limited by the recently observed imprint such voids leave on CMB. The averaging over all directions, performed by numerical ray tracing, is non-perturbative and includes the supernovas inside the voids. Voids aligning along a certain direction produce a cumulative gravitational lensing correction that increases with their number. Such corrections are destroyed by the averaging over all directions, even in non-randomized simple cubic void lattices. At low redshifts, the average correction is not zero but decays with the peculiar velocities and redshift. Its upper bound is provided by the maximal average correction which assumes no random cancelations between different voids. It is described well by a linear perturbation formula and, for the voids considered, is 20% of the correction corresponding to the maximal peculiar velocity. The average correction calculated in random and simple cubic void lattices is severely damped below the predicted maximal one after a single void diameter. That is traced to cancellations between the corrections from the fronts and backs of different voids. All that implies that voids cannot imitate the effect of dark energy unless they have radii and peculiar velocities much larger than the currently observed. The results obtained allow one to readily predict the redshift above which the direction-averaged fluctuation in the Hubble diagram falls below a required precision and suggest a method to extract the background Hubble constant from low redshift data without the need to correct for peculiar velocities.

Using numerical ray tracing, the paper studies how the average distance modulus in an inhomogeneous universe differs from its homogeneous counterpart. The averaging is over all directions from a fixed observer not over all possible observers (cosmic), thus it is more directly applicable to our observations. Unlike previous studies, the averaging is exact, non-perturbative, and includes all possible non-linear effects. The inhomogeneous universes are represented by Sweese-cheese models containing random and simple cubic lattices of mass-compensated voids. The Earth observer is in the homogeneous cheese which has an Einstein - de Sitter metric. For the first time, the averaging is widened to include the supernovas inside the voids by assuming the probability for supernova emission from any comoving volume is proportional to the rest mass in it. Despite the well known argument for photon flux conservation, the average distance modulus correction at low redshifts is not zero due to the peculiar velocities. A formula for the maximum possible average correction as a function of redshift is derived and shown to be in excellent agreement with the numerical results. The actual average correction calculated in random and simple cubic void lattices is severely damped below the predicted maximal average. That is traced to cancelations between the corrections coming from the fronts and backs of different voids at the same redshift from the observer. The calculated correction at low redshifts allows one to readily predict the redshift at which the averaged fluctuation in the Hubble diagram is below a required precision and suggests a method to extract the background Hubble constant from low redshift data without the need to correct for peculiar velocities.

The ever-expanding depth and quality of photometric and spectroscopic observations of stellar populations increase the need for theoretical models in regions of age-composition parameter space that are largely unexplored at present. Stellar evolution models that employ the most advanced physics and cover a wide range of compositions are needed to extract the most information from current observations of both resolved and unresolved stellar populations. The Dartmouth Stellar Evolution Database is a collection of stellar evolution tracks and isochrones that spans a range of [Fe/H] from -2.5 to +0.5, [alpha/Fe] from -0.2 to +0.8 (for [Fe/H] <=0) or +0.2 (for [Fe/H] >0), and initial He mass fractions from Y=0.245 to 0.40. Stellar evolution tracks were computed for masses between 0.1 and 4 Msun, allowing isochrones to be generated for ages as young as 250 Myr. For the range in masses where the core He flash occurs, separate He-burning tracks were computed starting from the zero age horizontal branch. The tracks and isochrones have been transformed to the observational plane in a variety of photometric systems including standard UBV(RI)c, Stromgren uvby, SDSS ugriz, 2MASS JHKs, and HST ACS-WFC and WFPC2. The Dartmouth Stellar Evolution Database is accessible through a website at http://stellar.dartmouth.edu/~models/ where all tracks, isochrones, and additional files can be downloaded.

We model a spiral galaxy by a thin axially symmetric disk that includes both visible and dark matter. The surface mass density of the disk is calculated directly from the rotational velocity curve without extra assumptions. We simplify the standard application of the model. Since most velocity curves are known out to some radius, r_{max}, we extrapolate them by attaching a Keplerian tail. The numerical procedure and the extrapolation are tested with a known toy mass density and shown to reconstruct it with a good precision if r_{max} includes a sufficient part of the velocity curve. Mass density curves are calculated for Milky Way and NGC 3198. We vary the extent of the flat part of the velocity curves from 30 kpc to 200 kpc and show that does not affect appreciably the calculated mass density inside r_{max}=30 kpc. The reconstructed masses for Milky Way are 15 x 10^10 solar masses inside the visible disk and 23 x 10^10 solar masses inside 30 kpc. For NGC 3198, the reconstructed mass inside the visible disk is 6.5 x 10^10 solar masses and 11 x 10^10 solar masses inside 30 kpc. The total galactic masses are roughly proportional to the extent of the flat part of the velocity curves which is currently unknown. The high light-to-mass ratios obtained for the visible disks of the galaxies - 11 solar units for Milky way and 9.3 for NGC 3198 - suggest presence of dark matter. The method is also applied to NGC 3031 - a spiral galaxy with a declining velocity curve in which case it is able to reconstruct both the mass density curve and the total mass (14 x 10^10 solar masses).

The time dependence of the dark energy density can be an important clue to the nature of dark energy in the universe. We show that future supernova data from dedicated telescopes (such as SNAP), when combined with data of nearby supernovae, can be used to determine how the dark energy density $\rho_X(z)$ depends on redshift, if $\rho_X(z)$ is not too close to a constant. For quantitative comparison, we have done an extensive study of a number of dark energy models. Based on these models we have simulated data sets in order to show that we can indeed reconstruct the correct sign of the time dependence of the dark energy density, outside of a degeneracy region centered on $1+w_0 = -w_1 z_{max}/3$ (where $z_{max}$ is the maximum redshift of the survey, e.g., $z_{max}=1.7$ for SNAP). We emphasize that, given the same data, one can obtain much more information about the dark energy density directly (and its time dependence) than about its equation of state.