The Cold Classical Kuiper Belt, a class of small bodies in undisturbed orbits beyond Neptune, are primitive objects preserving information about Solar System formation. The New Horizons spacecraft flew past one of these objects, the 36 km long contact binary (486958) Arrokoth (2014 MU69), in January 2019. Images from the flyby show that Arrokoth has no detectable rings, and no satellites (larger than 180 meters diameter) within a radius of 8000 km, and has a lightly-cratered smooth surface with complex geological features, unlike those on previously visited Solar System bodies. The density of impact craters indicates the surface dates from the formation of the Solar System. The two lobes of the contact binary have closely aligned poles and equators, constraining their accretion mechanism.

Many of the known hot Jupiters are formally unstable to tidal orbital decay. The only hot Jupiter for which orbital decay has been directly detected is WASP-12, for which transit timing measurements spanning more than a decade have revealed that the orbital period is decreasing at a rate of $dP/dt\approx 10^{-9}$, corresponding to a reduced tidal quality factor of about $2\times 10^5$. Here, we present a compilation of transit-timing data for WASP-12 and eleven other systems which are especially favorable for detecting orbital decay: KELT-16; WASP-18, 19, 43, 72, 103, 114, and 122; HAT-P-23; HATS-18; and OGLE-TR-56. For most of these systems we present new data that extend the time baseline over which observations have been performed. None of the systems besides WASP-12 displays convincing evidence for period changes, with typical upper limits on $dP/dt$ on the order of $10^{-9}$ or $10^{-10}$, and lower limits on the reduced tidal quality factor on the order of $10^5$. One possible exception is WASP-19, which shows a statistically significant trend, although it may be a spurious effect of starspot activity. Further observations are encouraged.

WASP-12b is a transiting hot Jupiter on a 1.09-day orbit around a late-F star. Since the planet's discovery in 2008, the time interval between transits has been decreasing by $29\pm 2$ msec year$^{-1}$. This is a possible sign of orbital decay, although the previously available data left open the possibility that the planet's orbit is slightly eccentric and is undergoing apsidal precession. Here, we present new transit and occultation observations that provide more decisive evidence for orbital decay, which is favored over apsidal precession by a $\Delta\mathrm{BIC}$ of 22.3 or Bayes factor of 70,000. We also present new radial-velocity data that rule out the R{\o}mer effect as the cause of the period change. This makes WASP-12 the first planetary system for which we can be confident that the orbit is decaying. The decay timescale for the orbit is $P/\dot{P} = 3.25\pm 0.23$ Myr. Interpreting the decay as the result of tidal dissipation, the modified stellar tidal quality factor is $Q'_\star = 1.8 \times10^{5}$.

As the NASA Transiting Exoplanet Survey Satellite (TESS) fulfills its primary mission, it is executing an unprecedented survey of almost the entire sky: TESS's approved extended mission will likely extend sky coverage to ~94%, including ~60% of the ecliptic. In an accompanying note we demonstrated that `digital tracking' techniques can be used to efficiently `shift-and-stack' TESS full frame images (FFIs) and showed that combining ~1,300 exposures from a TESS sector gives a 50% detection threshold of $I_C\sim 22.0\pm0.5$, raising the possibility that TESS could discover the hypothesized Planet Nine. In this note, we estimate the yield for such a survey and demonstrate that this technique has the potential to discover hundreds of Kuiper Belt Objects, Scattered Disk Objects and Centaurs in TESS FFI data.

As the NASA Transiting Exoplanet Survey Satellite (TESS) fulfills its primary mission it is executing an unprecedented all-sky survey with the potential to discover distant planets in our own solar system, as well as hundreds of Transneptunian Objects (TNOs) and Centaurs. We demonstrate that shift-and-stack techniques can be used to efficiently search the Full-Frame Image (FFI) data from the TESS mission and survey the entire sky for outer Solar System objects down to $\sim22^{nd}$ magnitude.

We consider the potential for the Transiting Exoplanet Survey Satellite (TESS) to detect transit timing variations (TTVs) during both its nominal and extended mission phases. Building on previous estimates of the overall yield of planetary systems from the TESS mission, we predict that during its nominal two-year mission, TESS will observe measurable TTVs in $\sim30$ systems, from which $\mathcal{O}(10)$ planet will get precise mass measurements from TTVs alone, $\sim 5$ planets will have significant constraints placed on their masses from TTVs, and over a dozen systems will be singly-transiting TTV systems. We consider a number of different extended mission scenarios, and predict that in a typical scenario, an extended mission will allow TESS to increase the number of systems with measurable TTVs to a total of $\sim90$, from which $\sim 15$ planets will have precise mass measurements, another $\sim 15$ will have significant constraints placed on their masses, and $\sim60$ will be singly-transiting TTV systems. We also describe how follow-up transit observations of multi-planet systems discovered by the TESS mission can be optimally planned to maximize TTV mass and eccentricity constraints.

Kepler-9, discovered by Holman et al. 2010, was the first system with multiple confirmed transiting planets and the first system to clearly show long-anticipated transit timing variations (TTVs). We review the historical circumstances behind the discovery and characterization of these planets and the publication of Holman et al. 2010. It was the first major novel exoplanet discovery of the Kepler Space Telescope mission. The Kepler pipeline identified two Saturn-radius candidates (called Kepler Objects of Interest or KOIs): KOI-377.01 with a 19-day period and KOI-377.02 with a 39-day period. Even with only 9 transits for KOI-377.01 and 6 of KOI-377.02, the transit times were completely inconsistent with a linear ephemeris and showed strongly anti-correlated variations in transit times. Holman et al. 2010 were able to readily show that these objects were planetary mass, confirming them as bona fide planets Kepler-9b and Kepler-9c. As a multi-transiting system exhibiting strong TTVs, the relative planetary properties (e.g., mass ratio, radius ratio) were strongly constrained, opening a new chapter in comparative planetology. KOI-377.03, a small planet with a 1.5-day period, was not initially discovered by the Kepler pipeline, but was identified during the analysis of the other planets and was later confirmed as Kepler-9d through the BLENDER technique by Torres et al. 2011. Holman et al. 2010 included significant dynamical analysis to characterize Kepler-9's particular TTVs: planets near resonance show large amplitude anti-correlated TTVs with a period corresponding to the rotation of the line of conjunctions and an additional "chopping" signal due to the changing positions of the planets. We also review the updated properties of this system and propose ideas for future investigations.

We describe the digest2 software package, a fast, short-arc orbit classifier for small Solar System bodies. The digest2 algorithm has been serving the community for more than 13 years. The code provides a score, D2, which represents a pseudo-probability that a tracklet belongs to a given Solar System orbit type. digest2 is primarily used as a classifier for Near-Earth Object (NEO) candidates, to identify those to be prioritized for follow-up observation. We describe the historical development of digest2 and demonstrate its use on real and synthetic data. We find that digest2 can accurately and precisely distinguish NEOs from non- NEOs. At the time of detection, 14% of NEO tracklets and 98.5% of non-NEOs tracklets have D2 below the critical value of D2 = 65. 94% of our simulated NEOs achieved the maximum D2 = 100 and 99.6% of NEOs achieved $D2 \ge 65$ at least once during the simulated 10-year timeframe. We demonstrate that D2 varies as a function of time, rate of motion, magnitude and sky-plane location, and show that NEOs tend to have lower D2 at low Solar elongations close to the ecliptic. We use our findings to recommend future development directions for the digest2 code.

Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3\pi$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$\sigma$ point source limiting sensitivities in the stacked 3$\pi$ Steradian Survey in $grizy_{P1}$ are (23.3, 23.2, 23.1, 22.3, 21.4) respectively. The upper bound on the systematic uncertainty in the photometric calibration across the sky is 7-12 millimag depending on the bandpass. The systematic uncertainty of the astrometric calibration using the Gaia frame comes from a comparison of the results with Gaia: the standard deviation of the mean and median residuals ($ \Delta ra, \Delta dec $) are (2.3, 1.7) milliarcsec, and (3.1, 4.8) milliarcsec respectively. The Pan-STARRS system and the design of the PS1 surveys are described and an overview of the resulting image and catalog data products and their basic characteristics are described together with a summary of important results. The images, reduced data products, and derived data products from the Pan-STARRS1 surveys are available to the community from the Mikulski Archive for Space Telescopes (MAST) at STScI.

The existence of Planet Nine has been suggested to explain the pericenter clustering of extreme trans-Neptunian objects (TNOs). However, the underlying dynamics involving Planet Nine, test particles and Neptune is rich, and it remains unclear which dynamical processes lead to the alignment and how they depends on the properties of Planet Nine. Here, we investigate the secular interactions between an eccentric outer perturber and TNOs starting in a near-coplanar configuration. We find that a large number of TNOs could survive outside of mean motion resonances at 4Gyr, which differs from previous results obtained in the exact coplanar case with Neptune being treated as a quadrupole potential. In addition, secular dynamics leads to the orbital clustering seen in N-body simulations. We find that a near coplanar Planet Nine can flip TNO orbital planes, and when this happens, the geometrical longitudes of pericenter of the TNOs librate around $180^\circ$ during the flip. Orbital precession caused by the inner giant planets can suppress the flips while keeping the longitude of pericenter librating when $30\lesssim r_p \lesssim 80$ AU $\&$ $a\gtrsim 250$ AU. This results in the alignment of the pericenter of the low inclination TNOs ($i \lesssim 40^\circ$). We find the anti-aligned population and the flipped orbits could be produced by an eccentric ($e_9 \gtrsim 0.4$) outer planet of $\sim 10M_{\oplus}$ in a wide $a_9 \gtrsim 400 \sim 800$ AU orbit. Future surveys on the high inclination TNOs will help further constrain the properties of possible outer planets.

In this paper we investigate systems previously identified to exhibit transit timing variations (TTVs) in Kepler data, with the goal of predicting the expected improvements to the mass and eccentricity constraints that will arise from combining Kepler data with future data from the TESS mission. We advocate for the use of the Kullback-Leibler (KL) divergence as a means to quantify improvements in the measured constraints. Compared to the original Kepler data, the TESS data will have a lower signal-to-noise ratio, rendering some of the planetary transits undetectable, and lowering the accuracy with which the transit mid-time can be estimated. Despite these difficulties, out of the 55 systems (containing 143 planets) investigated, we predict that the collection of short-cadence data by TESS will be of significant value (i.e. it will improve the mass uncertainty such that the KL divergence is > 0.1) for approximately 6 - 14 planets during the nominal mission, with the range primarily driven by the uncertain precision with which transit mid-times will be recovered from TESS data. In an extended mission this would increase to a total of approximately 12 - 25 planets.

We studied the Near-Earth Asteroid (NEA) candidates posted on the Minor Planet Center's Near-Earth Object Confirmation Page (NEOCP) between years 2013 and 2016. Out of more than 17,000 NEA candidates, while the majority became either new discoveries or were associated with previously known objects, about 11% were unable to be followed-up or confirmed. We further demonstrate that of the unconfirmed candidates, 926+/-50 are likely to be NEAs, representing 18% of discovered NEAs in that period. Only 11% (~93) of the unconfirmed NEA candidates were large (having absolute magnitude H<22). To identify the reasons why these NEAs were not recovered, we analyzed those from the most prolific asteroid surveys: Pan-STARRS, the Catalina Sky Survey, the Dark Energy Survey, and the Space Surveillance Telescope. We examined the influence of plane-of-sky positions and rates of motion, brightnesses, submission delays, and computed absolute magnitudes, as well as correlations with the phase of the moon and seasonal effects. We find that delayed submission of newly discovered NEA candidate to the NEOCP drove a large fraction of the unconfirmed NEA candidates. A high rate of motion was another significant contributing factor. We suggest that prompt submission of suspected NEA discoveries and rapid response to fast moving targets and targets with fast growing ephemeris uncertainty would allow better coordination among dedicated follow-up observers, decrease the number of unconfirmed NEA candidates, and increase the discovery rate of NEAs.

We present a novel approach to the minor planet linking problem. Our heliocentric transformation-and-propagation algorithm clusters tracklets at common epochs, allowing for the efficient identification of tracklets that represent the same minor planet. This algorithm scales as $O(N log N)$, with the number of tracklets $N$, a significant advance over standard methods, which scale as $O(N^3)$. This overcomes one of the primary computational bottlenecks faced by current and future asteroid surveys. We apply our algorithm to the Minor Planet Center's Isolated Tracklet File, establishing orbits for $\sim41,000$ new minor planets.

Observations of clustering among the orbits of the most distant trans-Neptunian objects (TNOs) has inspired interest in the possibility of an undiscovered ninth planet lurking in the outskirts of the solar system. Numerical simulations by a number of authors have demonstrated that, with appropriate choices of planet mass and orbit, such a planet can maintain clustering in the orbital elements of the population of distant TNOs, similar to the observed sample. However, many aspects of the rich underlying dynamical processes induced by such a distant eccentric perturber have not been fully explored. We report the results of our investigation of the dynamics of coplanar test-particles that interact with a massive body on an circular orbit (Neptune) and a massive body on a more distant, highly eccentric orbit (the putative Planet Nine). We find that a detailed examination of our idealized simulations affords tremendous insight into the rich test-particle dynamics that are possible. In particular, we find that chaos and resonance overlap plays an important role in particles' dynamical evolution. We develop a simple mapping model that allows us to understand, in detail, the web of overlapped mean-motion resonances explored by chaotically evolving particles. We also demonstrate that gravitational interactions with Neptune can have profound effects on the orbital evolution of particles. Our results serve as a starting point for a better understanding of the dynamical behavior observed in more complicated simulations that can be used to constrain the mass and orbit of Planet 9.

The Large Synoptic Survey Telescope (LSST) is uniquely equipped to search for Solar System bodies due to its unprecedented combination of depth and wide field coverage. Over a ten-year period starting in 2022, LSST will generate the largest catalog of Solar System objects to date. The main goal of the LSST Solar System Science Collaboration (SSSC) is to facilitate the efforts of the planetary community to study the planets and small body populations residing within our Solar System using LSST data. To prepare for future survey cadence decisions and ensure that interesting and novel Solar System science is achievable with LSST, the SSSC has identified and prioritized key Solar System research areas for investigation with LSST in this roadmap. The ranked science priorities highlighted in this living document will inform LSST survey cadence decisions and aid in identifying software tools and pipelines needed to be developed by the planetary community as added value products and resources before the planned start of LSST science operations.

We report the discovery of a $H_r = 3.4\pm0.1$ dwarf planet candidate by the Pan-STARRS Outer Solar System Survey. 2010 JO$_{179}$ is red with $(g-r)=0.88 \pm 0.21$, roughly round, and slowly rotating, with a period of $30.6$ hr. Estimates of its albedo imply a diameter of 600--900~km. Observations sampling the span between 2005--2016 provide an exceptionally well-determined orbit for 2010 JO$_{179}$, with a semi-major axis of $78.307\pm0.009$ au, distant orbits known to this precision are rare. We find that 2010 JO$_{179}$ librates securely within the 21:5 mean-motion resonance with Neptune on hundred-megayear time scales, joining the small but growing set of known distant dwarf planets on metastable resonant orbits. These imply a substantial trans-Neptunian population that shifts between stability in high-order resonances, the detached population, and the eroding population of the scattering disk.

We present new transit and occultation times for the hot Jupiter WASP-12b. The data are compatible with a constant period derivative: $\dot{P}=-29 \pm 3$ ms yr$^{-1}$ and $P/\dot{P}= 3.2$ Myr. However, it is difficult to tell whether we have observed orbital decay, or a portion of a 14-year apsidal precession cycle. If interpreted as decay, the star's tidal quality parameter $Q_\star$ is about $2\times 10^5$. If interpreted as precession, the planet's Love number is $0.44\pm 0.10$. Orbital decay appears to be the more parsimonious model: it is favored by $\Delta\chi^2=5.5$ despite having two fewer free parameters than the precession model. The decay model implies that WASP-12 was discovered within the final $\sim$0.2% of its existence, which is an unlikely coincidence but harmonizes with independent evidence that the planet is nearing disruption. Precession does not invoke any temporal coincidence, but does require some mechanism to maintain an eccentricity of $\approx$0.002 in the face of rapid tidal circularization. To distinguish unequivocally between decay and precession will probably require a few more years of monitoring. Particularly helpful will be occultation timing in 2019 and thereafter.

Hi'iaka is the larger outer satellite of the dwarf planet Haumea. Using relative photometry from the Hubble Space Telescope and Magellan and a phase dispersion minimization analysis, we have identified the rotation period of Hi'iaka to be ~9.8 hrs (double-peaked). This is ~120 times faster than its orbital period, creating new questions about the formation of this system and possible tidal evolution. The rapid rotation suggests that Hi'iaka could have a significant obliquity and spin precession that could be visible in light curves within a few years. We then turn to an investigation of what we learn about the (presently unclear) formation of the Haumea system and family based on this unexpectedly rapid rotation rate. We explore the importance of the initial semi-major axis and rotation period in tidal evolution theory and find they strongly influence the time required to despin to synchronous rotation, relevant to understanding a wide variety of satellite and binary systems. We find that despinning tides do not necessarily lead to synchronous spin periods for Hi'iaka, even if it formed near the Roche limit. Therefore the short rotation period of Hi'iaka does not rule out significant tidal evolution. Hi'iaka's spin period is also consistent with formation near its current location and spin up due to Haumea-centric impactors.

Roughly 1000 white dwarfs are known to be polluted with planetary material, and the progenitors of this material are typically assumed to be asteroids. The dynamical architectures which perturb asteroids into white dwarfs are still unknown, but may be crucially dependent on moons liberated from parent planets during post-main-sequence gravitational scattering. Here, we trace the fate of these exomoons, and show that they more easily achieve deep radial incursions towards the white dwarf than do scattered planets. Consequently, moons are likely to play a significant role in white dwarf pollution, and in some cases may be the progenitors of the pollution itself.

In this work we report the detection of seven Neptune Trojans (NTs) in the Pan-STARRS 1 (PS1) survey. Five of these are new discoveries, consisting of four L4 Trojans and one L5 Trojan. Our orbital simulations show that the L5 Trojan stably librates for only several million years. This suggests that the L5 Trojan must be of recent capture origin. On the other hand, all four new L4 Trojans stably occupy the 1:1 resonance with Neptune for more than 1 Gyr. They can, therefore, be of primordial origin. Our survey simulation results show that the inclination width of the Neptune Trojan population should be between $7^{\circ}$ and $27^{\circ}$ at $>$ 95% confidence, and most likely $\sim 11^{\circ}$. In this paper, we describe the PS1 survey, the Outer Solar System pipeline, the confirming observations, and the orbital/physical properties of the new Neptune Trojans.

Although the majority of Centaurs are thought to have originated in the scattered disk, with the high-inclination members coming from the Oort cloud, the origin of the high inclination component of trans-Neptunian objects (TNOs) remains uncertain. We report the discovery of a retrograde TNO, which we nickname "Niku", detected by the Pan-STARRS 1 Outer Solar System Survey. Our numerical integrations show that the orbital dynamics of Niku are very similar to that of 2008 KV$_{42}$ (Drac), with a half-life of $\sim 500$ Myr. Comparing similar high inclination TNOs and Centaurs ($q > 10$ AU, $a < 100$ AU and $i > 60^\circ$), we find that these objects exhibit a surprising clustering of ascending node, and occupy a common orbital plane. This orbital configuration has high statistical significance: 3.8-$\sigma$. An unknown mechanism is required to explain the observed clustering. This discovery may provide a pathway to investigate a possible reservoir of high-inclination objects.

The new discoveries of circumbinary planetary systems shed light on the understanding of planetary system formation. Learning the architectural properties of these systems is essential for constraining the different formation mechanisms. We first revisit the stability limit of circumbinary planets. Next, we focus on eclipsing stellar binaries and obtain an analytical expression for the transit probability in a realistic setting, where finite observation period and planetary orbital precession are included. Then, understanding of the architectural properties of the currently observed transiting systems is refined, based on Bayesian analysis and a series of hypothesis tests. We find 1) it is not a selection bias that the innermost planets reside near the stability limit for eight of the nine observed systems, and this is consistent with a log uniform distribution of the planetary semi-major axis; 2) it is not a selection bias that the planetary and stellar orbits are nearly coplanar ($\lesssim 3^\circ$), and this together with previous studies may imply an occurrence rate of circumbinary planets similar to that of single star systems; 3) the dominance of observed circumbinary systems with only one transiting planet may be caused by selection effects; 4) formation mechanisms involving Lidov-Kozai oscillations, which may produce misalignment and large separation between planet and stellar binaries, are consistent with the lack of transiting circumbinary planets around short-period stellar binaries, in agreement with previous studies. As a consequence of 4), eclipse timing variations may better suit the detection of planets in such configurations.

Kepler has discovered hundreds of systems with multiple transiting exoplanets which hold tremendous potential both individually and collectively for understanding the formation and evolution of planetary systems. Many of these systems consist of multiple small planets with periods less than ~50 days known as Systems with Tightly-spaced Inner Planets, or STIPs. One especially intriguing STIP, Kepler-80 (KOI-500), contains five transiting planets: f, d, e, b, and c with periods of 1.0, 3.1, 4.6, 7.1, 9.5 days, respectively. We provide measurements of transit times and a transit timing variation (TTV) dynamical analysis. We find that TTVs cannot reliably detect eccentricities for this system, though mass estimates are not affected. Restricting the eccentricity to a reasonable range, we infer masses for the outer four planets (d, e, b, and c) to be $6.75^{+0.69}_{-0.51}$, $4.13^{+0.81}_{-0.95}$, $6.93^{+1.05}_{-0.70}$, and $6.74^{+1.23}_{-0.86}$ Earth masses, respectively. The similar masses but different radii are consistent with terrestrial compositions for d and e and $\sim$2% H/He envelopes for b and c. We confirm that the outer four planets are in a rare dynamical configuration with four interconnected three-body resonances that are librating with few degree amplitudes. We present a formation model that can reproduce the observed configuration by starting with a multi-resonant chain and introducing dissipation. Overall, the information-rich Kepler-80 planets provide an important perspective into exoplanetary systems.

We significantly constrain the sky position, distance, and mass of a possible additional, distant planet in the solar system by examining its influence on the distance between Earth and the Cassini Spacecraft. Our preferred region is approximately centered on (RA, Dec) = ($40\arcdeg$, $-15\arcdeg$), extending approximately 20 degrees in all directions.

Previous studies indicate that more than a quarter of all white dwarf (WD) atmospheres are polluted by remnant planetary material, with some WDs being observed to accrete the mass of Pluto in 10^6 years. The short sinking timescale for the pollutants indicate that the material must be frequently replenished. Moons may contribute decisively to this pollution process if they are liberated from their parent planets during the post-main-sequence evolution of the planetary systems. Here, we demonstrate that gravitational scattering events among planets in WD systems easily triggers moon ejection. Repeated close encounters within tenths of a planetary Hill radii are highly destructive to even the most massive, close-in moons. Consequently, scattering increases both the frequency of perturbing agents in WD systems, as well as the available mass of polluting material in those systems, thereby enhancing opportunities for collision and fragmentation and providing more dynamical pathways for smaller bodies to reach the WD. Moreover, during intense scattering, planets themselves have pericenters with respect to the WD of only a fraction of an AU, causing extreme Hill-sphere contraction, and the liberation of moons into WD-grazing orbits. Many of our results are directly applicable to exomoons orbiting planets around main sequence stars.

We use astrometry of Pluto and other TNOs to constrain the sky location, distance, and mass of the possible additional planet (Planet Nine) hypothesized by Batygin and Brown (2016). We find that over broad regions of the sky, the inclusion of a massive, distant planet degrades the fits to the observations. However, in other regions, the fits are significantly improved by the addition of such a planet. Our best fits suggest a planet that is either more massive or closer than argued for by Batygin and Brown (2016) based on the orbital distribution of distant trans-neptunian objects (or by Fienga et al. (2016) based on range measured to the Cassini spacecraft). The trend to favor larger and closer perturbing planets is driven by the residuals to the astrometry of Pluto, remeasured from photographic plates using modern stellar catalogs (Buie and Folkner 2015), which show a clear trend in declination, over the course of two decades, that drive a preference for large perturbations. Although this trend may be the result of systematic errors of unknown origin in the observations, a possible resolution is that the declination trend may be due to perturbations from a body, additional to Planet Nine , that is closer to Pluto, but less massive than, Planet Nine .

K2-19 (EPIC201505350) is an interesting planetary system in which two transiting planets with radii ~ 7 $R_{Earth}$ (inner planet b) and ~ 4 $R_{Earth}$ (outer planet c) have orbits that are nearly in a 3:2 mean-motion resonance. Here, we present results of ground-based follow-up observations for the K2-19 planetary system. We have performed high-dispersion spectroscopy and high-contrast adaptive-optics imaging of the host star with the HDS and HiCIAO on the Subaru 8.2m telescope. We find that the host star is relatively old (>8 Gyr) late G-type star ($T_{eff}$ ~ 5350 K, $M_s$ ~ 0.9 $M_{Sun}$, and $R_{s}$ ~ 0.9 $R_{Sun}$). We do not find any contaminating faint objects near the host star which could be responsible for (or dilute) the transit signals. We have also conducted transit follow-up photometry for the inner planet with KeplerCam on the FLWO 1.2m telescope, TRAPPISTCAM on the TRAPPIST 0.6m telescope, and MuSCAT on the OAO 1.88m telescope. We confirm the presence of transit-timing variations, as previously reported by Armstrong and coworkers. We model the observed transit-timing variations of the inner planet using the synodic chopping formulae given by Deck & Agol (2015). We find two statistically indistinguishable solutions for which the period ratios ($P_{c}/P_{b}$) are located slightly above and below the exact 3:2 commensurability. Despite the degeneracy, we derive the orbital period of the inner planet $P_b$ ~ 7.921 days and the mass of the outer planet $M_c$ ~ 20 $M_{Earth}$. Additional transit photometry (especially for the outer planet) as well as precise radial-velocity measurements would be helpful to break the degeneracy and to determine the mass of the inner planet.

The T Tauri star PTFO 8-8695 exhibits periodic fading events that have been interpreted as the transits of a giant planet on a precessing orbit. Here we present three tests of the planet hypothesis. First, we sought evidence for the secular changes in light-curve morphology that are predicted to be a consequence of orbital precession. We observed 28 fading events spread over several years, and did not see the expected changes. Instead we found that the fading events are not strictly periodic. Second, we attempted to detect the planet's radiation, based on infrared observations spanning the predicted times of occultations. We ruled out a signal of the expected amplitude. Third, we attempted to detect the Rossiter-McLaughlin effect by performing high-resolution spectroscopy throughout a fading event. No effect was seen at the expected level, ruling out most (but not all) possible orientations for the hypothetical planetary orbit. Our spectroscopy also revealed strong, time-variable, high-velocity H{\alpha} and Ca H & K emission features. All these observations cast doubt on the planetary hypothesis, and suggest instead that the fading events represent starspots, eclipses by circumstellar dust, or occultations of an accretion hotspot.

We present the discovery of a transiting exoplanet candidate in the K2 Field-1 with an orbital period of 9.1457 hr: K2-22b. The highly variable transit depths, ranging from $\sim$0\% to 1.3\%, are suggestive of a planet that is disintegrating via the emission of dusty effluents. We characterize the host star as an M-dwarf with $T_{\rm eff} \simeq 3800$ K. We have obtained ground-based transit measurements with several 1-m class telescopes and with the GTC. These observations (1) improve the transit ephemeris; (2) confirm the variable nature of the transit depths; (3) indicate variations in the transit shapes; and (4) demonstrate clearly that at least on one occasion the transit depths were significantly wavelength dependent. The latter three effects tend to indicate extinction of starlight by dust rather than by any combination of solid bodies. The K2 observations yield a folded light curve with lower time resolution but with substantially better statistical precision compared with the ground-based observations. We detect a significant "bump" just after the transit egress, and a less significant bump just prior to transit ingress. We interpret these bumps in the context of a planet that is not only likely streaming a dust tail behind it, but also has a more prominent leading dust trail that precedes it. This effect is modeled in terms of dust grains that can escape to beyond the planet's Hill sphere and effectively undergo `Roche lobe overflow,' even though the planet's surface is likely underfilling its Roche lobe by a factor of 2.

The GJ876 system was among the earliest multi-planetary detections outside of the Solar System, and has long been known to harbor a resonant pair of giant planets. Subsequent characterization of the system revealed the presence of an additional Neptune mass object on an external orbit, locked in a three body Laplace mean motion resonance with the previously known planets. While this system is currently the only known extrasolar example of a Laplace resonance, it differs from the Galilean satellites in that the orbital motion of the planets is known to be chaotic. In this work, we present a simple perturbative model that illuminates the origins of stochasticity inherent to this system and derive analytic estimates of the Lyapunov time as well as the chaotic diffusion coefficient. We then address the formation of the multi-resonant structure within a protoplanetary disk and show that modest turbulent forcing in addition to dissipative effects is required to reproduce the observed chaotic configuration. Accordingly, this work places important constraints on the typical formation environments of planetary systems and informs the attributes of representative orbital architectures that arise from extended disk-driven evolution.

On timescales that greatly exceed an orbital period, typical planetary orbits evolve in a stochastic yet stable fashion. On even longer timescales, however, planetary orbits can spontaneously transition from bounded to unbound chaotic states. Large-scale instabilities associated with such behavior appear to play a dominant role in shaping the architectures of planetary systems, including our own. Here we show how such transitions are possible, focusing on the specific case of the long-term evolution of Mercury. We develop a simple analytical model for Mercury's dynamics and elucidate the origins of its short term stochastic behavior as well as of its sudden progression to unbounded chaos. Our model allows us to estimate the timescale on which this transition is likely to be triggered, i.e. the dynamical lifetime of the Solar System as we know it. The formulated theory is consistent with the results of numerical simulations and is broadly applicable to extrasolar planetary systems dominated by secular interactions. These results constitute a significant advancement in our understanding of the processes responsible for sculpting of the dynamical structures of generic planetary systems.

We report the discovery of one newly confirmed planet ($P=66.06$ days, $R_{\rm{P}}=2.68\pm0.17R_\oplus$) and mass determinations of two previously validated Kepler planets, Kepler-289 b ($P=34.55$ days, $R_{\rm{P}}=2.15\pm0.10R_\oplus$) and Kepler-289-c ($P=125.85$ days, $R_{\rm{P}}=11.59\pm0.10R_\oplus$), through their transit timing variations (TTVs). We also exclude the possibility that these three planets reside in a $1:2:4$ Laplace resonance. The outer planet has very deep ($\sim1.3%$), high signal-to-noise transits, which puts extremely tight constraints on its host star's stellar properties via Kepler's Third Law. The star PH3 is a young ($\sim1$ Gyr as determined by isochrones and gyrochronology), Sun-like star with $M_*=1.08\pm0.02M_\odot$, $R_*=1.00\pm0.02R_\odot$, and $T_{\rm{eff}}=5990\pm38$ K. The middle planet's large TTV amplitude ($\sim5$ hours) resulted either in non-detections or inaccurate detections in previous searches. A strong chopping signal, a shorter period sinusoid in the TTVs, allows us to break the mass-eccentricity degeneracy and uniquely determine the masses of the inner, middle, and outer planets to be $M=7.3\pm6.8M_\oplus$, $4.0\pm0.9M_\oplus$, and $M=132\pm17M_\oplus$, which we designate PH3 b, c, and d, respectively. Furthermore, the middle planet, PH3 c, has a relatively low density, $\rho=1.2\pm0.3$ g/cm$^3$ for a planet of its mass, requiring a substantial H/He atmosphere of $2.1^{+0.8}_{-0.3}%$ by mass, and joins a growing population of low-mass, low-density planets.

The Transiting Exoplanet Survey Satellite (TESS) will search for planets transiting bright and nearby stars. TESS has been selected by NASA for launch in 2017 as an Astrophysics Explorer mission. The spacecraft will be placed into a highly elliptical 13.7-day orbit around the Earth. During its two-year mission, TESS will employ four wide-field optical CCD cameras to monitor at least 200,000 main-sequence dwarf stars with I = 4-13 for temporary drops in brightness caused by planetary transits. Each star will be observed for an interval ranging from one month to one year, depending mainly on the star's ecliptic latitude. The longest observing intervals will be for stars near the ecliptic poles, which are the optimal locations for follow-up observations with the James Webb Space Telescope. Brightness measurements of preselected target stars will be recorded every 2 min, and full frame images will be recorded every 30 min. TESS stars will be 10-100 times brighter than those surveyed by the pioneering Kepler mission. This will make TESS planets easier to characterize with follow-up observations. TESS is expected to find more than a thousand planets smaller than Neptune, including dozens that are comparable in size to the Earth. Public data releases will occur every four months, inviting immediate community-wide efforts to study the new planets. The TESS legacy will be a catalog of the nearest and brightest stars hosting transiting planets, which will endure as highly favorable targets for detailed investigations.

The Hubble Space Telescope (HST) and the Kepler space mission observed a large number of planetary transits showing anomalies due to starspot eclipses, with more such observations expected in the near future by the K2 mission and the Transiting Exoplanet Survey Satellite (TESS). To facilitate analysis of this phenomenon, we present spotrod, a model for planetary transits of stars with an arbitrary limb darkening law and a number of homogeneous, circular spots on their surface. A free, open source implementation written in C, ready to use in Python, is available for download. We analyze Kepler observations of the planetary host star HAT-P-11, and study the size and contrast of more than two hundred starspots. We find that the flux ratio of spots ranges at least from 0.6 to 0.9, corresponding to an effective temperature approximately 100 to 450 K lower than the stellar surface, although it is possible that some spots are darker than 0.5. The largest detected spots have a radius less than approximately 0.2 stellar radii.

We report on the orbital architectures of Kepler systems having multiple planet candidates identified in the analysis of data from the first six quarters of Kepler data and reported by Batalha et al. (2013). These data show 899 transiting planet candidates in 365 multiple-planet systems and provide a powerful means to study the statistical properties of planetary systems. Using a generic mass-radius relationship, we find that only two pairs of planets in these candidate systems (out of 761 pairs total) appear to be on Hill-unstable orbits, indicating ~96% of the candidate planetary systems are correctly interpreted as true systems. We find that planet pairs show little statistical preference to be near mean-motion resonances. We identify an asymmetry in the distribution of period ratios near first-order resonances (e.g., 2:1, 3:2), with an excess of planet pairs lying wide of resonance and relatively few lying narrow of resonance. Finally, based upon the transit duration ratios of adjacent planets in each system, we find that the interior planet tends to have a smaller transit impact parameter than the exterior planet does. This finding suggests that the mode of the mutual inclinations of planetary orbital planes is in the range 1.0-2.2 degrees, for the packed systems of small planets probed by these observations.

The Kozai-Lidov mechanism can be applied to a vast variety of astrophysical systems involving hierarchical three-body systems. Here, we study the Kozai-Lidov mechanism systematically in the test particle limit at the octupole level of approximation. We investigate the chaotic and quasiperiodic orbital evolution by studying surfaces of section and the Lyapunov exponents. We find that the resonances introduced by the octupole level of approximation cause orbits to flip from prograde to retrograde and back as well as cause significant eccentricity excitation, and the chaotic behaviors occur when the mutual inclination between the inner and the outer binary is high. We characterize the parameter space that allows large amplitude oscillations in eccentricity and inclination.

A number of planet-host stars have been observed to rotate with a period equal to an integer multiple of the orbital period of their close planet. We expand this list by analyzing Kepler data of HAT-P-11 and finding a period ratio of 6:1. In particular, we present evidence for a long-lived spot on the stellar surface that is eclipsed by the planet in the same position four times, every sixth transit. We also identify minima in the out-of-transit lightcurve and confirm that their phase with respect to the stellar rotation is mostly stationary for the 48-month timeframe of the observations, confirming the proposed rotation period. For comparison, we apply our methods to Kepler-17 and confirm the findings of Bonomo & Lanza (2012) that the period ratio is not exactly 8:1 in that system. Finally, we provide a hypothesis on how interactions between a star and its planet could possibly result in an observed commensurability for systems where the stellar differential rotation profile happens to include a period at some latitude which is commensurable to the planetary orbit.

Transit timing variations (TTVs) have proven to be a powerful technique for confirming Kepler planet candidates, for detecting non-transiting planets, and for constraining the masses and orbital elements of multi-planet systems. These TTV applications often require the numerical integration of orbits for computation of transit times (as well as impact parameters and durations); frequently tens of millions to billions of simulations are required when running statistical analyses of the planetary system properties. We have created a fast code for transit timing computation, TTVFast, which uses a symplectic integrator with a Keplerian interpolator for the calculation of transit times (Nesvorny et al. 2013). The speed comes at the expense of accuracy in the calculated times, but the accuracy lost is largely unnecessary, as transit times do not need to be calculated to accuracies significantly smaller than the measurement uncertainties on the times. The time step can be tuned to give sufficient precision for any particular system. We find a speed-up of at least an order of magnitude relative to dynamical integrations with high precision using a Bulirsch-Stoer integrator.

P/2011 S1 (Gibbs) is an outer solar system comet or active Centaur with a similar orbit to that of the famous 29P/Schwassmann-Wachmann 1. P/2011 S1 (Gibbs) has been observed by the Pan-STARRS 1 (PS1) sky survey from 2010 to 2012. The resulting data allow us to perform multi-color studies of the nucleus and coma of the comet. Analysis of PS1 images reveals that P/2011 S1 (Gibbs) has a small nucleus $< 4$ km radius, with colors $g_{P1}-r_{P1} = 0.5 \pm 0.02$, $r_{P1}-i_{P1} = 0.12 \pm 0.02$ and $i_{P1}-z_{P1} = 0.46 \pm 0.03$. The comet remained active from 2010 to 2012, with a model-dependent mass-loss rate of $\sim100$ kg s$^{-1}$. The mass-loss rate per unit surface area of P/2011 S1 (Gibbs) is as high as that of 29P/Schwassmann-Wachmann 1, making it one of the most active Centaurs. The mass-loss rate also varies with time from $\sim 40$ kg s$^{-1}$ to 150 kg s$^{-1}$. Due to its rather circular orbit, we propose that P/2011 S1 (Gibbs) has 29P/Schwassmann-Wachmann 1-like outbursts that control the outgassing rate. The results indicate that it may have a similar surface composition to that of 29P/Schwassmann-Wachmann 1. Our numerical simulations show that the future orbital evolution of P/2011 S1 (Gibbs) is more similar to that of the main population of Centaurs than to that of 29P/Schwassmann-Wachmann 1. The results also demonstrate that P/2011 S1 (Gibbs) is dynamically unstable and can only remain near its current orbit for roughly a thousand years.

We report on the masses, sizes, and orbits of the planets orbiting 22 Kepler stars. There are 49 planet candidates around these stars, including 42 detected through transits and 7 revealed by precise Doppler measurements of the host stars. Based on an analysis of the Kepler brightness measurements, along with high-resolution imaging and spectroscopy, Doppler spectroscopy, and (for 11 stars) asteroseismology, we establish low false-positive probabilities for all of the transiting planets (41 of 42 have a false-positive probability under 1%), and we constrain their sizes and masses. Most of the transiting planets are smaller than 3X the size of Earth. For 16 planets, the Doppler signal was securely detected, providing a direct measurement of the planet's mass. For the other 26 planets we provide either marginal mass measurements or upper limits to their masses and densities; in many cases we can rule out a rocky composition. We identify 6 planets with densities above 5 g/cc, suggesting a mostly rocky interior for them. Indeed, the only planets that are compatible with a purely rocky composition are smaller than ~2 R_earth. Larger planets evidently contain a larger fraction of low-density material (H, He, and H2O).

There have been previous hints that the transiting planet WASP-3 b is accompanied by a second planet in a nearby orbit, based on small deviations from strict periodicity of the observed transits. Here we present 17 precise radial velocity measurements and 32 transit light curves that were acquired between 2009 and 2011. These data were used to refine the parameters of the host star and transiting planet. This has resulted in reduced uncertainties for the radii and masses of the star and planet. The radial-velocity data and the transit times show no evidence for an additional planet in the system. Therefore, we have determined the upper limit on the mass of any hypothetical second planet, as a function of its orbital period.

We perform numerical integrations of four-body (star, planet, planet, satellite) systems to investigate the stability of satellites in planetary Systems with Tightly-packed Inner Planets (STIPs). We find that the majority of closely-spaced stable two-planet systems can stably support satellites across a range of parameter-space which is only slightly decreased compared to that seen for the single-planet case. In particular, circular prograde satellites remain stable out to $\sim 0.4 R_H$ (where $R_H$ is the Hill Radius) as opposed to $\sim 0.5 R_H$ in the single-planet case. A similarly small restriction in the stable parameter-space for retrograde satellites is observed, where planetary close approaches in the range 2.5 to 4.5 mutual Hill radii destabilize most satellites orbits only if $a\sim 0.65 R_H$. In very close planetary pairs (e.g. the 12:11 resonance) the addition of a satellite frequently destabilizes the entire system, causing extreme close-approaches and the loss of satellites over a range of circumplanetary semi-major axes. The majority of systems investigated stably harbored satellites over a wide parameter-space, suggesting that STIPs can generally offer a dynamically stable home for satellites, albeit with a slightly smaller stable parameter-space than the single-planet case. As we demonstrate that multi-planet systems are not a priori poor candidates for hosting satellites, future measurements of satellite occurrence rates in multi-planet systems versus single-planet systems could be used to constrain either satellite formation or past periods of strong dynamical interaction between planets.

Motivated by the population of multi-planet systems with orbital period ratios 1<P2/P1<2, we study the long-term stability of packed two planet systems. The Hamiltonian for two massive planets on nearly circular and nearly coplanar orbits near a first order mean motion resonance can be reduced to a one degree of freedom problem (Sessin & Ferraz Mello (1984), Wisdom (1986), Henrard et al. (1986)). Using this analytically tractable Hamiltonian, we apply the resonance overlap criterion to predict the onset of large scale chaotic motion in close two planet systems. The reduced Hamiltonian has only a weak dependence on the planetary mass ratio, and hence the overlap criterion is independent of the planetary mass ratio at lowest order. Numerical integrations confirm that the planetary mass ratio has little effect on the structure of the chaotic phase space for close orbits in the low eccentricity (e <~0.1) regime. We show numerically that orbits in the chaotic web produced primarily by first order resonance overlap eventually experience large scale erratic variation in semimajor axes and are Lagrange unstable. This is also true of the orbits in this overlap region which are Hill stable. As a result, we can use the first order resonance overlap criterion as an effective stability criterion for pairs of observed planets. We show that for low mass (<~10 M_Earth) planetary systems with initially circular orbits the period ratio at which complete overlap occurs and widespread chaos results lies in a region of parameter space which is Hill stable. Our work indicates that a resonance overlap criterion which would apply for initially eccentric orbits needs to take into account second order resonances. Finally, we address the connection found in previous work between the Hill stability criterion and numerically determined Lagrange instability boundaries in the context of resonance overlap.

We present HST STIS observations of two occultations of the transiting exoplanet HAT-P-1b. By measuring the planet to star flux ratio near opposition, we constrain the geometric albedo of the planet, which is strongly linked to its atmospheric temperature gradient. An advantage of HAT-P-1 as a target is its binary companion ADS 16402 A, which provides an excellent photometric reference, simplifying the usual steps in removing instrumental artifacts from HST time-series photometry. We find that without this reference star, we would need to detrend the lightcurve with the time of the exposures as well as the first three powers of HST orbital phase, and this would introduce a strong bias in the results for the albedo. However, with this reference star, we only need to detrend the data with the time of the exposures to achieve the same per-point scatter, therefore we can avoid most of the bias associated with detrending. Our final result is a 2 sigma upper limit of 0.64 for the geometric albedo of HAT-P-1b between 577 and 947 nm.

We describe the Pan-STARRS Moving Object Processing System (MOPS), a modern software package that produces automatic asteroid discoveries and identifications from catalogs of transient detections from next-generation astronomical survey telescopes. MOPS achieves > 99.5% efficiency in producing orbits from a synthetic but realistic population of asteroids whose measurements were simulated for a Pan-STARRS4-class telescope. Additionally, using a non-physical grid population, we demonstrate that MOPS can detect populations of currently unknown objects such as interstellar asteroids. MOPS has been adapted successfully to the prototype Pan-STARRS1 telescope despite differences in expected false detection rates, fill-factor loss and relatively sparse observing cadence compared to a hypothetical Pan-STARRS4 telescope and survey. MOPS remains >99.5% efficient at detecting objects on a single night but drops to 80% efficiency at producing orbits for objects detected on multiple nights. This loss is primarily due to configurable MOPS processing limits that are not yet tuned for the Pan-STARRS1 mission. The core MOPS software package is the product of more than 15 person-years of software development and incorporates countless additional years of effort in third-party software to perform lower-level functions such as spatial searching or orbit determination. We describe the high-level design of MOPS and essential subcomponents, the suitability of MOPS for other survey programs, and suggest a road map for future MOPS development.

We present the discovery of a long-term stable L5 (trailing) Neptune Trojan in data acquired to search for candidate Trans-Neptunian objects for the New Horizons spacecraft to fly by during an extended post-Pluto mission. This Neptune Trojan, 2011 HM102, has the highest inclination (29.4 degrees) of any known member of this population. It is intrinsically brighter than any single L5 Jupiter Trojan at H~8.18. We have determined its gri colors (a first for any L5 Neptune Trojan), which we find to be similar to the moderately red colors of the L4 Neptune Trojans, suggesting similar surface properties for members of both Trojan clouds. We also present colors derived from archival data for two L4 Neptune Trojans (2006 RJ103 and 2007 VL305), better refining the overall color distribution of the population. In this document we describe the discovery circumstances, our physical characterization of 2011 HM102, and this object's implications for the Neptune Trojan population overall. Finally, we discuss the prospects for detecting 2011 HM102 from the New Horizons spacecraft during their close approach in mid- to late-2013.

We present 10 new transit light curves of the transiting hot Jupiter HAT-P-13b, obtained during two observational seasons by three different telescopes. When combined with 12 previously published light curves, we have a sample consisting of 22 transit light curves, spanning 1,041 days across four observational seasons. We use this sample to examine the recently observed large-amplitude transit timing variations (P\'al et al. 2011), and give refined system parameters. We find that the transit times are consistent with a linear ephemeris, with the exception of a single transit time, from UT 2009 Nov 5, for which the measured mid transit time significantly deviates from our linear ephemeris. The nature of this deviation is not clear, and the rest of the data does not show any significant transit timing variation.

Asteroid detections in astronomical images may appear as trails due to a combination of their apparent rate of motion and exposure duration. Nearby asteroids in particular typically have high apparent rates of motion and acceleration. Their recovery, especially on their discovery apparition, depends upon obtaining good astrometry from the trailed detections. We present an analytic function describing a trailed detection under the assumption of a Gaussian point spread function (PSF) and constant rate of motion. We have fit the function to both synthetic and real trailed asteroid detections from the Pan-STARRS1 survey telescope to obtain accurate astrometry and photometry. For short trails our trailing function yields the same astrometric and photometry accuracy as a functionally simpler 2-d Gaussian but the latter underestimates the length of the trail - a parameter that can be important for measuring the object's rate of motion and assessing its cometary activity. For trails longer than about 10 pixels (> 3xPSF) our trail fitting provides 3-times better astrometric accuracy and up to 2 magnitudes improvement in the photometry. The trail fitting algorithm can be implemented at the source detection level for all detections to provide trail length and position angle that can be used to reduce the false tracklet rate.

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.

Unlike Trojans, horseshoe coorbitals are not generally considered to be long-term stable (Dermott and Murray, 1981; Murray and Dermott, 1999). As the lifetime of Earth's and Venus's horseshoe coorbitals is expected to be about a Gyr, we investigated the possible contribution of late-escaping inner planet coorbitals to the lunar Late Heavy Bombardment. Contrary to analytical estimates, we do not find many horseshoe objects escaping after first 100 Myr. In order to understand this behaviour, we ran a second set of simulations featuring idealized planets on circular orbits with a range of masses. We find that horseshoe coorbitals are generally long lived (and potentially stable) for systems with primary-to-secondary mass ratios larger than about 1200. This is consistent with results of Laughlin and Chambers (2002) for equal-mass pairs or coorbital planets and the instability of Jupiter's horseshoe companions (Stacey and Connors, 2008). Horseshoe orbits at smaller mass ratios are unstable because they must approach within 5 Hill radii of the secondary. In contrast, tadpole orbits are more robust and can remain stable even when approaching within 4 Hill radii of the secondary.

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.

We report homogeneous spectroscopic determinations of the effective temperature, metallicity, and projected rotational velocity for the host stars of 56 transiting planets. Our analysis is based primarily on the Stellar Parameter Classification (SPC) technique. We investigate systematic errors by examining subsets of the data with two other methods that have often been used in previous studies (SME and MOOG). The SPC and SME results, both based on comparisons between synthetic spectra and actual spectra, show strong correlations between temperature, [Fe/H], and log g when solving for all three quantities simultaneously. In contrast the MOOG results, based on a more traditional curve-of-growth approach, show no such correlations. To combat the correlations and improve the accuracy of the temperatures and metallicities, we repeat the SPC analysis with a constraint on log g based on the mean stellar density that can be derived from the analysis of the transit light curves. Previous studies that have not taken advantage of this constraint have been subject to systematic errors in the stellar masses and radii of up to 20% and 10%, respectively, which can be larger than other observational uncertainties, and which also cause systematic errors in the planetary mass and radius.

We analyze the deviations of transit times from a linear ephemeris for the Kepler Objects of Interest (KOI) through Quarter six (Q6) of science data. We conduct two statistical tests for all KOIs and a related statistical test for all pairs of KOIs in multi-transiting systems. These tests identify several systems which show potentially interesting transit timing variations (TTVs). Strong TTV systems have been valuable for the confirmation of planets and their mass measurements. Many of the systems identified in this study should prove fruitful for detailed TTV studies.

The Sun's equator and the planets' orbital planes are nearly aligned, which is presumably a consequence of their formation from a single spinning gaseous disk. For exoplanetary systems this well-aligned configuration is not guaranteed: dynamical interactions may tilt planetary orbits, or stars may be misaligned with the protoplanetary disk through chaotic accretion, magnetic interactions or torques from neighbouring stars. Indeed, isolated 'hot Jupiters' are often misaligned and even orbiting retrograde. Here we report an analysis of transits of planets over starspots on the Sun-like star Kepler-30, and show that the orbits of its three planets are aligned with the stellar equator. Furthermore, the orbits are aligned with one another to within a few degrees. This configuration is similar to that of our Solar System, and contrasts with the isolated hot Jupiters. The orderly alignment seen in the Kepler-30 system suggests that high obliquities are confined to systems that experienced disruptive dynamical interactions. Should this be corroborated by observations of other coplanar multi-planet systems, then star-disk misalignments would be ruled out as the explanation for the high obliquities of hot Jupiters, and dynamical interactions would be implicated as the origin of hot Jupiters.

Transit timing variations provide a powerful tool for confirming and characterizing transiting planets, as well as detecting non-transiting planets. We report the results an updated TTV analysis for 1481 planet candidates (Borucki et al. 2011; Batalha et al. 2012) based on transit times measured during the first sixteen months of Kepler observations. We present 39 strong TTV candidates based on long-term trends (2.8% of suitable data sets). We present another 136 weaker TTV candidates (9.8% of suitable data sets) based on excess scatter of TTV measurements about a linear ephemeris. We anticipate that several of these planet candidates could be confirmed and perhaps characterized with more detailed TTV analyses using publicly available Kepler observations. For many others, Kepler has observed a long-term TTV trend, but an extended Kepler mission will be required to characterize the system via TTVs. We find that the occurrence rate of planet candidates that show TTVs is significantly increased (~68%) for planet candidates transiting stars with multiple transiting planet candidate when compared to planet candidates transiting stars with a single transiting planet candidate.

In the Solar system the planets' compositions vary with orbital distance, with rocky planets in close orbits and lower-density gas giants in wider orbits. The detection of close-in giant planets around other stars was the first clue that this pattern is not universal, and that planets' orbits can change substantially after their formation. Here we report another violation of the orbit-composition pattern: two planets orbiting the same star with orbital distances differing by only 10%, and densities differing by a factor of 8. One planet is likely a rocky `super-Earth', whereas the other is more akin to Neptune. These planets are thirty times more closely spaced--and have a larger density contrast--than any adjacent pair of planets in the Solar system.

We report on the long-term dynamical evolution of the two-planet Kepler-36 system, which we studied through numerical integrations of initial conditions that are consistent with observations of the system. The orbits are chaotic with a Lyapunov time of only ~10 years. The chaos is a consequence of a particular set of orbital resonances, with the inner planet orbiting 34 times for every 29 orbits of the outer planet. The rapidity of the chaos is due to the interaction of the 29:34 resonance with the nearby first order 6:7 resonance, in contrast to the usual case in which secular terms in the Hamiltonian play a dominant role. Only one contiguous region of phase space, accounting for ~4.5% of the sample of initial conditions studied, corresponds to planetary orbits that do not show large scale orbital instabilities on the timescale of our integrations (~200 million years). The long-lived subset of the allowed initial conditions are those that satisfy the Hill stability criterion by the largest margin. Any successful theory for the formation of this system will need to account for why its current state is so close to unstable regions of phase space.

We present multi-band optical photometry of 94 spectroscopically-confirmed Type Ia supernovae (SN Ia) in the redshift range 0.0055 to 0.073, obtained between 2006 and 2011. There are a total of 5522 light curve points. We show that our natural system SN photometry has a precision of roughly 0.03 mag or better in BVr'i', 0.06 mag in u', and 0.07 mag in U for points brighter than 17.5 mag and estimate that it has a systematic uncertainty of 0.014, 0.010, 0.012, 0.014, 0.046, and 0.073 mag in BVr'i'u'U, respectively. Comparisons of our standard system photometry with published SN Ia light curves and comparison stars reveal mean agreement across samples in the range of ~0.00-0.03 mag. We discuss the recent measurements of our telescope-plus-detector throughput by direct monochromatic illumination by Cramer et al (in prep.). This technique measures the whole optical path through the telescope, auxiliary optics, filters, and detector under the same conditions used to make SN measurements. Extremely well-characterized natural-system passbands (both in wavelength and over time) are crucial for the next generation of SN Ia photometry to reach the 0.01 mag accuracy level. The current sample of low-z SN Ia is now sufficiently large to remove most of the statistical sampling error from the dark energy error budget. But pursuing the dark-energy systematic errors by determining highly-accurate detector passbands, combining optical and near-infrared (NIR) photometry and spectra, using the nearby sample to illuminate the population properties of SN Ia, and measuring the local departures from the Hubble flow will benefit from larger, carefully measured nearby samples.

We present the results of a search for planetary companions orbiting near hot Jupiter planet candidates (Jupiter-size candidates with orbital periods near 3 days) identified in the Kepler data through its sixth quarter of science operations. Special emphasis is given to companions between the 2:1 interior and exterior mean-motion resonances. A photometric transit search excludes companions with sizes ranging from roughly 2/3 to 5 times the size of the Earth, depending upon the noise properties of the target star. A search for dynamically induced deviations from a constant period (transit timing variations or TTVs) also shows no significant signals. In contrast, comparison studies of warm Jupiters (with slightly larger orbits) and hot Neptune-size candidates do exhibit signatures of additional companions with these same tests. These differences between hot Jupiters and other planetary systems denote a distinctly different formation or dynamical history.

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.

Eighty planetary systems of two or more planets are known to orbit stars other than the Sun. For most, the data can be sufficiently explained by non-interacting Keplerian orbits, so the dynamical interactions of these systems have not been observed. Here we present 4 sets of lightcurves from the Kepler spacecraft, which each show multiple planets transiting the same star. Departure of the timing of these transits from strict periodicity indicates the planets are perturbing each other: the observed timing variations match the forcing frequency of the other planet. This confirms that these objects are in the same system. Next we limit their masses to the planetary regime by requiring the system remain stable for astronomical timescales. Finally, we report dynamical fits to the transit times, yielding possible values for the planets' masses and eccentricities. As the timespan of timing data increases, dynamical fits may allow detailed constraints on the systems' architectures, even in cases for which high-precision Doppler follow-up is impractical.

New transiting planet candidates are identified in sixteen months (May 2009 - September 2010) of data from the Kepler spacecraft. Nearly five thousand periodic transit-like signals are vetted against astrophysical and instrumental false positives yielding 1,091 viable new planet candidates, bringing the total count up to over 2,300. Improved vetting metrics are employed, contributing to higher catalog reliability. Most notable is the noise-weighted robust averaging of multi-quarter photo-center offsets derived from difference image analysis which identifies likely background eclipsing binaries. Twenty-two months of photometry are used for the purpose of characterizing each of the new candidates. Ephemerides (transit epoch, T_0, and orbital period, P) are tabulated as well as the products of light curve modeling: reduced radius (Rp/R*), reduced semi-major axis (d/R*), and impact parameter (b). The largest fractional increases are seen for the smallest planet candidates (197% for candidates smaller than 2Re compared to 52% for candidates larger than 2Re) and those at longer orbital periods (123% for candidates outside of 50-day orbits versus 85% for candidates inside of 50-day orbits). The gains are larger than expected from increasing the observing window from thirteen months (Quarter 1-- Quarter 5) to sixteen months (Quarter 1 -- Quarter 6). This demonstrates the benefit of continued development of pipeline analysis software. The fraction of all host stars with multiple candidates has grown from 17% to 20%, and the paucity of short-period giant planets in multiple systems is still evident. The progression toward smaller planets at longer orbital periods with each new catalog release suggests that Earth-size planets in the Habitable Zone are forthcoming if, indeed, such planets are abundant.

We present the discovery of the Kepler-20 planetary system, which we initially identified through the detection of five distinct periodic transit signals in the Kepler light curve of the host star 2MASSJ19104752+4220194. We find a stellar effective temperature Teff=5455+-100K, a metallicity of [Fe/H]=0.01+-0.04, and a surface gravity of log(g)=4.4+-0.1. Combined with an estimate of the stellar density from the transit light curves we deduce a stellar mass of Mstar=0.912+-0.034 Msun and a stellar radius of Rstar=0.944^{+0.060}_{-0.095} Rsun. For three of the transit signals, our results strongly disfavor the possibility that these result from astrophysical false positives. We conclude that the planetary scenario is more likely than that of an astrophysical false positive by a factor of 2e5 (Kepler-20b), 1e5 (Kepler-20c), and 1.1e3 (Kepler-20d), sufficient to validate these objects as planetary companions. For Kepler-20c and Kepler-20d, the blend scenario is independently disfavored by the achromaticity of the transit: From Spitzer data gathered at 4.5um, we infer a ratio of the planetary to stellar radii of 0.075+-0.015 (Kepler-20c) and 0.065+-0.011 (Kepler-20d), consistent with each of the depths measured in the Kepler optical bandpass. We determine the orbital periods and physical radii of the three confirmed planets to be 3.70d and 1.91^{+0.12}_{-0.21} Rearth for Kepler-20b, 10.85 d and 3.07^{+0.20}_{-0.31} Rearth for Kepelr-20c, and 77.61 d and 2.75^{+0.17}_{-0.30} Rearth for Kepler-20d. From multi-epoch radial velocities, we determine the masses of Kepler-20b and Kepler-20c to be 8.7\+-2.2 Mearth and 16.1+-3.5 Mearth, respectively, and we place an upper limit on the mass of Kepler-20d of 20.1 Mearth (2 sigma).

We present a statistical analysis that demonstrates that the overwhelming majority of Kepler candidate multiple transiting systems (multis) indeed represent true, physically-associated transiting planets. Binary stars provide the primary source of false positives among Kepler planet candidates, implying that false positives should be nearly randomly-distributed among Kepler targets. In contrast, true transiting planets would appear clustered around a smaller number of Kepler targets if detectable planets tend to come in systems and/or if the orbital planes of planets encircling the same star are correlated. There are more than one hundred times as many Kepler planet candidates in multi-candidate systems as would be predicted from a random distribution of candidates, implying that the vast majority are true planets. Most of these multis are multiple planet systems orbiting the Kepler target star, but there are likely cases where (a) the planetary system orbits a fainter star, and the planets are thus significantly larger than has been estimated, or (b) the planets orbit different stars within a binary/multiple star system. We use the low overall false positive rate among Kepler multis, together with analysis of Kepler spacecraft and ground-based data, to validate the closely-packed Kepler-33 planetary system, which orbits a star that has evolved somewhat off of the main sequence. Kepler-33 hosts five transiting planets with periods ranging from 5.67 to 41 days.

We present a new method for confirming transiting planets based on the combination of transit timingn variations (TTVs) and dynamical stability. Correlated TTVs provide evidence that the pair of bodies are in the same physical system. Orbital stability provides upper limits for the masses of the transiting companions that are in the planetary regime. This paper describes a non-parametric technique for quantifying the statistical significance of TTVs based on the correlation of two TTV data sets. We apply this method to an analysis of the transit timing variations of two stars with multiple transiting planet candidates identified by Kepler. We confirm four transiting planets in two multiple planet systems based on their TTVs and the constraints imposed by dynamical stability. An additional three candidates in these same systems are not confirmed as planets, but are likely to be validated as real planets once further observations and analyses are possible. If all were confirmed, these systems would be near 4:6:9 and 2:4:6:9 period commensurabilities. Our results demonstrate that TTVs provide a powerful tool for confirming transiting planets, including low-mass planets and planets around faint stars for which Doppler follow-up is not practical with existing facilities. Continued Kepler observations will dramatically improve the constraints on the planet masses and orbits and provide sensitivity for detecting additional non-transiting planets. If Kepler observations were extended to eight years, then a similar analysis could likely confirm systems with multiple closely spaced, small transiting planets in or near the habitable zone of solar-type stars.

We present a method to confirm the planetary nature of objects in systems with multiple transiting exoplanet candidates. This method involves a Fourier-Domain analysis of the deviations in the transit times from a constant period that result from dynamical interactions within the system. The combination of observed anti-correlations in the transit times and mass constraints from dynamical stability allow us to claim the discovery of four planetary systems Kepler-25, Kepler-26, Kepler-27, and Kepler-28, containing eight planets and one additional planet candidate.

Since the discovery of the first extrasolar giant planets around Sun-like stars, evolving observational capabilities have brought us closer to the detection of true Earth analogues. The size of an exoplanet can be determined when it periodically passes in front of (transits) its parent star, causing a decrease in starlight proportional to its radius. The smallest exoplanet hitherto discovered has a radius 1.42 times that of the Earth's radius (R Earth), and hence has 2.9 times its volume. Here we report the discovery of two planets, one Earth-sized (1.03R Earth) and the other smaller than the Earth (0.87R Earth), orbiting the star Kepler-20, which is already known to host three other, larger, transiting planets. The gravitational pull of the new planets on the parent star is too small to measure with current instrumentation. We apply a statistical method to show that the likelihood of the planetary interpretation of the transit signals is more than three orders of magnitude larger than that of the alternative hypothesis that the signals result from an eclipsing binary star. Theoretical considerations imply that these planets are rocky, with a composition of iron and silicate. The outer planet could have developed a thick water vapour atmosphere.

A search of the time-series photometry from NASA's Kepler spacecraft reveals a transiting planet candidate orbiting the 11th magnitude G5 dwarf KIC 10593626 with a period of 290 days. The characteristics of the host star are well constrained by high-resolution spectroscopy combined with an asteroseismic analysis of the Kepler photometry, leading to an estimated mass and radius of 0.970 +/- 0.060 MSun and 0.979 +/- 0.020 RSun. The depth of 492 +/- 10ppm for the three observed transits yields a radius of 2.38 +/- 0.13 REarth for the planet. The system passes a battery of tests for false positives, including reconnaissance spectroscopy, high-resolution imaging, and centroid motion. A full BLENDER analysis provides further validation of the planet interpretation by showing that contamination of the target by an eclipsing system would rarely mimic the observed shape of the transits. The final validation of the planet is provided by 16 radial velocities obtained with HIRES on Keck 1 over a one year span. Although the velocities do not lead to a reliable orbit and mass determination, they are able to constrain the mass to a 3{\sigma} upper limit of 124 MEarth, safely in the regime of planetary masses, thus earning the designation Kepler-22b. The radiative equilibrium temperature is 262K for a planet in Kepler-22b's orbit. Although there is no evidence that Kepler-22b is a rocky planet, it is the first confirmed planet with a measured radius to orbit in the Habitable Zone of any star other than the Sun.

We report the detection of three transiting planets around a Sunlike star, which we designate Kepler-18. The transit signals were detected in photometric data from the Kepler satellite, and were confirmed to arise from planets using a combination of large transit-timing variations, radial-velocity variations, Warm-Spitzer observations, and statistical analysis of false-positive probabilities. The Kepler-18 star has a mass of 0.97M_sun, radius 1.1R_sun, effective temperature 5345K, and iron abundance [Fe/H]= +0.19. The planets have orbital periods of approximately 3.5, 7.6 and 14.9 days. The innermost planet "b" is a "super-Earth" with mass 6.9 \pm 3.4M_earth, radius 2.00 \pm 0.10R_earth, and mean density 4.9 \pm 2.4 g cm^-3. The two outer planets "c" and "d" are both low-density Neptune-mass planets. Kepler-18c has a mass of 17.3 \pm 1.9M_earth, radius 5.49 \pm 0.26R_earth, and mean density 0.59 \pm 0.07 g cm^-3, while Kepler-18d has a mass of 16.4 \pm 1.4M_earth, radius 6.98 \pm 0.33R_earth, and mean density 0.27 \pm 0.03 g cm^-3. Kepler-18c and Kepler-18d have orbital periods near a 2:1 mean-motion resonance, leading to large and readily detected transit timing variations.

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 the discovery of the Kepler-19 planetary system, which we first identified from a 9.3-day periodic transit signal in the Kepler photometry. From high-resolution spectroscopy of the star, we find a stellar effective temperature Teff=5541 \pm 60 K, a metallicity [Fe/H]=-0.13 \pm 0.06, and a surface gravity log(g)=4.59 \pm 0.10. We combine the estimate of Teff and [Fe/H] with an estimate of the stellar density derived from the photometric light curve to deduce a stellar mass of M_star = 0.936 \pm 0.040 M_Sun and a stellar radius of R_star = 0.850 \pm 0.018 R_Sun. We rule out the possibility that the transits result from an astrophysical false positive by first identifying the subset of stellar blends that reproduce the precise shape of the light curve. We conclude that the planetary scenario is more than three orders of magnitude more likely than a blend. The blend scenario is independently disfavored by the achromaticity of the transit: we measure a transit depth with Spitzer at 4.5 {\mu}m of 547+113-110 ppm, consistent with the depth measured in the Kepler optical bandpass of 567\pm6 ppm. We determine a physical radius of the planet Kepler-19b of R_p = 2.209 \pm 0.048 R_Earth. From radial-velocity observations of the star, we find an upper limit on the planet mass of 20.3 M_Earth, corresponding to a maximum density of 10.4 g cm^-3. We report a significant sinusoidal deviation of the transit times from a predicted linear ephemeris, which we conclude is due to an additional perturbing body in the system. We cannot uniquely determine the orbital parameters of the perturber, as various dynamical mechanisms match the amplitude, period, and shape of the transit timing signal and satisfy the host star's radial velocity limits. However, the perturber in these mechanisms has period <160 days and mass <6 M_Jup, confirming its planetary nature as Kepler-19c. [Abridged]

We present the faintest mid ecliptic latitude survey in the second part of HST archival search for outer Solar System bodies. We report the discovery of 28 new trans-Neptunian objects and 1 small centaur (R ~ 2km) in the band 5{\circ} - 20{\circ} off the ecliptic. The inclination distribution of these excited ob jects is consistent with the distribution derived from brighter ecliptic surveys. We suggest that the size and inclination distribution should be estimated consistently using suitable surveys with calibrated search algorithms and reliable orbital information.

Although OGLE-TR-56b was the second transiting exoplanet discovered, only one light curve, observed in 2006, has been published besides the discovery data. We present twenty-one light curves of nineteen different transits observed between July 2003 and July 2009 with the Magellan Telescopes and Gemini South. The combined analysis of the new light curves confirms a slightly inflated planetary radius relative to model predictions, with R_p = 1.378 +/- 0.090 R_J. However, the values found for the transit duration, semimajor axis, and inclination values differ significantly from the previous result, likely due to systematic errors. The new semimajor axis and inclination, a = 0.01942 +/- 0.00015 AU and i = 73.72 +/- 0.18 degrees, are smaller than previously reported, while the total duration, T_14 = 7931 +/- 38 s, is 18 minutes longer. The transit midtimes have errors from 23 s to several minutes, and no evidence is seen for transit midtime or duration variations. Similarly, no change is seen in the orbital period, implying a nominal stellar tidal decay factor of Q_* = 10^7, with a three-sigma lower limit of 10^5.7.

About one-third of the ~1200 transiting planet candidates detected in the first four months of \ik data are members of multiple candidate systems. There are 115 target stars with two candidate transiting planets, 45 with three, 8 with four, and one each with five and six. We characterize the dynamical properties of these candidate multi-planet systems. The distribution of observed period ratios shows that the vast majority of candidate pairs are neither in nor near low-order mean motion resonances. Nonetheless, there are small but statistically significant excesses of candidate pairs both in resonance and spaced slightly too far apart to be in resonance, particularly near the 2:1 resonance. We find that virtually all candidate systems are stable, as tested by numerical integrations that assume a nominal mass-radius relationship. Several considerations strongly suggest that the vast majority of these multi-candidate systems are true planetary systems. Using the observed multiplicity frequencies, we find that a single population of planetary systems that matches the higher multiplicities underpredicts the number of singly-transiting systems. We provide constraints on the true multiplicity and mutual inclination distribution of the multi-candidate systems, revealing a population of systems with multiple super-Earth-size and Neptune-size planets with low to moderate mutual inclinations.

We report the discovery of Kepler-15b, a new transiting exoplanet detected by NASA's Kepler mission. The transit signal with a period of 4.94 days was detected in the quarter 1 (Q1) Kepler photometry. For the first time, we have used the High-Resolution-Spectrograph (HRS) at the Hobby-Eberly Telescope (HET) to determine the mass of a Kepler planet via precise radial velocity (RV) measurements. The 24 HET/HRS radial velocities (RV) and 6 additional measurements from the FIES spectrograph at the Nordic Optical Telescope (NOT) reveal a Doppler signal with the same period and phase as the transit ephemeris. We used one HET/HRS spectrum of Kepler-15 taken without the iodine cell to determine accurate stellar parameters. The host star is a metal-rich ([Fe/H]=0.36+/-0.07) G-type main sequence star with T_eff=5515+/-124 K. The amplitude of the RV-orbit yields a mass of the planet of 0.66+/-0.1 M_Jup. The planet has a radius of 0.96+/-0.06 R_Jup and a mean bulk density of 0.9+/-0.2 g/cm^3. The planetary radius resides on the lower envelope for transiting planets with similar mass and irradiation level. This suggests significant enrichment of the planet with heavy elements. We estimate a heavy element mass of 30-40 M_Earth within Kepler-15b.

We have detected transits of the innermost planet "e" orbiting 55 Cnc (V=6.0), based on two weeks of nearly continuous photometric monitoring with the MOST space telescope. The transits occur with the period (0.74 d) and phase that had been predicted by Dawson & Fabrycky, and with the expected duration and depth for the crossing of a Sun-like star by a hot super-Earth. Assuming the star's mass and radius to be 0.963_{-0.029}^{+0.051} M_sun and 0.943 +/- 0.010 R_sun, the planet's mass, radius, and mean density are 8.63 +/- 0.35 Mearth, 2.00 +/- 0.14 Rearth, and 5.9_{-1.1}^{+1.5} g/cm^3. The mean density is comparable to that of Earth, despite the greater mass and consequently greater compression of the interior of 55 Cnc e. This suggests a rock-iron composition supplemented by a significant mass of water, gas, or other light elements. Outside of transits, we detected a sinusoidal signal resembling the expected signal due to the changing illuminated phase of the planet, but with a full range (168 +/- 70 ppm) too large to be reflected light or thermal emission. This signal has no straightforward interpretation and should be checked with further observations. The host star of 55 Cnc e is brighter than that of any other known transiting planet, which will facilitate future investigations.

The architectures of multiple planet systems can provide valuable constraints on models of planet formation, including orbital migration, and excitation of orbital eccentricities and inclinations. NASA's Kepler mission has identified 1235 transiting planet candidates (Borcuki et al 2011). The method of transit timing variations (TTVs) has already confirmed 7 planets in two planetary systems (Holman et al. 2010; Lissauer et al. 2011a). We perform a transit timing analysis of the Kepler planet candidates. We find that at least ~12% of planet candidates currently suitable for TTV analysis show evidence suggestive of TTVs, representing at least ~65 TTV candidates. In all cases, the time span of observations must increase for TTVs to provide strong constraints on planet masses and/or orbits, as expected based on n-body integrations of multiple transiting planet candidate systems (assuming circular and coplanar orbits). We find that the fraction of planet candidates showing TTVs in this data set does not vary significantly with the number of transiting planet candidates per star, suggesting significant mutual inclinations and that many stars with a single transiting planet should host additional non-transiting planets. We anticipate that Kepler could confirm (or reject) at least ~12 systems with multiple transiting planet candidates via TTVs. Thus, TTVs will provide a powerful tool for confirming transiting planets and characterizing the orbital dynamics of low-mass planets. If Kepler observations were extended to at least six years, then TTVs would provide much more precise constraints on the dynamics of systems with multiple transiting planets and would become sensitive to planets with orbital periods extending into the habitable zone of solar-type stars.

We present the discovery of a hot Jupiter transiting an F star in a close visual (0.3" sky projected angular separation) binary system. The dilution of the host star's light by the nearly equal magnitude stellar companion (~ 0.5 magnitudes fainter) significantly affects the derived planetary parameters, and if left uncorrected, leads to an underestimate of the radius and mass of the planet by 10% and 60%, respectively. Other published exoplanets, which have not been observed with high-resolution imaging, could similarly have unresolved stellar companions and thus have incorrectly derived planetary parameters. Kepler-14b (KOI-98) has a period of P = 6.790 days and correcting for the dilution, has a mass of Mp = 8.40 +0.19-0.18 MJ and a radius of Rp = 1.136 +0.073-0.054 RJ, yielding a mean density of rho = 7.1 +- 1.1 g cm-3.

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 has recently announced the discovery of Kepler-10 b, the smallest exoplanet discovered to date and the first rocky planet found by the spacecraft. A second, 45-day period transit-like signal present in the photometry from the first eight months of data could not be confirmed as being caused by a planet at the time of that announcement. Here we apply the light-curve modeling technique known as BLENDER to explore the possibility that the signal might be due to an astrophysical false positive (blend). To aid in this analysis we report the observation of two transits with the Spitzer Space Telescope at 4.5 {\mu}m. When combined they yield a transit depth of 344 \pm 85 ppm that is consistent with the depth in the Kepler passband (376 \pm 9 ppm, ignoring limb darkening), which rules out blends with an eclipsing binary of a significantly different color than the target. Using these observations along with other constraints from high resolution imaging and spectroscopy we are able to exclude the vast majority of possible false positives. We assess the likelihood of the remaining blends, and arrive conservatively at a false alarm rate of 1.6 \times 10-5 that is small enough to validate the candidate as a planet (designated Kepler-10 c) with a very high level of confidence. The radius of this object is measured to be Rp = 2.227+0.052 -0.057 Earth radii. Kepler-10 c represents another example (with Kepler-9 d and Kepler-11 g) of statistical "validation" of a transiting exoplanet, as opposed to the usual "confirmation" that can take place when the Doppler signal is detected or transit timing variations are measured. It is anticipated that many of Kepler's smaller candidates will receive a similar treatment since dynamical confirmation may be difficult or impractical with the sensitivity of current instrumentation.

We present transit photometry of three exoplanets, TrES-4b, HAT-P-3b, and WASP-12b, allowing for refined estimates of the systems' parameters. TrES-4b and WASP-12b were confirmed to be "bloated" planets, with radii of 1.706 +/- 0.056 R_Jup and 1.736 +/- 0.092 R_Jup, respectively. These planets are too large to be explained with standard models of gas giant planets. In contrast, HAT-P-3b has a radius of 0.827 +/- 0.055 R_Jup, smaller than a pure hydrogen-helium planet and indicative of a highly metal-enriched composition. Analyses of the transit timings revealed no significant departures from strict periodicity. For TrES-4, our relatively recent observations allow for improvement in the orbital ephemerides, which is useful for planning future observations.

We present photometry of 4 transits of the exoplanet WASP-4b, each with a precision of approximately 500 ppm and a time sampling of 40-60s. We have used the data to refine the estimates of the system parameters and ephemerides. During two of the transits we observed a short-lived, low-amplitude anomaly that we interpret as the occultation of a starspot by the planet. We also find evidence for a pair of similar anomalies in previously published photometry. The recurrence of these anomalies suggests that the stellar rotation axis is nearly aligned with the orbital axis, or else the star spot would not have remained on the transit chord. By analyzing the timings of the anomalies we find the sky-projected stellar obliquity to be -1_{-12}^{+14} degrees. This result is consistent with (and more constraining than) a recent observation of the Rossiter-McLaughlin effect. It suggests that the planet migration mechanism preserved the initially low obliquity, or else that tidal evolution has realigned the system. Future applications of this method using data from the Corot and Kepler missions will allow spin-orbit alignment to be probed for many other exoplanets.

In this letter we present an overview of the rich population of systems with multiple candidate transiting planets found in the first four months of Kepler data. The census of multiples includes 115 targets that show 2 candidate planets, 45 with 3, 8 with 4, and 1 each with 5 and 6, for a total of 170 systems with 408 candidates. When compared to the 827 systems with only one candidate, the multiples account for 17 percent of the total number of systems, and a third of all the planet candidates. We compare the characteristics of candidates found in multiples with those found in singles. False positives due to eclipsing binaries are much less common for the multiples, as expected. Singles and multiples are both dominated by planets smaller than Neptune; 69 +2/-3 percent for singles and 86 +2/-5 percent for multiples. This result, that systems with multiple transiting planets are less likely to include a transiting giant planet, suggests that close-in giant planets tend to disrupt the orbital inclinations of small planets in flat systems, or maybe even to prevent the formation of such systems in the first place.

We report the distribution of planets as a function of planet radius (R_p), orbital period (P), and stellar effective temperature (Teff) for P < 50 day orbits around GK stars. These results are based on the 1,235 planets (formally "planet candidates") from the Kepler mission that include a nearly complete set of detected planets as small as 2 Earth radii (Re). For each of the 156,000 target stars we assess the detectability of planets as a function of R_p and P. We also correct for the geometric probability of transit, R*/a. We consider first stars within the "solar subset" having Teff = 4100-6100 K, logg = 4.0-4.9, and Kepler magnitude Kp < 15 mag. We include only those stars having noise low enough to permit detection of planets down to 2 Re. We count planets in small domains of R_p and P and divide by the included target stars to calculate planet occurrence in each domain. Occurrence of planets varies by more than three orders of magnitude and increases substantially down to the smallest radius (2 Re) and out to the longest orbital period (50 days, ~0.25 AU) in our study. For P < 50 days, the radius distribution is given by a power law, df/dlogR= k R^\alpha. This rapid increase in planet occurrence with decreasing planet size agrees with core-accretion, but disagrees with population synthesis models. We fit occurrence as a function of P to a power law model with an exponential cutoff below a critical period P_0. For smaller planets, P_0 has larger values, suggesting that the "parking distance" for migrating planets moves outward with decreasing planet size. We also measured planet occurrence over Teff = 3600-7100 K, spanning M0 to F2 dwarfs. The occurrence of 2-4 Re planets in the Kepler field increases with decreasing Teff, making these small planets seven times more abundant around cool stars than the hottest stars in our sample. [abridged]

On 1 February 2011 the Kepler Mission released data for 156,453 stars observed from the beginning of the science observations on 2 May through 16 September 2009. There are 1235 planetary candidates with transit like signatures detected in this period. These are associated with 997 host stars. Distributions of the characteristics of the planetary candidates are separated into five class-sizes; 68 candidates of approximately Earth-size (radius < 1.25 Earth radii), 288 super-Earth size (1.25 Earth radii < radius < 2 Earth radii), 662 Neptune-size (2 Earth radii < radius < 6 Earth radii), 165 Jupiter-size (6 Earth radii < radius < 15 Earth radii), and 19 up to twice the size of Jupiter (15 Earth radii < radius < 22 Earth radii). In the temperature range appropriate for the habitable zone, 54 candidates are found with sizes ranging from Earth-size to larger than that of Jupiter. Five are less than twice the size of the Earth. Over 74% of the planetary candidates are smaller than Neptune. The observed number versus size distribution of planetary candidates increases to a peak at two to three times Earth-size and then declines inversely proportional to area of the candidate. Our current best estimates of the intrinsic frequencies of planetary candidates, after correcting for geometric and sensitivity biases, are 6% for Earth-size candidates, 7% for super-Earth size candidates, 17% for Neptune-size candidates, and 4% for Jupiter-size candidates. Multi-candidate, transiting systems are frequent; 17% of the host stars have multi-candidate systems, and 33.9% of all the candidates are part of multi-candidate systems.

We present time series photometry and constraints on additional planets in five of the exoplanetary systems studied by the EPOCh (Extrasolar Planet Observation and Characterization) component of the NASA EPOXI mission: HAT-P-4, TrES-3, TrES-2, WASP-3, and HAT-P-7. We conduct a search of the high-precision time series for photometric transits of additional planets. We find no candidate transits with significance higher than our detection limit. From Monte Carlo tests of the time series using putative periods from 0.5 days to 7 days, we demonstrate the sensitivity to detect Neptune-sized companions around TrES-2, sub-Saturn-sized companions in the HAT-P-4, TrES-3, and WASP-3 systems, and Saturn-sized companions around HAT-P-7. We investigate in particular our sensitivity to additional transits in the dynamically favorable 3:2 and 2:1 exterior resonances with the known exoplanets: if we assume coplanar orbits with the known planets, then companions in these resonances with HAT-P-4b, WASP-3b, and HAT-P-7b would be expected to transit, and we can set lower limits on the radii of companions in these systems. In the nearly grazing exoplanetary systems TrES-3 and TrES-2, additional coplanar planets in these resonances are not expected to transit. However, we place lower limits on the radii of companions that would transit if the orbits were misaligned by 2.0 degrees and 1.4 degrees for TrES-3 and TrES-2, respectively.

NASA's Kepler Mission uses transit photometry to determine the frequency of earth-size planets in or near the habitable zone of Sun-like stars. The mission reached a milestone toward meeting that goal: the discovery of its first rocky planet, Kepler-10b. Two distinct sets of transit events were detected: 1) a 152 +/- 4 ppm dimming lasting 1.811 +/- 0.024 hours with ephemeris T[BJD]=2454964.57375+N*0.837495 days and 2) a 376 +/- 9 ppm dimming lasting 6.86 +/- 0.07 hours with ephemeris T[BJD]=2454971.6761+N*45.29485 days. Statistical tests on the photometric and pixel flux time series established the viability of the planet candidates triggering ground-based follow-up observations. Forty precision Doppler measurements were used to confirm that the short-period transit event is due to a planetary companion. The parent star is bright enough for asteroseismic analysis. Photometry was collected at 1-minute cadence for >4 months from which we detected 19 distinct pulsation frequencies. Modeling the frequencies resulted in precise knowledge of the fundamental stellar properties. Kepler-10 is a relatively old (11.9 +/- 4.5 Gyr) but otherwise Sun-like Main Sequence star with Teff=5627 +/- 44 K, Mstar=0.895 +/- 0.060 Msun, and Rstar=1.056 +/- 0.021 Rsun. Physical models simultaneously fit to the transit light curves and the precision Doppler measurements yielded tight constraints on the properties of Kepler-10b that speak to its rocky composition: Mpl=4.56 +/- 1.29 Mearth, Rpl=1.416 +/- 0.036 Rearth, and density=8.8 +/- 2.9 gcc. Kepler-10b is the smallest transiting exoplanet discovered to date.

The Kepler spacecraft has been monitoring the light from 150,000 stars in its primary quest to detect transiting exoplanets. Here we report on the detection of an eclipsing stellar hierarchical triple, identified in the Kepler photometry. KOI-126 (A,(B, C)), is composed of a low-mass binary (masses M_B = 0.2413+/-0.0030 M_Sun, M_C = 0.2127+/-0.0026 M_Sun; radii R_B = 0.2543+/-0.0014 R_Sun, R_C = 0.2318+/-0.0013 R_Sun; orbital period P_1 = 1.76713+/-0.00019 days) on an eccentric orbit about a third star (mass M_A = 1.347+/-0.032 M_Sun; radius R_A = 2.0254+/-0.0098 R_Sun; period of orbit around the low-mass binary P_2 = 33.9214+/-0.0013 days; eccentricity of that orbit e_2 = 0.3043+/-0.0024). The low-mass pair probe the poorly sampled fully-convective stellar domain offering a crucial benchmark for theoretical stellar models.

Doppler planet searches have discovered that giant planets follow orbits with a wide range of orbital eccentricities, revolutionizing theories of planet formation. The discovery of hundreds of exoplanet candidates by NASA's Kepler mission enables astronomers to characterize the eccentricity distribution of small exoplanets. Measuring the eccentricity of individual planets is only practical in favorable cases that are amenable to complementary techniques (e.g., radial velocities, transit timing variations, occultation photometry). Yet even in the absence of individual eccentricities, it is possible to study the distribution of eccentricities based on the distribution of transit durations (relative to the maximum transit duration for a circular orbit). We analyze the transit duration distribution of Kepler planet candidates. We find that for host stars with T_eff > 5100 K we cannot invert this to infer the eccentricity distribution at this time due to uncertainties and possible systematics in the host star densities. With this limitation in mind, we compare the observed transit duration distribution with models to rule out extreme distributions. If we assume a Rayleigh eccentricity distribution for Kepler planet candidates, then we find best-fits with a mean eccentricity of 0.1-0.25 for host stars with T_eff < 5100 K. We compare the transit duration distribution for different subsets of Kepler planet candidates and discuss tentative trends with planetary radius and multiplicity. High-precision spectroscopic follow-up observations for a large sample of host stars will be required to confirm which trends are real and which are the results of systematic errors in stellar radii. Finally, we identify planet candidates that must be eccentric or have a significantly underestimated stellar radius.

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.

We present optical photometry of 16 transits of the super-Earth GJ 1214b, allowing us to refine the system parameters and search for additional planets via transit timing. Starspot-crossing events are detected in two light curves, and the star is found to be variable by a few percent. Hence, in our analysis, special attention is given to systematic errors that result from star spots. The planet-to-star radius ratio is 0.11610+/-0.00048, subject to a possible upward bias by a few percent due to the unknown spot coverage. Even assuming this bias to be negligible, the mean density of planet can be either 3.03+/-0.50 g cm^{-3} or 1.89+/-0.33 g cm^{-3}, depending on whether the stellar radius is estimated from evolutionary models or from an empirical mass-luminosity relation combined with the light curve parameters. One possible resolution is that the orbit is eccentric (e approximately equal to 0.14), which would favor the higher density, and hence a much thinner atmosphere for the planet. The transit times were found to be periodic within about 15s, ruling out the existence of any other super-Earths with periods within a factor-of-two of the known planet.

We present observations of the Rossiter-McLaughlin effect for two exoplanetary systems, revealing the orientations of their orbits relative to the rotation axes of their parent stars. HAT-P-4b is prograde, with a sky-projected spin-orbit angle of lambda = -4.9 +/- 11.9 degrees. In contrast, HAT-P-14b is retrograde, with lambda = 189.1 +/- 5.1 degrees. These results conform with a previously noted pattern among the stellar hosts of close-in giant planets: hotter stars have a wide range of obliquities and cooler stars have low obliquities. This, in turn, suggests that three-body dynamics and tidal dissipation are responsible for the short-period orbits of many exoplanets. In addition, our data revealed a third body in the HAT-P-4 system, which could be a second planet or a companion star.

Light curves from the Kepler Mission contain valuable information on the nature of the phenomena producing the transit-like signals. To assist in exploring the possibility that they are due to an astrophysical false positive, we describe a procedure (BLENDER) to model the photometry in terms of a "blend" rather than a planet orbiting a star. A blend may consist of a background or foreground eclipsing binary (or star-planet pair) whose eclipses are attenuated by the light of the candidate and possibly other stars within the photometric aperture. We apply BLENDER to the case of Kepler-9, a target harboring two previously confirmed Saturn-size planets (Kepler-9b and Kepler-9c) showing transit timing variations, and an additional shallower signal with a 1.59-day period suggesting the presence of a super-Earth-size planet. Using BLENDER together with constraints from other follow-up observations we are able to rule out all blends for the two deeper signals, and provide independent validation of their planetary nature. For the shallower signal we rule out a large fraction of the false positives that might mimic the transits. The false alarm rate for remaining blends depends in part (and inversely) on the unknown frequency of small-size planets. Based on several realistic estimates of this frequency we conclude with very high confidence that this small signal is due to a super-Earth-size planet (Kepler-9d) in a multiple system, rather than a false positive. The radius is determined to be 1.64 (+0.19/-0.14) R(Earth), and current spectroscopic observations are as yet insufficient to establish its mass.

As part of the NASA EPOXI Mission of Opportunity, we observed seven known transiting extrasolar planet systems in order to construct time series photometry of extremely high phase coverage and precision. Here we present the results for four "hot-Jupiter systems" with near-solar stars - HAT-P-4, TrES-3, TrES-2 and WASP-3. We observe ten transits of HAT-P-4, estimating the planet radius Rp = 1.332 \pm 0.052 RJup, the stellar radius R \star = 1.602 \pm 0.061 R \odot, the inclination i = 89.67 \pm 0.30 degrees and the transit duration from first to fourth contact T = 255.6 \pm 1.9 minutes. For TrES-3, we observe seven transits, and find Rp = 1.320 \pm 0.057 RJup, R\star = 0.817 \pm 0.022 R\odot, i = 81.99 \pm 0.30 degrees and T = 81.9 \pm 1.1 minutes. We also note a long term variability in the TrES-3 light curve, which may be due to star spots. We observe nine transits of TrES-2, and find Rp = 1.169 \pm 0.034 RJup, R\star = 0.940 \pm 0.026 R\odot, i = 84.15 \pm 0.16 degrees and T = 107.3 \pm 1.1 minutes. Finally we observe eight transits of WASP-3, finding Rp = 1.385 \pm 0.060 RJup, R\star = 1.354 \pm 0.056 R\odot, i = 84.22 \pm 0.81 degrees and T = 167.3 \pm 1.3 minutes. We present refined orbital periods and times of transit for each target. We state 95% confidence upper limits on the secondary eclipse depths in our broadband visible bandpass centered on 650 nm. These limits are 0.073% for HAT-P-4, 0.062% for TrES-3, 0.16% for TrES-2 and 0.11% for WASP-3. We combine the TrES-3 secondary eclipse information with the existing published data and confirm that the atmosphere likely does not have a temperature inversion.

We present and discuss five candidate exoplanetary systems identified with the Kepler spacecraft. These five systems show transits from multiple exoplanet candidates. Should these objects prove to be planetary in nature, then these five systems open new opportunities for the field of exoplanets and provide new insights into the formation and dynamical evolution of planetary systems. We discuss the methods used to identify multiple transiting objects from the Kepler photometry as well as the false-positive rejection methods that have been applied to these data. One system shows transits from three distinct objects while the remaining four systems show transits from two objects. Three systems have planet candidates that are near mean motion commensurabilities---two near 2:1 and one just outside 5:2. We discuss the implications that multitransiting systems have on the distribution of orbital inclinations in planetary systems, and hence their dynamical histories; as well as their likely masses and chemical compositions. A Monte Carlo study indicates that, with additional data, most of these systems should exhibit detectable transit timing variations (TTV) due to gravitational interactions---though none are apparent in these data. We also discuss new challenges that arise in TTV analyses due to the presence of more than two planets in a system.

We present observations of three distinct transits of HD 17156b obtained with the Fine Guidance Sensors (FGS) on board the Hubble Space Telescope} (HST). We analyzed both the transit photometry and previously published radial velocities to find the planet-star radius ratio R_p/R_s = 0.07454 +/- 0.00035, inclination i=86.49 +0.24/-0.20 deg, and scaled semi-major axis a/R = 23.19 +0.32/-0.27. This last value translates directly to a mean stellar density determination of 0.522 +0.021/-0.018 g cm^-3. Analysis of asteroseismology observations by the companion paper of Gilliland et al. (2009) provides a consistent but significantly refined measurement of the stellar mean density. We compare stellar isochrones to this density estimate and find M_s = 1.275 +/- 0.018 M_sun and a stellar age of $3.37 +0.20/-0.47 Gyr. Using this estimate of M_s and incorporating the density constraint from asteroseismology, we model both the photometry and published radial velocities to estimate the planet radius R_p= 1.0870 +/- 0.0066 Jupiter radii and the stellar radius R_s = 1.5007 +/- 0.0076 R_sun. The planet radius is larger than that found in previous studies and consistent with theoretical models of a solar-composition gas giant of the same mass and equilibrium temperature. For the three transits, we determine the times of mid-transit to a precision of 6.2 s, 7.6 s, and 6.9 s, and the transit times for HD 17156 do not show any significant departures from a constant period. The joint analysis of transit photometry and asteroseismology presages similar studies that will be enabled by the NASA Kepler Mission.

We present ground-based optical observations of the September 2009 and January 2010 transits of HD 80606b. Based on 3 partial light curves of the September 2009 event, we derive a midtransit time of T_c [HJD] = 2455099.196 +- 0.026, which is about 1 sigma away from the previously predicted time. We observed the January 2010 event from 9 different locations, with most phases of the transit being observed by at least 3 different teams. We determine a midtransit time of T_c [HJD] = 2455210.6502 +- 0.0064, which is within 1.3 sigma of the time derived from a Spitzer observation of the same event.

Bernstein et al. (2004) found that the population of faint (R>26) trans-Neptunian objects (TNOs) known at that time was dominated by "Classical" objects, which have low inclinations (i<5 degrees) and distances 40--45 AU. Since those observations, the number of faint TNOs whose orbits are sufficiently well known to be classified as "Classical" or "Excited" has grown from seven to 39. We analyze the dynamical classifications of faint TNOs known today and find that this population is dominated by Excited objects. We discuss some implications of this result.

We discovered evidence for a possible additional 0.75 R_Earth transiting planet in the NASA EPOXI observations of the known M dwarf exoplanetary system GJ 436. Based on an ephemeris determined from the EPOXI data, we predicted a transit event in an extant Spitzer Space Telescope 8 micron data set of this star. Our subsequent analysis of those Spitzer data confirmed the signal of the predicted depth and at the predicted time, but we found that the transit depth was dependent on the aperture used to perform the photometry. Based on these suggestive findings, we gathered new Warm Spitzer Observations of GJ 436 at 4.5 microns spanning a time of transit predicted from the EPOXI and Spitzer 8 micron candidate events. The 4.5 micron data permit us to rule out a transit at high confidence, and we conclude that the earlier candidate transit signals resulted from correlated noise in the EPOXI and Spitzer 8 micron observations. In the course of this investigation, we developed a novel method for correcting the intrapixel sensitivity variations of the 3.6 and 4.5 micron channels of the Infrared Array Camera (IRAC) instrument. We demonstrate the sensitivity of Warm Spitzer observations of M dwarfs to confirm sub-Earth sized planets. Our analysis will inform similar work that will be undertaken to use Warm Spitzer observations to confirm rocky planets discovered by the Kepler mission.

We report photometric and radial velocity observations of the XO-4 transiting planetary system, conducted with the FLWO 1.2m telescope and the 8.2m Subaru Telescope. Based on the new light curves, the refined transit ephemeris of XO-4b is $P = 4.1250828 \pm 0.0000040$ days and $T_c [BJD_TDB] = 2454485.93323 \pm 0.00039$. We measured the Rossiter-McLaughlin effect of XO-4b and estimated the sky-projected angle between the stellar spin axis and the planetary orbital axis to be $\lambda = -46.7^{\circ} ^{+8.1^{\circ}}_{-6.1^{\circ}}$. This measurement of $\lambda$ is less robust than in some other cases because the impact parameter of the transit is small, causing a strong degeneracy between $\lambda$ and the projected stellar rotational velocity. Nevertheless, our finding of a spin-orbit misalignment suggests that the migration process for XO-4b involved few-body dynamics rather than interaction with a gaseous disk. In addition, our result conforms with the pattern reported by Winn et al. (2010, ApJL, 718, L145) that high obliquities are preferentially found for stars with effective temperatures hotter than 6250~K.

We introduce a novel search technique that can identify trans-neptunian objects in three to five exposures of a pointing within a single Hubble Space Telescope orbit. The process is fast enough to allow the discovery of candidates soon after the data are available. This allows sufficient time to schedule follow up observations with HST within a month. We report the discovery of 14 slow-moving objects found within 5\circ of the ecliptic in archival data taken with the Wide Field Channel of the Advanced Camera for Surveys. The luminosity function of these objects is consistent with previous ground-based and space-based results. We show evidence that the size distribution of both high and low inclination populations is similar for objects smaller than 100 km, as expected from collisional evolution models, while their size distribution differ for brighter objects. We suggest the two populations formed in different parts of the protoplanetary disk and after being dynamically mixed have collisionally evolved together. Among the objects discovered there is an equal mass binary with an angular separation ~ 0."53.

Among other things, studies of the formation and evolution of planetary systems currently draw on two important observational resources: the precise characterization available for planets that transit their parent stars and the frequency and nature of systems with multiple planets. Thus far, the study of transiting exoplanets has focused almost exclusively on systems with only one planet, except for considering the influence of additional planets on the transit light curve, mostly through transit timing variations (TTVs). This work considers systems where multiple planets are seen to transit the same star and concludes that such "multi-transiting" systems will be the most information-rich planetary systems besides our own solar system. Five new candidate multi-transiting systems from \emph{Kepler} have been announced in Steffen et al. 2010, though these candidates have not yet been fully confirmed as planets. In anticipation of the likely confirmation of multi-transiting systems, we discuss the value of these systems in detail. For example, proper interpretation of transit timing variations is significantly improved in multi-transiting systems. The true mutual inclination, a valuable probe of planetary formation, can also be well determined in certain systems, especially through Rossiter-McLaughlin measurements of each planet. In addition, such systems may undergo predictable and observable mutual events, where one planet crosses over the other, which allow for unique constraints on various physical and orbital parameters, particularly the mutual inclination.

We report the discovery of an eclipsing companion to NLTT 41135, a nearby M5 dwarf that was already known to have a wider, slightly more massive common proper motion companion, NLTT 41136, at 2.4 arcsec separation. Analysis of combined-light and radial velocity curves of the system indicates that NLTT 41135B is a 31-34 +/- 3 MJup brown dwarf (where the range depends on the unknown metallicity of the host star) on a circular orbit. The visual M-dwarf pair appears to be physically bound, so the system forms a hierarchical triple, with masses approximately in the ratio 8:6:1. The eclipses are grazing, preventing an unambiguous measurement of the secondary radius, but follow-up observations of the secondary eclipse (e.g. with the James Webb Space Telescope) could permit measurements of the surface brightness ratio between the two objects, and thus place constraints on models of brown dwarfs.

We present time series photometry of the M dwarf transiting exoplanet system GJ 436 obtained with the the EPOCh (Extrasolar Planet Observation and Characterization) component of the NASA EPOXI mission. We conduct a search of the high-precision time series for additional planets around GJ 436, which could be revealed either directly through their photometric transits, or indirectly through the variations these second planets induce on the transits of the previously known planet. In the case of GJ 436, the presence of a second planet is perhaps indicated by the residual orbital eccentricity of the known hot Neptune companion. We find no candidate transits with significance higher than our detection limit. From Monte Carlo tests of the time series, we rule out transiting planets larger than 1.5 R_Earth interior to GJ 436b with 95% confidence, and larger than 1.25 R_Earth with 80% confidence. Assuming coplanarity of additional planets with the orbit of GJ 436b, we cannot expect that putative planets with orbital periods longer than about 3.4 days will transit. However, if such a planet were to transit, we rule out planets larger than 2.0 R_Earth with orbital periods less than 8.5 days with 95% confidence. We also place dynamical constraints on additional bodies in the GJ 436 system. Our analysis should serve as a useful guide for similar analyses for which radial velocity measurements are not available, such as those discovered by the Kepler mission. These dynamical constraints on additional planets with periods from 0.5 to 9 days rule out coplanar secular perturbers as small as 10 M_Earth and non-coplanar secular perturbers as small as 1 M_Earth in orbits close in to GJ 436b. We present refined estimates of the system parameters for GJ 436. We also report a sinusoidal modulation in the GJ 436 light curve that we attribute to star spots. [Abridged]

We present new spectroscopic and photometric observations of the transiting exoplanetary system WASP-3. Spectra obtained during two separate transits exhibit the Rossiter-McLaughlin (RM) effect and allow us to estimate the sky-projected angle between the planetary orbital axis and the stellar rotation axis, lambda = 3.3^{+2.5}_{-4.4} degrees. This alignment between the axes suggests that WASP-3b has a low orbital inclination relative to the equatorial plane of its parent star. During our first night of spectroscopic measurements, we observed an unexpected redshift briefly exceeding the expected sum of the orbital and RM velocities by 140 m/s. This anomaly could represent the occultation of material erupting from the stellar photosphere, although it is more likely to be an artifact caused by moonlight scattered into the spectrograph.

We have analyzed the first 3.75 years of data from TAOS, the Taiwanese American Occultation Survey. TAOS monitors bright stars to search for occultations by Kuiper Belt Objects (KBOs). This dataset comprises 5e5 star-hours of multi-telescope photometric data taken at 4 or 5 Hz. No events consistent with KBO occultations were found in this dataset. We compute the number of events expected for the Kuiper Belt formation and evolution models of Pan & Sari (2005), Kenyon & Bromley (2004), Benavidez & Campo Bagatin (2009), and Fraser (2009). A comparison with the upper limits we derive from our data constrains the parameter space of these models. This is the first detailed comparison of models of the KBO size distribution with data from an occultation survey. Our results suggest that the KBO population is comprised of objects with low internal strength and that planetary migration played a role in the shaping of the size distribution.

A decade ago, the detection of the first transiting extrasolar planet provided a direct constraint on its composition and opened the door to spectroscopic investigations of extrasolar planetary atmospheres. As such characterization studies are feasible only for transiting systems that are both nearby and for which the planet-to-star radius ratio is relatively large, nearby small stars have been surveyed intensively. Doppler studies and microlensing have uncovered a population of planets with minimum masses of 1.9-10 times the Earth's mass (M_Earth), called super-Earths. The first constraint on the bulk composition of this novel class of planets was afforded by CoRoT-7b, but the distance and size of its star preclude atmospheric studies in the foreseeable future. Here we report observations of the transiting planet GJ 1214b, which has a mass of 6.55 M_Earth and a radius 2.68 times Earth's radius (R_Earth), indicating that it is intermediate in stature between Earth and the ice giants of the Solar System. We find that the planetary mass and radius are consistent with a composition of primarily water enshrouded by a hydrogen-helium envelope that is only 0.05% of the mass of the planet. The atmosphere is probably escaping hydrodynamically, indicating that it has undergone significant evolution during its history. As the star is small and only 13 parsecs away, the planetary atmosphere is amenable to study with current observatories.

The highly irradiated transiting exoplanet, HAT-P-7b, currently provides one of the best opportunities for studying planetary emission in the optical and infrared wavelengths. We observe six near-consecutive secondary eclipses of HAT-P-7b at optical wavelengths with the EPOXI spacecraft. We place an upper limit on the relative eclipse depth of 0.055% (95% confidence). We also analyze Spitzer observations of the same target in the infrared, obtaining secondary eclipse depths of 0.098+/-0.017%, 0.159+/-0.022%, 0.245+/-0.031% and 0.225+/-0.052% in the 3.6, 4.5, 5.8 and 8.0 micron IRAC bands respectively. We combine these measurements with the recently published Kepler secondary eclipse measurement, and generate atmospheric models for the day-side of the planet that are consistent with both the optical and infrared measurements. The data are best fit by models with a temperature inversion, as expected from the high incident flux. The models predict a low optical albedo of ~< 0.13, with subsolar abundances of Na, K, TiO and VO. We also find that the best fitting models predict that 10% of the absorbed stellar flux is redistributed to the night side of the planet, which is qualitatively consistent with the inefficient day-night redistribution apparent in the Kepler phase curve. Models without thermal inversions fit the data only at the 1.25 sigma level, and also require an overabundance of methane, which is not expected in the very hot atmosphere of HAT-P-7b. We also analyze the eight transits of HAT-P-7b present in the EPOXI dataset and improve the constraints on the system parameters, finding a period of P = 2.2047308+/-0.0000025 days, a stellar radius of R* = 1.824+/-0.089Rsun, a planetary radius of Rp = 1.342+/-0.068RJup and an inclination of i = 85.7+3.5-2.2 deg.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

We show that the exoplanet HAT-P-7b has an extremely tilted orbit, with a true angle of at least 86 degrees with respect to its parent star's equatorial plane, and a strong possibility of retrograde motion. We also report evidence for an additional planet or companion star. The evidence for the unparalleled orbit and the third body is based on precise observations of the star's apparent radial velocity. The anomalous radial velocity due to rotation (the Rossiter-McLaughlin effect) was found to be a blueshift during the first half of the transit and a redshift during the second half, an inversion of the usual pattern, implying that the angle between the sky-projected orbital and stellar angular momentum vectors is 182.5 +/- 9.4 degreees. The third body is implicated by excess radial-velocity variation of the host star over 2 yr. Some possible explanations for the tilted orbit are a close encounter with another planet, the Kozai effect, and resonant capture by an inward-migrating outer planet.

We present photometric and spectroscopic observations of the 2009 Feb. 2 transit of the exoplanet XO-3b. The new data show that the planetary orbital axis and stellar rotation axis are misaligned, as reported earlier by Hebrard and coworkers. We find the angle between the sky projections of the two axes to be 37.3 +/- 3.7 deg, as compared to the previously reported value of 70 +/- 15 deg. The significance of this discrepancy is unclear because there are indications of systematic effects. XO-3b is the first exoplanet known to have a highly inclined orbit relative to the equatorial plane of its parent star, and as such it may fulfill the predictions of some scenarios for the migration of massive planets into close-in orbits. We revisit the statistical analysis of spin-orbit alignment in hot-Jupiter systems. Assuming the stellar obliquities to be drawn from a single Rayleigh distribution, we find the mode of the distribution to be 13^{+5}_{-2} deg. However, it remains the case that a model representing two different migration channels--in which some planets are drawn from a perfectly-aligned distribution and the rest are drawn from an isotropic distribution--is favored over a single Rayleigh distribution.

We conducted a search for occultations of bright stars by Kuiper Belt Objects (KBOs) to estimate the density of sub-km KBOs in the sky. We report here the first results of this occultation survey of the outer solar system conducted in June 2007 and June/July 2008 at the MMT Observatory using Megacam, the large MMT optical imager. We used Megacam in a novel shutterless continuous--readout mode to achieve high precision photometry at 200 Hz. We present an analysis of 220 star hours at signal-to-noise ratio of 25 or greater. The survey efficiency is greater than 10% for occultations by KBOs of diameter d>=0.7 km, and we report no detections in our dataset. We set a new 95% confidence level upper limit for the surface density \Sigma_N(d) of KBOs larger than 1 km: \Sigma_N(d>=1 km) <= 2.0e8 deg^-2, and for KBOs larger than 0.7 km \Sigma_N(d>= 0.7 km) <= 4.8e8 deg^-2.

We present photometry of six transits of the exoplanet XO-2b. By combining the light-curve analysis with theoretical isochrones to determine the stellar properties, we find the planetary radius to be 0.996 +0.031/-0.018 rjup and the planetary mass to be 0.565 +/- 0.054 mjup. These results are consistent with those reported previously, and are also consistent with theoretical models for gas giant planets. The mid-transit times are accurate to within 1 min and are consistent with a constant period. However, the period we derive differs by 2.5 sigma from the previously published period. More data are needed to tell whether the period is actually variable (as it would be in the presence of an additional body) or if the timing errors have been underestimated.

The transiting exoplanet HD 149026b is an important case for theories of planet formation and planetary structure, because the planet's relatively small size has been interpreted as evidence for a highly metal-enriched composition. We present observations of 4 transits with the Near Infrared Camera and Multi-Object Spectrometer on the Hubble Space Telescope, within a wavelength range of 1.1--2.0 $\mu$m. Analysis of the light curve gives the most precise estimate yet of the stellar mean density, $\rho_\star = 0.497^{+0.042}_{-0.057}$ g cm$^{-3}$. By requiring agreement between the observed stellar properties (including $\rho_\star$) and stellar evolutionary models, we refine the estimate of the stellar radius: $R_\star = 1.541^{+0.046}_{-0.042}$ $R_\sun$. We also find a deeper transit than has been measured at optical and mid-infrared wavelengths. Taken together, these findings imply a planetary radius of $R_p = 0.813^{+0.027}_{-0.025}$ $R_{\rm Jup}$, which is larger than earlier estimates. Models of the planetary interior still require a metal-enriched composition, although the required degree of metal enrichment is reduced. It is also possible that the deeper NICMOS transit is caused by wavelength-dependent absorption by constituents in the planet's atmosphere, although simple model atmospheres do not predict this effect to be strong enough to account for the discrepancy. We use the 4 newly-measured transit times to compute a refined transit ephemeris.

We present photometry of two transits of the giant planet WASP-4b with a photometric precision of 400-800 parts per million and a time sampling of 25-40 seconds. The two midtransit times are determined to within 6 seconds. Together with previously published times, the data are consistent with a constant orbital period, giving no compelling evidence for period variations that would be produced by a satellite or additional planets. Analysis of the new photometry, in combination with stellar-evolutionary modeling, gives a planetary mass and radius of 1.237 +/- 0.064 M_jup and 1.365 +/- 0.021 R_jup. The planet is 15% larger than expected based on previously published models of solar-composition giant planets. With data of the quality presented here, the detection of transits of a "super-Earth" of radius 1.75 R_earth would have been possible.

This paper presents multi-band photometric follow-up observations of the Neptune-mass transiting planet GJ 436b, consisting of 5 new ground-based transit light curves obtained in May 2007. Together with one already published light curve we have at hand a total of 6 light curves, spanning 29 days. The analysis of the data yields an orbital period P = 2.64386+-0.00003 days, mid-transit time T_c [HJD] =2454235.8355+-0.0001, planet mass M_p = 23.1+-0.9 M_{\earth} = 0.073+-0.003 M_{Jup}, planet radius R_p = 4.2+-0.2 R_{\earth} = 0.37+-0.01 R_{Jup} and stellar radius R_s = 0.45+-0.02 R_{\sun}. Our typical precision for the mid transit timing for each transit is about 30 seconds. We searched the data for a possible signature of a second planet in the system through transit timing variations (TTV) and variation of the impact parameter. The analysis could not rule out a small, of the order of a minute, TTV and a long-term modulation of the impact parameter, of the order of +0.2 year^{-1}.

We present the results of a deep (15 ~< r ~< 23), 20 night survey for transiting planets in the intermediate age open cluster M37 (NGC 2099) using the Megacam wide-field mosaic CCD camera on the 6.5m MMT. We do not detect any transiting planets among the ~1450 observed cluster members. We do, however, identify a ~ 1 R_J candidate planet transiting a ~ 0.8 Msun Galactic field star with a period of 0.77 days. The source is faint (V = 19.85 mag) and has an expected velocity semi-amplitude of K ~ 220 m/s (M/M_J). We conduct Monte Carlo transit injection and recovery simulations to calculate the 95% confidence upper limit on the fraction of cluster members and field stars with planets as a function of planetary radius and orbital period. Assuming a uniform logarithmic distribution in orbital period, we find that < 1.1%, < 2.7% and < 8.3% of cluster members have 1.0 R_J planets within Extremely Hot Jupiter (EHJ, 0.4 < T < 1.0 day), Very Hot Jupiter (VHJ, 1.0 < T < 3.0 days) and Hot Jupiter (HJ, 3.0 < T < 5.0 days) period ranges respectively. For 0.5 R_J planets the limits are < 3.2%, and < 21% for EHJ and VHJ period ranges, while for 0.35 R_J planets we can only place an upper limit of < 25% on the EHJ period range. For a sample of 7814 Galactic field stars, consisting primarily of FGKM dwarfs, we place 95% upper limits of < 0.3%, < 0.8% and < 2.7% on the fraction of stars with 1.0 R_J EHJ, VHJ and HJ assuming the candidate planet is not genuine. If the candidate is genuine, the frequency of ~ 1.0 R_J planets in the EHJ period range is 0.002% < f_EHJ < 0.5% with 95% confidence. We place limits of < 1.4%, < 8.8% and < 47% for 0.5 R_J planets, and a limit of < 16% on 0.3 R_J planets in the EHJ period range. This is the first transit survey to place limits on the fraction of stars with planets as small as Neptune.

The nearby transiting system GJ 436b offers a unique opportunity to probe the structure and atmosphere of an extra-solar "hot Neptune". In this Letter, we present the main results of observations covering two transit events with the NICMOS camera on the Hubble Space Telescope. The data consist in high-cadence time series of grism spectra covering the 1.1-1.9 micron spectral range. We find Rpl=4.04 +- 0.10 R_earth and Rstar= 0.446 +- 0.011 Rsun for the planet and star radius, confirming and improving earlier measurements with ground-based photometry and a Spitzer lightcurve at 8 microns, as opposed to a much higher value obtained with the Fine Guidance Sensor on the Hubble Space Telescope. We measure no departure from strict periodicity in the transits to the level of ~7 seconds. This strongly disfavours the proposed explanation of the orbital eccentricity of GJ 436b in terms of the perturbation by another close-by planet. We measure a flat transmission spectrum at the level of a few parts per 10 000 in flux, with no significant signal in the 1.4-micron water band to a level comparable to the maximum amplitude of the effect predicted by planetary atmosphere models.

We present Spitzer 8 micron transit observations of the extrasolar planet HD 149026b. At this wavelength, transit light curves are weakly affected by stellar limb-darkening, allowing for a simpler and more accurate determination of planetary parameters. We measure a planet-star radius ratio of Rp/Rs = 0.05158 +/- 0.00077, and in combination with ground-based data and independent constraints on the stellar mass and radius, we derive an orbital inclination of i = 85.4 +0.9/-0.8 degrees and a planet radius of Rp = 0.755 +/- 0.040 R_jup. These measurements further support models in which the planet is greatly enriched in heavy elements.

Photometry is presented of the Dec. 25, 2007 transit of HD 17156b, which has the longest orbital period and highest orbital eccentricity of all the known transiting exoplanets. New measurements of the stellar radial velocity are also presented. All the data are combined and integrated with stellar-evolutionary modeling to derive refined system parameters. The planet's mass and radius are found to be 3.212_{-0.082}^{+0.069} Jupiter masses and 1.023_{-0.055}^{+0.070} Jupiter radii. The corresponding stellar properties are 1.263_{-0.047}^{+0.035} solar masses and 1.446_{-0.067}^{+0.099} solar radii. The planet is smaller by 1 sigma than a theoretical solar-composition gas giant with the same mass and equilibrium temperature, a possible indication of heavy-element enrichment. The midtransit time is measured to within 1 min, and shows no deviation from a linear ephemeris (and therefore no evidence for orbital perturbations from other planets). We provide ephemerides for future transits and superior conjunctions. There is an 18% chance that the orbital plane is oriented close enough to edge-on for secondary eclipses to occur at superior conjunction. Observations of secondary eclipses would reveal the thermal emission spectrum of a planet that experiences unusually large tidal heating and insolation variations.

We present the results of an archival search for Trans-neptunian objects (TNOs) in an ecliptic field observed with Subaru in 2002. The depth of the search allowed us to find 20 new TNOs with magnitudes between R=24 and 27. We fit a double power law model to the data; the most likely values for the bright and faint power law exponents are alpha_1=0.73_{-0.09}^{+0.08} and alpha_2=0.20_{-0.14}^{+0.12}; the differential number density at R=23 is sigma_{23}=1.46_{-0.12}^{+0.14} and the break magnitude is R_{eq}=25.0_{-0.6}^{+0.8}. This is the most precise measurement of the break in the TNO luminosity function to date. The break in the size distribution corresponds to a diameter of D = 90+-30 km assuming a 4% albedo.

We report new spectroscopic and photometric observations of the parent stars of the recently discovered transiting planets TrES-3 and TrES-4. A detailed abundance analysis based on high-resolution spectra yields [Fe/H] $= -0.19\pm 0.08$, $T_\mathrm{eff} = 5650\pm 75$ K, and $\log g = 4.4\pm 0.1$ for TrES-3, and [Fe/H] $= +0.14\pm 0.09$, $T_\mathrm{eff} = 6200\pm 75$ K, and $\log g = 4.0\pm0.1$ for TrES-4. The accuracy of the effective temperatures is supported by a number of independent consistency checks. The spectroscopic orbital solution for TrES-3 is improved with our new radial-velocity measurements of that system, as are the light-curve parameters for both systems based on newly acquired photometry for TrES-3 and a reanalysis of existing photometry for TrES-4. We have redetermined the stellar parameters taking advantage of the strong constraint provided by the light curves in the form of the normalized separation $a/R_\star$ (related to the stellar density) in conjunction with our new temperatures and metallicities. The masses and radii we derive are $M_\star=0.928_{-0.048}^{+0.028} M_{\sun}$,$R_\star = 0.829_{-0.022}^{+0.015} R_{\sun}$, and $M_\star = 1.404_{-0.134}^{+0.066} M_{\sun}$, $R_\star=1.846_{-0.087}^{+0.096} R_{\sun}$ for TrES-3 and TrES-4, respectively. With these revised stellar parameters we obtain improved values for the planetary masses and radii. We find $M_p = 1.910_{-0.080}^{+0.075} M_\mathrm{Jup}$, $R_p=1.336_{-0.036}^{+0.031} R_\mathrm{Jup}$ for TrES-3, and $M_p=0.925 \pm 0.082 M_\mathrm{Jup}$, $R_p=1.783_{-0.086}^{+0.093} R_\mathrm{Jup}$ for TrES-4. We confirm TrES-4 as the planet with the largest radius among the currently known transiting hot Jupiters.

In the course of conducting a deep (14.5 ~< r ~< 23), 20 night survey for transiting planets in the rich ~550 Myr old open cluster M37 we have measured the rotation periods of 575 stars which lie near the cluster main sequence, with masses 0.2 Msun ~< M ~< 1.3 Msun. This is the largest sample of rotation periods for a cluster older than 500 Myr. Using this rich sample we investigate a number of relations between rotation period, color and the amplitude of photometric variability. Stars with M >~ 0.8 Msun show a tight correlation between period and mass with heavier stars rotating more rapidly. There is a group of 4 stars with P > 15 days that fall well above this relation, which, if real, would present a significant challenge to theories of stellar angular momentum evolution. Below 0.8 Msun the stars continue to follow the period-mass correlation but with a broad tail of rapid rotators that expands to shorter periods with decreasing mass. We combine these results with observations of other open clusters to test the standard theory of lower-main sequence stellar angular momentum evolution. We find that the model reproduces the observations for solar mass stars, but discrepancies are apparent for stars with 0.6 ~< M ~< 1.0 Msun. We also find that for late-K through early-M dwarf stars in this cluster rapid rotators tend to be bluer than slow rotators in B-V but redder than slow rotators in V-I_{C}. This result supports the hypothesis that the significant discrepancy between the observed and predicted temperatures and radii of low-mass main sequence stars is due to stellar activity.

EPOXI (EPOCh + DIXI) is a NASA Discovery Program Mission of Opportunity using the Deep Impact flyby spacecraft. The EPOCh (Extrasolar Planet Observation and Characterization) Science Investigation will gather photometric time series of known transiting exoplanet systems from January through August 2008. Here we describe the steps in the photometric extraction of the time series and present preliminary results of the first four EPOCh targets.

The NASA Discovery mission EPOXI, utilizing the Deep Impact flyby spacecraft, comprises two phases: EPOCh (Extrasolar Planet Observation and Characterization) and DIXI (Deep Impact eXtended Investigation). With EPOCh, we use the 30-cm high resolution visible imager to obtain ultraprecise photometric light curves of known transiting planet systems. We will analyze these data for evidence of additional planets, via transit timing variations or transits; for planetary moons or rings; for detection of secondary eclipses and the constraint of geometric planetary albedos; and for refinement of the system parameters. Over a period of four months, EPOCh observed four known transiting planet systems, with each system observed continuously for several weeks. Here we present an overview of EPOCh, including the spacecraft and science goals, and preliminary photometry results.

We present Spitzer 8 micron transit observations of the extrasolar planet system HD 149026. At this wavelength, transit light curves are weakly affected by stellar limb-darkening, allowing for a simpler and more accurate determination of planetary parameters. We measure a planet-star radius ratio of R_p/R_s = 0.05158 +/- 0.00077, and in combination with ground-based data and independent constraints on the stellar mass and radius, we derive an orbital inclination of i = 85.4 +0.9/-0.8 deg. and a planet radius of 0.755 +/- 0.040 Jupiter radii. These measurements further support models in which the planet is greatly enriched in heavy elements.

With 40 or more transiting exoplanets now known, the time is ripe to seek patterns and correlations among their observed properties, which may give important insights into planet formation, structure, and evolution. This task is made difficult by the widely different methodologies that have been applied to measure their properties in individual cases. Furthermore, in many systems our knowledge of the planet properties is limited by the knowledge of the properties of the parent stars. To address these difficulties we have undertaken the first comprehensive analysis of the data for 23 transiting planets using a uniform methodology. We revisit several of the recently proposed correlations, and find new ones involving the metallicity of the parent stars.

We present photometry of 13 transits of XO-3b, a massive transiting planet on an eccentric orbit. Previous data led to two inconsistent estimates of the planetary radius. Our data strongly favor the smaller radius, with increased precision: R_p = 1.217 +/- 0.073 R_Jup. A conflict remains between the mean stellar density determined from the light curve, and the stellar surface gravity determined from the shapes of spectral lines. We argue the light curve should take precedence, and revise the system parameters accordingly. The planetary radius is about 1 sigma larger than the theoretical radius for a hydrogen-helium planet of the given mass and insolation. To help in planning future observations, we provide refined transit and occultation ephemerides.

We present the results of a survey for trans-neptunian objects (TNOs) based on Subaru archival images, originally collected by Sheppard et al. (2005) as part of a search for irregular satellites of Uranus. The survey region covers 2.8 deg^2, centered on Uranus and observed near opposition on two adjacent nights. Our survey reaches half its maximum detection efficiency at R=25.69$\pm$0.01. The objects detected correspond to 82 TNOs, five Centaurs, and five irregular satellites. We model the cumulative number of TNOs brighter than a given apparent magnitude with both a single and double power law. The best fit single power law, with one object per square degree at magnitude $R_0 =22.6_{-0.4}^{+0.3}$ and a slope of $\alpha = 0.51_{-0.6}^{+0.5}$, is inconsistent with the results of similar searches with shallower limiting magnitudes. The best fit double power law, with a bright-end slope $\alpha_1 = 0.7_{-0.1}^{+0.2}$, a faint-end slope $\alpha_2=0.3_{-0.2}^{+0.2}$, a differential number density at R=23, $\sigma_{23} = 2.0_{-0.5}^{+0.5}$ and a magnitude break in the slope at $R_{eq} = 24.3_{-0.1}^{+0.8}$, is more likely than the single power law by a Bayes factor of ~26. This is the first survey with sufficient depth and areal coverage to identify the magnitude at which the break occurs without relying on the results of other surveys. ...

We present photometry of the exoplanet host star TrES-3 spanning six occultations (secondary eclipses) of its giant planet. No flux decrements were detected, leading to 99%-confidence upper limits on the planet-to-star flux ratio of 0.00024, 0.0005, and 0.00086 in the i, z, and R bands respectively. The corresponding upper limits on the planet's geometric albedo are 0.30, 0.62, and 1.07. The upper limit in the i band rules out the presence of highly reflective clouds, and is only a factor of 2-3 above the predicted level of thermal radiation from the planet.

We have performed an ecliptic imaging survey of the Kuiper belt with our deepest and widest field achieving a limiting flux of m(g') = 26.4, with a sky coverage of 3.0 square-degrees. This is the largest coverage of any other Kuiper belt survey to this depth. We detect 72 objects, two of which have been previously observed. We have improved the Bayesian maximum likelihood fitting technique presented in Gladman et al. (1998) to account for calibration and sky density variations and have used this to determine the luminosity function of the Kuiper belt. Combining our detections with previous surveys, we find the luminosity function is well represented by a single power-law with slope alpha = 0.65 +/- 0.05 and an on ecliptic sky density of 1 object per square-degree brighter than m(R)=23.42 +/- 0.13. Assuming constant albedos, this slope suggests a differential size-distribution slope of 4.25 +/- 0.25, which is steeper than the Dohnanyi slope of 3.5 expected if the belt is in a state of collisional equilibrium. We find no evidence for a roll-over or knee in the luminosity function and reject such models brightward of m(R) ~ 24.6.

We present refined values for the physical parameters of transiting exoplanets, based on a self-consistent and uniform analysis of transit light curves and the observable properties of the host stars. Previously it has been difficult to interpret the ensemble properties of transiting exoplanets, because of the widely different methodologies that have been applied in individual cases. Furthermore, previous studies often ignored an important constraint on the mean stellar density that can be derived directly from the light curve. The main contributions of this work are 1) a critical compilation and error assessment of all reported values for the effective temperature and metallicity of the host stars; 2) the application of a consistent methodology and treatment of errors in modeling the transit light curves; and 3) more accurate estimates of the stellar mass and radius based on stellar evolution models, incorporating the photometric constraint on the stellar density. We use our results to revisit some previously proposed patterns and correlations within the ensemble. We confirm the mass-period correlation, and we find evidence for a new pattern within the scatter about this correlation: planets around metal-poor stars are more massive than those around metal-rich stars at a given orbital period. Likewise, we confirm the proposed dichotomy of planets according to their Safronov number, and we find evidence that the systems with small Safronov numbers are more metal-rich on average. Finally, we confirm the trend that led to the suggestion that higher-metallicity stars harbor planets with a greater heavy-element content.

In the ongoing HATNet survey we have detected a giant planet, with radius 1.33 +/- 0.06 RJup and mass 1.06 +/- 0.12 MJup, transiting the bright (V = 10.5) star GSC 03239-00992. The planet is in a circular orbit with period 3.852985 +/- 0.000005 days and mid-transit epoch 2,454,035.67575 +/- 0.00028 (HJD). The parent star is a late F star with mass 1.29 +/- 0.06 Msun, radius 1.46 +/- 0.06 Rsun, Teff ~ 6570 +/- 80 K, [Fe=H] = -0.13 +/- 0.08 and age ~ 2.3+/-^{0.5}_{0.7}Gy. With this radius and mass, HAT-P-6b has somewhat larger radius than theoretically expected. We describe the observations and their analysis to determine physical properties of the HAT-P-6 system, and briefly discuss some implications of this finding.

We have conducted a deep ($15 \la r \la 23$), 20 night survey for transiting planets in the intermediate age open cluster M37 (NGC 2099) using the Megacam wide-field mosaic CCD camera on the 6.5m Multiple Mirror Telescope (MMT). In this paper we describe the observations and data reduction procedures for the survey and analyze the stellar content and dynamical state of the cluster. By combining high resolution spectroscopy with existing $BVI_{C}K_{S}$ and new $gri$ color magnitude diagrams we determine the fundamental cluster parameters: $t = 485 \pm 28$ Myr without overshooting ($t = 550 \pm 30 {\rm Myr}$ with overshooting), $E(B-V) = 0.227 \pm 0.038$, $(m-M)_{V} = 11.57 \pm 0.13$ and $[M/H] = +0.045 \pm 0.044$ which are in good agreement with, though more precise than, previous measurements. We determine the mass function down to $0.3 M_{\odot}$ and use this to estimate the total cluster mass of $3640 \pm 170 M_{\odot}$.

We have conducted a deep ($15 \la r \la 23$), 20 night survey for transiting planets in the intermediate age ($\sim 550 {\rm Myr}$) open cluster M37 (NGC 2099) using the Megacam wide-field mosaic CCD camera on the 6.5m Multiple Mirror Telescope (MMT). In this paper we present a catalog and light curves for 1445 variable stars; 1430 (99%) of these are new discoveries. We have discovered 20 new eclipsing binaries and 31 new short period ($P < 1 {\rm day}$) pulsating stars. The bulk of the variables are most likely rapidly rotating young low-mass stars, including a substantial number ($\ga 500$) that are members of the cluster. We identify and analyze five particularly interesting individual variables including a previously identified variable which we suggest is probably a hybrid $\gamma$ Doradus/$\delta$ Scuti pulsator, two possible quiescent cataclysmic variables, a detached eclipsing binary (DEB) with at least one $\gamma$ Doradus pulsating component (only the second such variable found in an eclipsing binary), and a low mass ($M_{P} \sim M_{S} \sim 0.6 M_{\odot}$) DEB that is a possible cluster member. A preliminary determination of the physical parameters for the DEB+$\gamma$ Doradus system yields $M_{P} = 1.58 \pm 0.04 M_{\odot}$, $M_{S} = 1.58 \pm 0.04 M_{\odot}$, $R_{P} = 1.39 \pm 0.07 R_{\odot}$ and $R_{S} = 1.38 \pm 0.07 R_{\odot}$.

We summarize the results of two experiments to address important issues related to the correlation between planet frequencies and properties and the metallicity of the hosts. Our results can usefully inform formation, structural, and evolutionary models of gas giant planets.

We present new photometry of HD 149026 spanning five transits of its "super-Neptune" planet. In combination with previous data, we improve upon the determination of the planet-to-star radius ratio: R_p/R_star = 0.0491^{+0.0018}_{-0.0005}. We find the planetary radius to be 0.71 +/- 0.05 R_Jup, in accordance with previous theoretical models invoking a high metal abundance for the planet. The limiting error is the uncertainty in the stellar radius. Although we find agreement among four different ways of estimating the stellar radius, the uncertainty remains at 7%. We also present a refined transit ephemeris and a constraint on the orbital eccentricity and argument of pericenter, e cos(omega) = -0.0014 +/- 0.0012, based on the measured interval between primary and secondary transits.

We present photometry of the G0 star HAT-P-1 during six transits of its close-in giant planet, and we refine the estimates of the system parameters. Relative to Jupiter's properties, HAT-P-1b is 1.20 +/- 0.05 times larger and its surface gravity is 2.7 +/- 0.2 times weaker. Although it remains the case that HAT-P-1b is among the least dense of the known sample of transiting exoplanets, its properties are in accord with previously published models of strongly irradiated, coreless, solar-composition giant planets. The times of the transits have a typical accuracy of 1 min and do not depart significantly from a constant period.

(Abridged) The discovery of extrasolar planetary systems revealed an unexpected diversity of planetary systems that has revolutionized planet formation theory. A strong program of theoretical research is essential to maximize both the discovery potential and the scientific returns of future observational programs, so as to achieve a deeper understanding of the formation and evolution of planetary systems. We outline three broad categories of theoretical research: detailed studies of specific planetary systems, testing planet formation models by comparing their predictions to the observed exoplanet population, and detailed modeling of specific physical processes. We describe how such theoretical research plays an important role in analyzing observations for a wide range detection methods and contributes to understanding the Earth's place in the universe and the potential for Earth-like life beyond our solar system. In this white paper, we suggest how to maximize the scientific return of future exoplanet observations. Our recommendations include a strong theory program, support for multiple observational programs that will study a diverse set of planets and stars, significant observing time devoted to follow-up observations, and healthy collaboration between observers and theorists.

Theoretical studies predict that Trojans are likely a frequent byproduct of planet formation and evolution. We examine the sensitivity of transit timing observations for detecting Trojan companions to transiting extrasolar planets. We demonstrate that this method offers the potential to detect terrestrial-mass Trojans using existing ground-based observatories. We compare the transit timing variation (TTV) method with other techniques for detecting extrasolar Trojans and outline the future prospects for this method.

Of the nearby transiting exoplanets that are amenable to detailed study, TrES-2 is both the most massive and has the largest impact parameter. We present z-band photometry of three transits of TrES-2. We improve upon the estimates of the planetary, stellar, and orbital parameters, in conjunction with the spectroscopic analysis of the host star by Sozzetti and co-workers. We find the planetary radius to be 1.222 +/- 0.038 R_Jup and the stellar radius to be 1.003 +/- 0.027 R_Sun. The quoted uncertainties include the systematic error due to the uncertainty in the stellar mass (0.980 +/- 0.062 M_Sun). The timings of the transits have an accuracy of 25s and are consistent with a uniform period, thus providing a baseline for future observations with the NASA Kepler satellite, whose field of view will include TrES-2.

We report on a spectroscopic determination of the atmospheric parameters and chemical abundance of the parent star of the recently discovered transiting planet {TrES-2}. A detailed LTE analysis of a set of \ion{Fe}{1} and \ion{Fe}{2} lines from our Keck spectra yields $T_\mathrm{eff} = 5850\pm 50$ K, $\log g = 4.4\pm 0.1$, and [Fe/H] $= -0.15\pm 0.10$. Several independent checks (e.g., additional spectroscopy, line-depth ratios) confirm the reliability of our spectroscopic $T_\mathrm{eff}$ estimate. The mass and radius of the star, needed to determine the properties of the planet, are traditionally inferred by comparison with stellar evolution models using $T_\mathrm{eff}$ and some measure of the stellar luminosity, such as the spectroscopic surface gravity (when a trigonometric parallax is unavailable, as in this case). We apply here a new method in which we use instead of $\log g$ the normalized separation $a/R_\star$ (related to the stellar density), which can be determined directly from the light curves of transiting planets with much greater precision. With the $a/R_\star$ value from the light curve analysis of Holman et al. \citeyearpar{holman07b} and our $T_\mathrm{eff}$ estimate we obtain $M_\star = 0.980\pm0.062 M_\odot$ and $R_\star = 1.000_{-0.033}^{+0.036} R_\odot$, and an evolutionary age of $5.1^{+2.7}_{-2.3}$ Gyr, in good agreement with other constraints based on the strength of the emission in the \ion{Ca}{2} H & K line cores, the Lithium abundance, and rotation. The new stellar parameters yield improved values for the planetary mass and radius of $M_p = 1.198 \pm 0.053 M_\mathrm{Jup}$ and $R_p = 1.220^{+0.045}_{-0.042} R_\mathrm{Jup}$, confirming that {TrES-2} is the most massive among the currently known nearby ($d\lesssim 300$ pc) transiting hot Jupiters. [Abridged]

In the coming decades, research in extrasolar planets aims to advance two goals: 1) detecting and characterizing low-mass planets increasingly similar to the Earth, and 2) improving our understanding of planet formation. We present a new planet detection method that is capable of making large advances towards both of these objectives and describe a modest network of telescopes that is able to make the requisite observations. In a system where a known planet transits its host star, a second planet in that system will cause the time between transits to vary. These transit timing variations can be used to infer the orbital elements and mass of the perturbing planet even if it has a mass that is smaller than the mass of the Earth. This detection technique complements other techniques because it is most sensitive in mean-motion resonances where, due to degeneracies, other techniques have reduced sensitivity. Small ground-based observatories have already exceeded the photometric precision necessary to detect sub-Earth mass planets. However, TTV planet searches are currently limited by the relatively small number of high-precision transit data and insufficient observing time on existing telescopes. These issues will be compounded as the number of known transiting planets suitable for TTV study will increase substantially in the near future. A relatively modest investment in a ground-based network of small ($\sim 0.5 {\rm m}$ telescopes could provide the needed coverage and so dramatically increase the effectiveness of transit timing observations.

We are conducting a transit survey of the open cluster M37 using the Megacam instrument on the 6.5 m Multiple-Mirror Telescope. We have obtained ~4500 images of this cluster over 18.5 nights and have achieved the precision necessary to detect planets smaller than Saturn. In this presentation we provide an overview of the project, describe the ongoing data reduction/analysis and present some of our preliminary results.

We present photometry of HD 189733 during eight transits of its close-in giant planet, and out-of-transit photometry spanning two years. Using the transit photometry, we determine the stellar and planetary radii and the photometric ephemeris. Outside of transits, there are quasiperiodic flux variations with a 13.4 day period that we attribute to stellar rotation. In combination with previous results, we derive upper limits on the orbital eccentricity, and on the true angle between the stellar rotation axis and planetary orbit (as opposed to the angle between the projections of those axes on the sky).

We present z band photometry of three consecutive transits of the exoplanet TrES-1, with an accuracy of 0.15% and a cadence of 40 seconds. We improve upon estimates of the system parameters, finding in particular that the planetary radius is 1.081 +/- 0.029 Jupiter radii, and the stellar radius is 0.811 +/- 0.020 solar radii. The uncertainties include both the statistical error and the systematic error arising from the uncertainty in the stellar mass. The transit times are determined to within about 15 seconds, and allow us to refine the estimate of the mean orbital period: P = 3.0300737 +/- 0.0000026 days. We find no evidence for star spots or other irregularities that have been previously reported.

We present spectroscopy of a transit of the exoplanet HD 189733b. By modeling the Rossiter-McLaughlin effect (the anomalous Doppler shift due to the partial eclipse of the rotating stellar surface), we find the angle between the sky projections of the stellar spin axis and orbit normal to be lambda = -1.4 +/- 1.1 deg. This is the third case of a ``hot Jupiter'' for which lambda has been measured. In all three cases lambda is small, ruling out random orientations with 99.96% confidence, and suggesting that the inward migration of hot Jupiters generally preserves spin-orbit alignment.

We present I and B photometry of five distinct transits of the exoplanet OGLE-TR-10b. By modeling the light curves, we find the planetary radius to be R_P = 1.06 +/- 0.08 R_Jup and the stellar radius to be R_S = 1.10 +/- 0.07 R_sun. The uncertainties are dominated by statistical errors in the photometry. Our estimate of the planetary radius is smaller than previous estimates that were based on lower-precision photometry, and hence the planet is not as anomalously large as was previously thought. We provide updated determinations of all the system parameters, including the transit ephemerides.

We present precise z-band photometric time series spanning times of transit of the two exoplanets recently discovered by the SuperWASP collaboration. We find planetary radii of 1.44 +/- 0.08 R_J and 1.04 +/- 0.06 R_J for WASP-1b and WASP-2b, respectively. These error estimates include both random errors in the photometry and also the uncertainty in the stellar masses. Our results are 5 times more precise than the values derived from the discovery data alone. Our measurement of the radius of WASP-2b agrees with previously published models of hot Jupiters that include both a 20-M_Earth core of solid material and the effects of stellar insolation. In contrast, we find that the models cannot account for the large size of WASP-1b, even if the planet has no core. Thus, we add WASP-1b to the growing list of hot Jupiters that are larger than expected. This suggests that ``inflated'' hot Jupiters are more common than previously thought, and that any purported explanations involving highly unusual circumstances are disfavored.

As part of our ongoing effort to measure exoplanet sizes and transit times with greater accuracy, we present I band observations of two transits of OGLE-TR-111b. The photometry has an accuracy of 0.15-0.20% and a cadence of 1-2 minutes. We derive a planetary radius of 1.067 +/- 0.054 Jupiter radii and a stellar radius of 0.831 +/- 0.031 solar radii. The uncertainties are dominated by errors in the photometry, rather than by systematic errors arising from uncertainties in the limb darkening function or the stellar mass. Both the stellar radius and the planetary radius are in agreement with theoretical expectations. The transit times are accurate to within 30 seconds, and allow us to refine the estimate of the mean orbital period: 4.0144479 +/- 0.0000041 days.

The Edgeworth-Kuiper belt encodes the dynamical history of the outer solar system. Kuiper belt objects (KBOs) bear witness to coagulation physics, the evolution of planetary orbits, and external perturbations from the solar neighborhood. We critically review the present-day belt's observed properties and the theories designed to explain them. Theories are organized according to a possible time-line of events. In chronological order, epochs described include (1) coagulation of KBOs in a dynamically cold disk, (2) formation of binary KBOs by fragmentary collisions and gravitational captures, (3) stirring of KBOs by Neptune-mass planets (``oligarchs''), (4) eviction of excess oligarchs, (5) continued stirring of KBOs by remaining planets whose orbits circularize by dynamical friction, (6) planetary migration and capture of Resonant KBOs, (7) creation of the inner Oort cloud by passing stars in an open stellar cluster, and (8) collisional comminution of the smallest KBOs. Recent work underscores how small, collisional, primordial planetesimals having low velocity dispersion permit the rapid assembly of ~5 Neptune-mass oligarchs at distances of 15-25 AU. We explore the consequences of such a picture. We propose that Neptune-mass planets whose orbits cross into the Kuiper belt for up to ~20 Myr help generate the high-perihelion members of the hot Classical disk and Scattered belt. By contrast, raising perihelia by sweeping secular resonances during Neptune's migration might fill these reservoirs too inefficiently when account is made of how little primordial mass might reside in bodies having sizes of order 100 km. These and other frontier issues in trans-Neptunian space are discussed quantitatively.

We present RIz photometry of four consecutive transits of the newly discovered exoplanet XO-1b. We improve upon the estimates of the transit parameters, finding the planetary radius to be R_P = 1.184 +0.028/-0.018 R_Jupiter and the stellar radius to be R_S = 0.928 +0.018/-0.013 R_Sun, assuming a stellar mass of M_S = 1.00 +/- 0.03 M_Sun. The uncertainties in the planetary and stellar radii are dominated by the uncertainty in the stellar mass. These uncertainties increase by a factor of 2-3 if a more conservative uncertainty of 0.10 M_Sun is assumed for the stellar mass. Our estimate of the planetary radius is smaller than that reported by McCullough et al. (2006) and yields a mean density that is comparable to that of TrES-1 and HD 189733b. The timings of the transits have an accuracy ranging from 0.2 to 2.5 minutes, and are marginally consistent with a uniform period.

The unusual flux variations of the pre-main-sequence binary star KH 15D have been attributed to occultations by a circumbinary disk. We test whether or not this theory is compatible with newly available data, including recent radial velocity measurements, CCD photometry over the past decade, and photographic photometry over the past 50 years. We find the model to be successful, after two refinements: a more realistic motion of the occulting feature, and a halo around each star that probably represents scattering by the disk. The occulting feature is exceptionally sharp-edged, raising the possibility that the dust in the disk has settled into a thin layer, and providing a tool for fine-scale mapping of the immediate environment of a T Tauri star. However, the window of opportunity is closing, as the currently visible star may be hidden at all orbital phases by as early as 2008.

We present the data acquisition strategy and characterization procedures for the Canada-France Ecliptic Plane Survey (CFEPS), a sub-component of the Canada-France-Hawaii Telescope Legacy Survey. The survey began in early 2003 and as of summer 2005 has covered 430 square degrees of sky within a few degrees of the ecliptic. Moving objects beyond the orbit of Uranus are detected to a magnitude limit of $m_R$=23 -- 24 (depending on the image quality). To track as large a sample as possible and avoid introducing followup bias, we have developed a multi-epoch observing strategy that is spread over several years. We present the evolution of the uncertainties in ephemeris position and orbital elements as the objects progress through the epochs. We then present a small 10-object sample that was tracked in this manner as part of a preliminary survey starting a year before the main CFEPS project. We describe the CFEPS survey simulator, to be released in 2006, which allows theoretical models of the Kuiper Belt to be compared with the survey discoveries since CFEPS has a well-documented pointing history with characterized detection efficiencies as a function of magnitude and rate of motion on the sky. Using the pre-survey objects we illustrate the usage of the simulator in modeling the classical Kuiper Belt.

We determine the stellar, planetary, and orbital properties of the transiting planetary system HD 209458, through a joint analysis of high-precision radial velocities, photometry, and timing of the secondary eclipse. Of primary interest is the strong detection of the Rossiter-McLaughlin effect, the alteration of photospheric line profiles that occurs because the planet occults part of the rotating surface of the star. We develop a new technique for modeling this effect, and use it to determine the inclination of the planetary orbit relative to the apparent stellar equator [lambda = (-4.4 +/- 1.4) degrees], and the line-of-sight rotation speed of the star [v*sin(I) = (4.70 +/- 0.16) km/s]. The uncertainty in these quantities has been reduced by an order of magnitude relative to the pioneering measurements by Queloz and collaborators. The small but nonzero misalignment is probably a relic of the planet formation epoch, because the expected timescale for tidal coplanarization is larger than the age of the star. Our determination of v*sin(I) is a rare case in which rotational line broadening has been isolated from other broadening mechanisms.

We present precise, ~1%, r-band relative photometry of the unusual solar system object (90377) Sedna. Our data consist of 143 data points taken over eight nights in October 2004 and January 2005. The RMS variability over the longest contiguous stretch of five nights of data spanning nine days is only 1.3%. This subset of data alone constrain the amplitude of any long-period variations with period P to be A<1% (P/20 days)^2. Over the course of any given 5-hour segment, the data exhibits significant linear trends not seen in a comparison star of similar magnitude, and in a few cases these segments show clear evidence for curvature at the level of a few millimagnitudes per hour^2. These properties imply that the rotation period of Sedna is O(10 hours), cannot be < 5 hours, and cannot be >10 days, unless the intrinsic light curve has significant and comparable power on multiple timescales, which is unlikely. A sinusoidal fit yields a period of P=(10.273 +/- 0.002) hours and semi-amplitude of A=(1.1 +/- 0.1)%. There are additional acceptable fits with flanking periods separated by ~3 minutes, as well as another class of fits with P ~ 18 hours, although these later fits appear less viable based on visual inspection. Our results indicate that the period of Sedna is likely consistent with typical rotation periods of solar system objects, thus obviating the need for a massive companion to slow its rotation.

We report g, V, and r photometric time series of HD 149026 spanning predicted times of transit of the Saturn-mass planetary companion, which was recently discovered by Sato and collaborators. We present a joint analysis of our observations and the previously reported photometry and radial velocities of the central star. We refine the estimate of the transit ephemeris to Tc [HJD] = 2453527.87455^{+0.00085}_{-0.00091} + N * 2.87598^{+0.00012}_{-0.00017}. Assuming that the star has a radius of 1.45 +/- 0.10 R_Sun and a mass of 1.30 +/- 0.10 M_Sun, we estimate the planet radius to be 0.726 +/- 0.064 R_Jup, which implies a mean density of 1.07^{+0.42}_{-0.30} g/cm^3. This density is significantly greater than that predicted for models which include the effects of stellar insolation and for which the planet has only a small core of solid material. Thus we confirm that this planet likely contains a large core, and that the ratio of core mass to total planet mass is more akin to that of Uranus and Neptune than that of either Jupiter or Saturn.

We present high-precision photometry of the hypervelocity star SDSS J090745.0+024507 (HVS), which has a Galactic rest-frame radial velocity of v=709 km/s, and so has likely been ejected from the supermassive black hole in the Galactic center. Our data were obtained on two nights using the MMT 6.5m telescope, and is supplemented by lower precision photometry obtained on four nights using the FLWO 1.2m telescope. The high-precision photometry indicates that the HVS is a short-period, low-amplitude variable, with period P=0.2-2 days and amplitude A = 2-10%. Together with the known effective temperature of T_eff ~ 10,500 K (spectral type B9), this variability implies that the HVS is a member of the class of slowly pulsating B-type main sequence stars, thus resolving the previously-reported two-fold degeneracy in the luminosity and distance of the star. The HVS has a heliocentric distance of 71 kpc, and an age of ~0.35 Gyr. The time of ejection from the center of the Galaxy is < 100 Myr, and thus the existence of the OS constitutes observational evidence of a population of young stars in the proximity of the central supermassive black hole ~0.1 Gyr ago. It is possible that the HVS was a member of a binary that was tidally disrupted by the central black hole; we discuss constraints on the properties of the companion's orbit.

Obliquity tides are a potentially important source of heat for extrasolar planets on close-in orbits. Although tidal dissipation will usually reduce the obliquity to zero, a nonzero obliquity can persist if the planet is in a Cassini state, a resonance between spin precession and orbital precession. Obliquity tides might be the cause of the anomalously large size of the transiting planet HD 209458b.

We present the results from a three night, time-series study of the open cluster NGC 6791 using the Megacam wide-field mosaic CCD camera on the 6.5m MMT telescope. The aim of this study was to demonstrate the ability to obtain very high precision photometry for a large number of stars. We achieved better than 1% precision for more than 8000 stars with 14.3 < R < 20.1 and sub-millimagnitude (as low as 0.36 mmag) precision for over 300 stars with 14.6 < R < 16.3 in the field of this cluster. We also discovered 8 new variable stars, including a possible delta-Scuti variable with an amplitude of 2%, 6 likely W UMa contact binaries, and a possible RS CVn star, and we identified 5 suspected low-amplitude variables, including one star with an amplitude of 3 mmag. We comment on the implications of this study for a ground-based survey for transiting planets as small as Neptune.

We present deep optical photometry of the afterglow of gamma-ray burst (GRB) 041006 and its associated hypernova obtained over 65 days after detection (55 R-band epochs on 10 different nights). Our early data (t<4 days) joined with published GCN data indicates a steepening decay, approaching F_nu ~t^{-0.6} at early times (<<1 day) and F_nu ~t^{-1.3} at late times. The break at t_b=0.16+-0.04 days is the earliest reported jet break among all GRB afterglows. During our first night, we obtained 39 exposures spanning 2.15 hours from 0.62 to 0.71 days after the burst that reveal a smooth afterglow, with an rms deviation of 0.024 mag from the local power-law fit, consistent with photometric errors. After t~4 days, the decay slows considerably, and the light curve remains approximately flat at R~24 mag for a month before decaying by another magnitude to reach R~25 mag two months after the burst. This ``bump'' is well-fitted by a k-corrected light curve of SN1998bw, but only if stretched by a factor of 1.38 in time. In comparison with the other GRB-related SNe bumps, GRB 041006 stakes out new parameter space for GRB/SNe, with a very bright and significantly stretched late-time SN light curve. Within a small sample of fairly well observed GRB/SN bumps, we see a hint of a possible correlation between their peak luminosity and their ``stretch factor'', broadly similar to the well-studied Phillips relation for the type Ia supernovae.

Future surveys for transiting extrasolar planets, including the space-based mission Kepler (Borucki et al 2003), are expected to detect hundreds of Jovian mass planets and tens of terrestrial mass planets. For many of these newly discovered planets, the intervals between successive transits will be measured with an accuracy of 0.1--100 minutes. We show that these timing measurements will allow for the detection of additional planets in the system (not necessarily transiting), via their gravitational interaction with the transiting planet. The transit time variations depend on the mass of the additional planet, and in some cases Earth-mass planets will produce a measurable effect.

[Condensed] We search 0.02 deg^2 for trans-Neptunian objects (TNOs) with m<=29.2 (diameter ~15 km) using the ACS on HST. Three new objects are discovered, roughly 25 times fewer than expected from extrapolation of the differential sky density Sigma(m) of brighter objects. The ACS and other recent TNO surveys show departures from a power law size distribution. Division of the TNO sample into ``classical Kuiper belt'' (CKB) and ``Excited'' samples reveals that Sigma(m) differs for the two populations at 96% confidence. A double power law adequately fits all data. Implications include: The total mass of the CKB is ~0.010 M_Earth, only a few times Pluto's mass, and is predominately in the form of ~100 km bodies. The mass of Excited objects is perhaps a few times larger. The Excited class has a shallower bright-end size distribution; the largest objects, including Pluto, comprise tens of percent of the total mass whereas the largest CKBOs are only ~2% of its mass. The predicted mass of the largest Excited body is close to the Pluto mass; the largest CKBO is ~60 times less massive. The deficit of small TNOs occurs for sizes subject to disruption by present-day collisions, suggesting extensive depletion by collisions. Both accretion and erosion appearing to have proceeded to more advanced stages in the Excited class than the CKB. The absence of distant TNOs implies that any distant (60 AU) population must have less than the CKB mass in the form of objects 40 km or larger. The CKB population is sparser than theoretical estimates of the required precursor population for short period comets, but the Excited population could be a viable precursor population.

We present BVR photometric colors of six Uranian and two Neptunian irregular satellites, collected using the Magellan Observatory (Las Campanas, Chile) and the Keck Observatory, (Manua Kea, Hawaii). The colors range from neutral to light red, and like the Jovian and the Saturnian irregulars (Grav et al. 2003) there is an apparent lack of the extremely red objects found among the Centaurs and Kuiper belt objects. The Uranian irregulars can be divided into three possible dynamical families, but the colors collected show that two of these dynamical families, the Caliban and Sycorax-clusters, have heterogeneous colors. Of the third possible family, the 168-degree cluster containing two objects with similar average inclinations but quite different average semi-major axis, only one object (U XXI Trinculo) was observed. The heterogeneous colors and the large dispersion of the average orbital elements leads us to doubt that they are collisional families. We favor single captures as a more likely scenario. The two neptunians observed (N II Nereid and S/2002 N1) both have very similar neutral, sun-like colors. Together with the high collisional probability between these two objects over the age of the solar system (Nesvorny et al. 2003, Holman et al. 2004), this suggests that S/2002 N1 be a fragment of Nereid, broken loose during a collision or cratering event with an undetermined impactor.

Precise-Doppler experiments suggest that a massive (m sin i=0.86 M_J) planet orbits at semimajor axis a=3.4 AU around Epsilon Eridani, a nearby star with a massive debris disk. The dynamical perturbations from such a planet would mold the distribution of dust around this star. We numerically integrated the orbits of dust grains in this system to predict the central dust cloud structure. For a supply of grains that begin in low-inclination, low-eccentricity orbits at 15 AU, the primary feature of the dust distribution is a pair of dense clumps, containing dust particles trapped in mean-motion resonances of the form n:1. These clumps appear to revolve around the star once every two planet revolutions. Future observations with the IRAM Plateau de Bure Interferometer, the SMA, or ALMA could detect these clumps, confirming the existence of the planet and revealing its location.

We present JHKs photometry of 10 Jovian and 4 Saturnian irregular satellites, taken with the Near-InfraRed Imager (NIRI) at the 8-m Gemini North Observatory on Mauna Kea, Hawaii. The observed objects have near-infrared colors consistent with C, P and D-type asteroids, although J XII Ananke and S IX Phoebe show weak indications of possible water features in the H filter. The four members of the Himalia-family have similar near-infrared colors, as do the two members of the Gallic family, S XX Paaliaq and S XXIX Siarnaq. From low resolution normalized reflectance spectra based on the broadband colors and covering 0.4 to 2.2 microns, the irregular satellites are identified as C-type (J VII Pasiphae, J VI Himalia and S IX Phoebe), P-type (J XII Ananke and J XVIII Themisto) and D-type (J IX Carme and J X Sinope), showing a diversity of origins of these objects.

We propose that the extraordinary ``winking star'' KH 15D is an eccentric pre-main-sequence binary that is gradually being occulted by an opaque screen. This model accounts for the periodicity, depth, duration, and rate of growth of the modern eclipses; the historical light curve from photographic plates; and the existing radial velocity measurements. It also explains the re-brightening events that were previously observed during eclipses, and the subsequent disappearance of these events. We predict the future evolution of the system and its full radial velocity curve. Given the small velocity of the occulting screen relative to the center of mass of the binary, the screen is probably associated with the binary, and may be the edge of a precessing circumbinary disk.

We determine the period, p = 11.52 \pm 0.14 h, and a light curve peak-to-peak amplitude, a = 0.029 \pm 0.003 magnitudes, of the Neptunian irregular satellite Nereid. If the light curve variation is due to albedo variations across the surface, rather than solely to the shape of Nereid variations, the rotation period would be a factor of two shorter. In either case, such a rotation period and light curve amplitude, together with Nereid's orbital period, p=360.14 days, imply that Nereid is almost certainly in a regular rotation state, rather than the chaotic rotation state suggested by Schaefer and Schaefer (1988,2000) and Dobrovolskis (1995). Assuming that Nereid is perfectly spherical, the albedo variation is 3% across the observed surface. Assuming a uniform geometric albedo, the observed cross sectional area varies by 3%. We caution that the lightcurve found in this paper only sets limits on the combination of albedo and physical irregularity and that we cannot determine the orientation of Nereid's spin axis from our data.

Using simple geometrical arguments, we paint an overview of the variety of resonant structures a single planet with moderate eccentricity (e < ~0.6) can create in a dynamically cold, optically thin dust disk. This overview may serve as a key for interpreting images of perturbed debris disks and inferring the dynamical properties of the planets responsible for the perturbations. We compare the resonant geometries found in the solar system dust cloud with observations of dust clouds around Vega and other stars and we discuss a new model for the asymmetries in the Epsilon Eridani cloud.

We present BVRI colors of 13 Jovian and 8 Saturnian irregular satellites obtained with the 2.56m Nordic Optical Telescope on La Palma, the 6.5m Magellan Baade Telescope on La Campanas, and the 6m MMT on Mt. Hopkins. The observations were performed between December 2001 to March 2002. Nearly all of the known irregular satellites can be divided into two distinct classes based on their colors. One, the grey color class, has the similar colors to the C-type asteroid, and the other, the light red color class, has colors similar to P/D-type asteroids. We also find at least one object, the Jovian irregular J XXIII Kalyke, that has colors similar to the red colored Centaurs/TNOs, although its classification is unsecure. We also find that there is a correlation between the physical properties and dynamical properties of the irregular satellites. Most of the dynamical clusters have homogeneous colors, which points to single homogeneous progenitors being cratered or fragmented as the source of each individual cluster. The heterogeneous colored clusters are most easily explained by assuming that there are several dynamical clusters in the area, rather than just one.

We report UBVRI observations of the optical afterglow of the gamma-ray burst GRB 021004. We observed significant (10-20%) deviations from a power law decay on several time scales, ranging from a few hours down to 20-30 minutes. We also observed a significant color change starting ~1.5 days after the burst, confirming the spectroscopic results already reported by Matheson et al. (2002). We discuss these results in the context of several models that have recently been proposed to account for the anomalous photometric behavior of this event.

ESPI has been proposed for direct imaging and spectral analysis of giant planets orbiting solar-type stars. ESPI extends the concept suggested by Nisenson and Papaliolios (2001) for a square aperture apodized telescope that has sufficient dynamic range to directly detect exo-planets. With a 1.5 M square mirror, ESPI can deliver high dynamic range imagery as close as 0.3 arcseconds to bright sources, permitting a sensitive search for exoplanets around nearby stars and a study of their characteristics in reflected light.

We report polarization measurements and photometry for the optical afterglow of the gamma-ray burst GRB 020405. We measured a highly significant 9.9% polarization (in V band) 1.3 days after the burst and argue that it is intrinsic to the GRB. The light curve decay is well fitted by a $t^{-1.72}$ power-law; we do not see any evidence for a break between 1.24 and 4.3 days after the burst. We discuss these measurements in the light of several models of GRB afterglows.

We numerically integrated the orbits of 1458 particles in the region of the classical Kuiper Belt (41 AU < a < 47 AU) to explore the role of dynamical instabilities in sculpting the inclination distribution of the classical Kuiper Belt Objects (KBOs). We find that the selective removal of low-inclination objects by overlapping secular resonances (nu_17 and nu_18) acts to raise the mean inclination of the surviving population of particles over 4 billion years of interactions with Jupiter, Saturn, Uranus and Neptune, though these long-term dynamical effects do not themselves appear to explain the discovery of KBOs with inclinations near 30 degrees. Our integrations also imply that after 3 billion years of interaction with the massive planets, high inclination KBOs more efficiently supply Neptune-encountering objects, the likely progenitors of short-period comets, Centaurs, and scattered KBOs. The secular resonances at low inclinations may indirectly cause this effect by weeding out objects unprotected by mean motion resonances during the first 3 billion years.

We present images of the Vega system obtained with the IRAM Plateau de Bure interferometer at 1.3 millimeters wavelength with sub-mJy sensitivity and $\sim2\farcs5$ resolution (about 20 AU). These observations clearly detect the stellar photosphere and two dust emission peaks offset from the star by $9\farcs5$ and $8\farcs0$ to the northeast and southwest, respectively. These offset emission peaks are consistent with the barely resolved structure visible in previous submillimeter images, and they account for a large fraction of the dust emission. The presence of two dust concentrations at the observed locations is plausibly explained by the dynamical influence of an unseen planet of a few Jupiter masses in a highly eccentric orbit that traps dust in principle mean motion resonances.

Our understanding of the Solar System has been revolutionized over the past decade by the finding that the orbits of the planets are inherently chaotic. In extreme cases, chaotic motions can change the relative positions of the planets around stars, and even eject a planet from a system. Moreover, the spin axis of a planet-Earth's spin axis regulates our seasons-may evolve chaotically, with adverse effects on the climates of otherwise biologically interesting planets. Some of the recently discovered extrasolar planetary systems contain multiple planets, and it is likely that some of these are chaotic as well.

The physical basis of chaos in the solar system is now better understood: in all cases investigated so far, chaotic orbits result from overlapping resonances. Perhaps the clearest examples are found in the asteroid belt. Overlapping resonances account for its Kirkwood gaps and were used to predict and find evidence for very narrow gaps in the outer belt. Further afield, about one new ``short-period'' comet is discovered each year. They are believed to come from the ``Kuiper Belt'' (at 40 AU or more) via chaotic orbits produced by mean-motion and secular resonances with Neptune. Finally, the planetary system itself is not immune from chaos. In the inner solar system, overlapping secular resonances have been identified as the possible source of chaos. For example, Mercury, in 10^{12} years, may suffer a close encounter with Venus or plunge into the Sun. In the outer solar system, three-body resonances have been identified as a source of chaos, but on an even longer time scale of 10^9 times the age of the solar system. On the human time scale, the planets do follow their orbits in a stately procession, and we can predict their trajectories for hundreds of thousands of years. That is because the mavericks, with shorter instability times, have long since been ejected. The solar system is not stable; it is just old!

We examine the eccentricity evolution of a system of two planets locked in a mean motion resonance, in which the outer planet loses energy and angular momentum. The sink of energy and angular momentum could be either a gas or planetesimal disk. We show that the eccentricity of both planetary bodies can grow to large values, particularly if the inner body does not directly exchange energy or angular momentum with the disk. We analytically calculate the eccentricity damping rate in the case of a single planet migrating through a planetesimal disk. We present the results of numerical integrations of two resonant planets showing rapid growth of eccentricity. We also present integrations in which a Jupiter-mass planet is forced to migrate inward through a system of 5-10 roughly Earth mass planets. The migrating planet can eject or accrete the smaller bodies; roughly 5% of the mass (averaged over all the integrations) accretes onto the central star. The results are discussed in the context of the currently known extrasolar planetary systems.

By telescopic tracking, we have established that the orbit of the trans-neptunian object (2000 CR$_{105}$) has a perihelion of $\simeq$44 AU, and is thus outside the domain controlled by strong gravitational close encounters with Neptune. Because this object is on a very large, eccentric orbit (with semimajor axis $a\simeq$216 AU and eccentricity $e\simeq$0.8) this object must have been placed on this orbit by a gravitational perturbation which is {\it not} direct gravitational scattering off of any of the giant planets (on their current orbits). The existence of this object may thus have profound cosmogonic implications for our understanding of the formation of the outer Solar System. We discuss some viable scenarios which could have produced it, including long-term diffusive chaos and scattering off of other massive bodies in the outer Solar System. This discovery implies that there must be a large population of trans-neptunian objects in an `extended scattered disk' with perihelia above the previously-discussed 38 AU boundary.

During orbital migration of a giant extrasolar planet via ejection of planetesimals (Murray et al.~1998), inner mean motion resonances can be strong enough to cause planetesimals to graze or impact the star. We integrate numerically the motions of particles which pass through the 3:1 or 4:1 mean motion resonances of a migrating Jupiter mass planet. We find that many particles can be trapped in the 3:1 or 4:1 resonances and pumped to high enough eccentricities that they impact the star. This implies that for a planet migrating a substantial fraction of its semi major axis, a significant fraction of its mass in planetesimals could impact the star. This process may be capable of enriching the metallicity of the star, and at a time when the star is no longer fully convective. Upon close approaches to the star the surfaces of these planetesimals will be sublimated. Orbital migration should cause continuing production of evaporating bodies, suggesting that this process should be detectable with searches for transient absorption lines in young stars. The remainder of the particles will not impact the star but can be subsequently ejected by the planet as it migrates further inwards. This allows the planet to migrate a substantial fraction of its initial semi-major axis via ejection of planetesimals.

Classical analytic theories of the solar system indicate that it is stable, but numerical integrations suggest that it is chaotic. This disagreement is resolved by a new analytic theory. The theory shows that the chaos among the Jovian planets results from the overlap of the components of a mean motion resonance among Jupiter, Saturn, and Uranus, and provides rough estimates of the Lyapunov time (10 million years) and the dynamical lifetime of Uranus (10^{18} years). The Jovian planets must have entered the resonance after all the gas and most of the planetesimals in the protoplanetary disk were removed.

A simple question of celestial mechanics is investigated: in what regions of phase space near a binary system can planets persist for long times? The planets are taken to be test particles moving in the field of an eccentric binary system. A range of values of the binary eccentricity and mass ratio is studied, and both the case of planets orbiting close to one of the stars, and that of planets outside the binary orbiting the system's center of mass, are examined. From the results, empirical expressions are developed for both 1) the largest orbit around each of the stars, and 2) the smallest orbit around the binary system as a whole, in which test particles survive the length of the integration (10^4 binary periods). The empirical expressions developed, which are roughly linear in both the mass ratio mu and the binary eccentricity e, are determined for the range 0.0 <= e <= 0.7-0.8 and 0.1 <= mu <= 0.9 in both regions, and can be used to guide searches for planets in binary systems. After considering the case of a single low-mass planet in binary systems, the stability of a mutually-interacting system of planets orbiting one star of a binary system is examined, though in less detail.

We examine the non-linear stability of the Wisdom-Holman (WH) mapping applied to the integration of perturbed, highly eccentric two-body orbits. We find that the method is unstable and introduces artificial chaos into the computed trajectories, unless the step size is chosen small enough to resolve pericenter. The origin of the instability is analyzed analytically using the Stark problem as a fiducial test case. We similarly examine the robustness of several alternative methods: a regularized WH map due to Mikkola (1997); the potential-splitting (PS) approach of Lee et al. (1997); and two types of Stark-based methods. Comparative simulations show the regularized WH map to be stable at high eccentricities, and an enhanced PS algorithm to perform well when close encounters are additionally present. We find Stark-based schemes to be of marginal use in N-body type integrations.

A planet orbiting in a disk of planetesimals can experience an instability in which it migrates to smaller orbital radii. Resonant interactions between the planet and planetesimals remove angular momentum from the planetesimals, increasing their eccentricities. Subsequently, the planetesimals either collide with or are ejected by the planet, reducing the semimajor axis of the planet. If the surface density of planetesimals exceeds a critical value, corresponding to 0.03 solar masses of gas inside the orbit of Jupiter, the planet will migrate inward a large distance. This instability may explain the presence of Jupiter-mass objects in small orbits around nearby stars.

This paper investigates the long-term orbital stability of small bodies near the central binary of the Alpha Centauri system. Test particles on circular orbits are integrated in the field of this binary for 32000 binary periods or approximately 2.5 Myr. In the region exterior to the binary, particles with semi-major axes less than roughly three times the binary's semi-major axis are unstable. Inside the binary, particles are unstable if further than 0.2 binary semimajor axes from the primary, with stablility closer in a strong function of inclination: orbits inclined near 90 degrees are unstable in as close as 0.01 binary semimajor axes from either star.