We confirm the planetary nature of TOI-1728b using a combination of ground-based photometry, near-infrared Doppler velocimetry and spectroscopy with the Habitable-zone Planet Finder.TOI-1728 is an old, inactive M0 star with \teff{} $= 3980^{+31}_{-32}$ K, which hosts a transiting super Neptune at an orbital period of $\sim$ 3.49 days. Joint fitting of the radial velocities and TESS and ground-based transits yields a planetary radius of $5.05_{-0.17}^{+0.16}$ R$_{\oplus}$, mass $26.78_{-5.13}^{+5.43}$ M$_{\oplus}$ and eccentricity $0.057_{-0.039}^{+0.054}$. We estimate the stellar properties, and perform a search for He 10830 \AA absorption during the transit of this planet and claim a null detection with an upper limit of 1.1$\%$ with 90\% confidence. A deeper level of He 10830 \AA ~ absorption has been detected in the planet atmosphere of GJ 3470b, a comparable gaseous planet. TOI-1728b is the largest super Neptune -- the intermediate subclass of planets between Neptune and the more massive gas-giant planets -- discovered around an M dwarf. With its relatively large mass and radius, TOI-1728 represents a valuable datapoint in the M-dwarf exoplanet mass-radius diagram, bridging the gap between the lighter Neptune-sized planets and the heavier Jovian planets known to orbit M-dwarfs. With a low bulk density of $1.14_{-0.24}^{+0.26}$ g/cm$^3$, and orbiting a bright host star (J $\sim 9.6$, V $\sim 12.4$), TOI-1728b is also a promising candidate for transmission spectroscopy both from the ground and from space, which can be used to constrain planet formation and evolutionary models.

We report on the validation of two planets orbiting the nearby (36pc) M2 dwarf TOI-1266 observed by the TESS mission. The inner planet is sub-Neptune-sized ($R=2.46 \pm 0.08 R_\oplus$) with an orbital period of 10.9 days. The outer planet has a radius of $1.67_{-0.11}^{+0.09} R_\oplus$ and resides in the exoplanet Radius Valley---the transition region between rocky and gaseous planets. With an orbital period of 18.8 days, the outer planet receives an insolation flux of 2.4 times that of Earth, similar to the insolation of Venus. Using precision near-infrared radial velocities with the Habitable-zone Planet Finder Spectrograph, we place upper mass limits of $15.9 M_\oplus$ and $6.4 M_\oplus$ at 95% confidence for the inner and outer planet, respectively. A more precise mass constraint of planet c, achievable with current RV instruments given the host star brightness (V=12.9, J=9.7), will yield further insights into the dominant processes sculpting the exoplanet Radius Valley.

The Kepler mission observed thousands of transiting exoplanet candidates around hundreds of thousands of FGK dwarf stars. He, Ford, & Ragozzine (2019) (arXiv:1907.07773) applied forward modelling to infer the distribution of intrinsic architectures of planetary systems, developed a clustered Poisson point process model for exoplanetary systems (SysSim) to reproduce the marginal distributions of the observed Kepler population, and showed that orbital periods and planet radii are clustered within a given planetary system. Here, we extend the clustered model to explore correlations between planetary systems and their host star properties. We split the sample of Kepler FGK dwarfs into two halves and model the fraction of stars with planets (between 0.5-10 $R_\oplus$ and 3--300 d), $f_{\rm swpa}$, as a linear function of the Gaia DR2 $b_p-r_p-E^*$ colour, where $E^*$ is a reddening correction. We find that the occurrence of these planetary systems rises significantly towards later type (redder or higher $b_p-r_p-E^*$) stars, with a slope of $df_{\rm swpa}/d(b_p-r_p-E^*) = 0.84_{-0.35}^{+0.37}$. The fraction of stars with planets increases from $f_{\rm swpa} = 0.34_{-0.12}^{+0.12}$ for F2V dwarfs to $f_{\rm swpa} = 0.92_{-0.15}^{+0.08}$ for mid K-dwarfs. About half ($f_{\rm swpa} = 0.57_{-0.10}^{+0.15}$) of all solar-type (G2V) dwarfs harbour a planetary system between 3 and 300 d. We considered three related models and find that the rise in $f_{\rm swpa}$ with $b_p-r_p-E^*$ is robust. While this linear $f_{\rm swpa}(b_p-r_p-E^*)$ model is simple, it can closely match the observed multiplicity distributions of both bluer and redder halves in our sample, suggesting that the architectures of planetary systems around stars of different spectral types may be similar aside from a shift in the overall fraction of planet hosting stars.

We present robust planet occurrence rates for Kepler planet candidates around M stars for planet radii $R_p = 0.5-4~\textrm{R}_\oplus$ and orbital periods $P = 0.5-256$ days using the approximate Bayesian computation (ABC) technique. This work incorporates the final Kepler DR25 planet candidate catalog and data products and augments them with updated stellar properties using Gaia DR2 and 2MASS PSC. We apply a set of selection criteria to select a sample of 1,746 Kepler M dwarf targets that host 89 associated planet candidates. These early M dwarfs and late K dwarfs were selected from cross-referenced targets using several photometric quality flags from Gaia DR2 and color-magnitude cuts using 2MASS magnitudes. We estimate a habitable zone occurrence rate of $f_{\textrm{M,HZ}} = 0.33^{+0.10}_{-0.12}$ for planets with $0.75-1.5$ R$_\oplus$ size. We caution that occurrence rate estimates for Kepler M stars are sensitive to the choice of prior due to the small sample of target stars and planet candidates. For example, we find an occurrence rate of $4.2^{+0.6}_{-0.6}$ or $8.4^{+1.2}_{-1.1}$ planets per M dwarf (integrating over $R_p = 0.5-4~\textrm{R}_\oplus$ and $P = 0.5-256$ days) for our two choices of prior. These occurrence rates are greater than those for FGK dwarf target when compared at the same range of orbital periods, but similar to occurrence rates when computed as a function of equivalent stellar insolation. Combining our result with recent studies of exoplanet architectures indicates that most, and potentially all, early-M dwarfs harbor planetary systems.

Young, rapidly-rotating M dwarfs exhibit prominent starspots, which create quasiperiodic signals in their photometric and Doppler spectroscopic measurements. The periodic Doppler signals can mimic radial velocity (RV) changes expected from orbiting exoplanets. Exoplanets can be distinguished from activity-induced false positives by the chromaticity and long-term incoherence of starspot signals, but these qualities are poorly constrained for fully-convective M stars. Coherent photometric starspot signals on M dwarfs may persist for hundreds of rotations, and the wavelength dependence of starspot RV signals may not be consistent between stars due to differences in their magnetic fields and active regions. We obtained precise multi-wavelength RVs of four rapidly-rotating M dwarfs (AD Leo, G 227-22, GJ 1245B, GJ 3959) using the near-infrared (NIR) Habitable-zone Planet Finder, and the optical Keck/HIRES spectrometer. Our RVs are complemented by photometry from Kepler, TESS, and the Las Cumbres Observatory (LCO) network of telescopes. We found that all four stars exhibit large spot-induced Doppler signals at their rotation periods, and investigated the longevity and optical-to-NIR chromaticity for these signals. The phase curves remain coherent much longer than is typical for Sunlike stars. Their chromaticity varies, and one star (GJ 3959) exhibits optical and NIR RV modulation consistent in both phase and amplitude. In general, though, we find that the NIR amplitudes are lower than their optical counterparts. We conclude that starspot modulation for rapidly-rotating M stars frequently remains coherent for hundreds of stellar rotations, and gives rise to Doppler signals that, due to this coherence, may be mistaken for exoplanets.

Recent and upcoming stabilized spectrographs are pushing the frontier for Doppler spectroscopy to detect and characterize low-mass planets. Specifications for these instruments are so impressive that intrinsic stellar variability is expected to limit their Doppler precision for most target stars (Fischer et al. 2016). To realize their full potential, astronomers must develop new strategies for distinguishing true Doppler shifts from intrinsic stellar variability. Stellar variability due to star spots, faculae and other rotationally-linked variability are particularly concerning, as the stellar rotation period is often included in the range of potential planet orbital periods. To robustly detect and accurately characterize low-mass planets via Doppler planet surveys, the exoplanet community must develop statistical models capable of jointly modeling planetary perturbations and intrinsic stellar variability. Towards this effort, this note presents simulations of extremely high resolution, solar-like spectra created with SOAP 2.0 (arXiv:1409.3594) that includes multiple evolving star spots. We anticipate this data set will contribute to future studies developing, testing, and comparing statistical methods for measuring physical radial velocities amid contamination by stellar variability.

Understanding the dynamics and kinematics of out-flowing atmospheres of hot and warm exoplanets is crucial to understanding the origins and evolutionary history of the exoplanets near the evaporation desert. Recently, ground based measurements of the meta-stable Helium atom's resonant absorption at 10830 \AA~has become a powerful probe of the base environment which is driving the outflow of exoplanet atmospheres. We report evidence for the He I 10830 \AA~in absorption (equivalent width $\sim$ $0.012 \pm 0.002$ \AA) in the exosphere of a warm Neptune orbiting the M-dwarf GJ 3470, during three transits using the Habitable Zone Planet Finder (HPF) near infrared spectrograph. This marks the first reported evidence for He I 10830 \AA\, atmospheric absorption for a planet orbiting an M-dwarf. Our detected absorption is broad and its blueshifted wing extends to -36 km/sec, the largest reported in the literature to date. We modelled the state of Helium atoms in the exosphere of GJ3470b based on assumptions on the UV and X-ray flux of GJ 3470, and found our measurement of flux-weighted column density of meta-stable state Helium $(N_{He^2_3S} = 2.4 \times 10^{10} \mathrm{cm^{-2}})$, derived from our transit observations, to be consistent with model, within its uncertainties. The methodology developed here will be useful to study and constrain the atmospheric outflow models of other exoplanets like GJ 3470b which are near the edge of the evaporation desert.

Planets with 2 $R_{\oplus}$ < $R$ < 3 $R_{\oplus}$ and orbital period $<$100 d are abundant; these sub-Neptune exoplanets are not well understood. For example, $Kepler$ sub-Neptunes are likely to have deep magma oceans in contact with their atmospheres, but little is known about the effect of the magma on the atmosphere. Here we study this effect using a basic model, assuming that volatiles equilibrate with magma at $T$ $\sim$ 3000 K. For our Fe-Mg-Si-O-H model system, we find that chemical reactions between the magma and the atmosphere and dissolution of volatiles into the magma are both important. Thus, magma matters. For H, most moles go into the magma, so the mass target for both H$_2$ accretion and H$_2$ loss models is weightier than is usually assumed. The known span of magma oxidation states can produce sub-Neptunes that have identical radius but with total volatile masses varying by 20-fold. Thus, planet radius is a proxy for atmospheric composition but not for total volatile content. This redox diversity degeneracy can be broken by measurements of atmosphere mean molecular weight. We emphasise H$_2$ supply by nebula gas, but also consider solid-derived H$_2$O. We find that adding H$_2$O to Fe probably cannot make enough H$_2$ to explain sub-Neptune radii because $>$10$^3$-km thick outgassed atmospheres have high mean molecular weight. The hypothesis of magma-atmosphere equilibration links observables such as atmosphere H$_2$O/H$_2$ ratio to magma FeO content and planet formation processes. Our model's accuracy is limited by the lack of experiments (lab and/or numerical) that are specific to sub-Neptunes; we advocate for such experiments.

We characterize the occurrence rate of planets, ranging in size from 0.5-16 R$_\oplus$, orbiting FGK stars with orbital periods from 0.5-500 days. Our analysis is based on results from the "DR25" catalog of planet candidates produced by NASA's Kepler mission and stellar radii from Gaia "DR2". We incorporate additional Kepler data products to accurately characterize the efficiency of planets being recognized as a "threshold crossing events" (TCE) by Kepler's Transiting Planet Search pipeline and labeled as a planet candidate by the robovetter. Using a hierarchical Bayesian model, we derive planet occurrence rates for a wide range of planet sizes and orbital periods. For planets with sizes $0.75-1.5$ R$_\oplus$ and orbital periods of 237-500 days, we find a rate of planets per FGK star of $<0.27$ ($84.13$th percentile). While the true rate of such planets could be lower by a factor of $\sim~2$ (primarily due to potential contamination of planet candidates by false alarms), the upper limits on the occurrence rate of such planets are robust to $\sim~10\%$. We recommend that mission concepts aiming to characterize potentially rocky planets in or near the habitable zone of sun-like stars prepare compelling science programs that would be robust for a true rate in the range $f_{R,P} = $ $0.03-0.40$ for $0.75-1.5$ R$_\oplus$ planets with orbital periods in 237-500 days, or a differential rate of $\Gamma_\oplus \equiv (d^2 f)/[d(\ln P)~d(\ln R_{p})] = $ $0.06-0.76$.

Transiting planets with radii 2-3 $R_\bigoplus$ are much more numerous than larger planets. We propose that this drop-off is so abrupt because at $R$ $\sim$ 3 $R_\bigoplus$, base-of-atmosphere pressure is high enough for the atmosphere to readily dissolve into magma, and this sequestration acts as a strong brake on further growth. The viability of this idea is demonstrated using a simple model. Our results support extensive magma-atmosphere equilibration on sub-Neptunes, with numerous implications for sub-Neptune formation and atmospheric chemistry.

We validate the discovery of a 2 Earth radii sub-Neptune-size planet around the nearby high proper motion M2.5-dwarf G 9-40 (EPIC 212048748), using high-precision near-infrared (NIR) radial velocity (RV) observations with the Habitable-zone Planet Finder (HPF), precision diffuser-assisted ground-based photometry with a custom narrow-band photometric filter, and adaptive optics imaging. At a distance of $d=27.9\mathrm{pc}$, G 9-40b is the second closest transiting planet discovered by K2 to date. The planet's large transit depth ($\sim$3500ppm), combined with the proximity and brightness of the host star at NIR wavelengths (J=10, K=9.2) makes G 9-40b one of the most favorable sub-Neptune-sized planet orbiting an M-dwarf for transmission spectroscopy with JWST, ARIEL, and the upcoming Extremely Large Telescopes. The star is relatively inactive with a rotation period of $\sim$29 days determined from the K2 photometry. To estimate spectroscopic stellar parameters, we describe our implementation of an empirical spectral matching algorithm using the high-resolution NIR HPF spectra. Using this algorithm, we obtain an effective temperature of $T_{\mathrm{eff}}=3404\pm73$K, and metallicity of $\mathrm{[Fe/H]}=-0.08\pm0.13$. Our RVs, when coupled with the orbital parameters derived from the transit photometry, exclude planet masses above $11.7 M_\oplus$ with 99.7% confidence assuming a circular orbit. From its radius, we predict a mass of $M=5.0^{+3.8}_{-1.9} M_\oplus$ and an RV semi-amplitude of $K=4.1^{+3.1}_{-1.6}\mathrm{m\:s^{-1}}$, making its mass measurable with current RV facilities. We urge further RV follow-up observations to precisely measure its mass, to enable precise transmission spectroscopic measurements in the future.

Observations of exoplanetary systems provide clues about the intrinsic distribution of planetary systems, their architectures, and how they formed. We develop a forward modelling framework for generating populations of planetary systems and "observed" catalogues by simulating the Kepler detection pipeline (SysSim). We compare our simulated catalogues to the Kepler DR25 catalogue of planet candidates, updated to include revised stellar radii from Gaia DR2. We constrain our model based on the observed 1D marginal distributions of orbital periods, period ratios, transit depths, transit depth ratios, transit durations, transit duration ratios, and transit multiplicities. Models assuming planets with independent periods and sizes do not adequately account for the properties of the multiplanet systems. Instead, a clustered point process model for exoplanet periods and sizes provides a significantly better description of the Kepler population, particularly the observed multiplicity and period ratio distributions. We find that $0.56^{+0.18}_{-0.15}$ of FGK stars have at least one planet larger than $0.5 R_\oplus$ between 3 and 300 d. Most of these planetary systems ($\sim 98\%$) consist of one or two clusters with a median of three planets per cluster. We find that the Kepler dichotomy is evidence for a population of highly inclined planetary systems and is unlikely to be solely due to a population of intrinsically single planet systems. We provide a large ensemble of simulated physical and observed catalogues of planetary systems from our models, as well as publicly available code for generating similar catalogues given user-defined parameters.

We present results from a data challenge posed to the radial velocity (RV) community: namely, to quantify the Bayesian "evidence" for n={0,1,2,3} planets in a set of synthetically generated RV datasets containing a range of planet signals. Participating teams were provided the same likelihood function and set of priors to use in their analysis. They applied a variety of methods to estimate Z, the marginal likelihood for each n-planet model, including cross-validation, the Laplace approximation, importance sampling, and nested sampling. We found the dispersion in Z across different methods grew with increasing n-planet models: ~3 for 0-planets, ~10 for 1-planet, ~100-1000 for 2-planets, and >10,000 for 3-planets. Most internal estimates of uncertainty in Z for individual methods significantly underestimated the observed dispersion across all methods. Methods that adopted a Monte Carlo approach by comparing estimates from multiple runs yielded plausible uncertainties. Finally, two classes of numerical algorithms (those based on importance and nested samplers) arrived at similar conclusions regarding the ratio of Zs for n and (n+1)-planet models. One analytic method (the Laplace approximation) demonstrated comparable performance. We express both optimism and caution: we demonstrate that it is practical to perform rigorous Bayesian model comparison for <=3-planet models, yet robust planet discoveries require researchers to better understand the uncertainty in Z and its connections to model selection.

This Astro2020 State of the Profession Consideration White Paper highlights the growth of astrostatistics and astroinformatics in astronomy, identifies key issues hampering the maturation of these new subfields, and makes recommendations for structural improvements at different levels that, if acted upon, will make significant positive impacts across astronomy.

To date a dozen transiting "Tatooines" or circumbinary planets (CBPs) have been discovered, by eye, in the data from the Kepler mission; by contrast, thousands of confirmed circumstellar planets orbiting around single stars have been detected using automated algorithms. Automated detection of CBPs is challenging because their transits are strongly aperiodic with irregular profiles. Here, we describe an efficient and automated technique for detecting circumbinary planets that transit their binary hosts in Kepler light curves. Our method accounts for large transit timing and duration variations (TTVs and TDVs), induced by binary reflex motion, in two ways: 1) We directly correct for large-scale TTVs and TDVs in the light curves by using Keplerian models to approximate binary and CBP orbits; and 2) We allow additional aperiodicities on the corrected light curves by employing the Quasi-periodic Automated Transit Search algorithm (QATS). We demonstrate that our method dramatically improves detection significance using simulated data and two previously identified CBP systems, Kepler-35 and Kepler-64.

We infer the number of planets-per-star as a function of orbital period and planet size using $Kepler$ archival data products with updated stellar properties from the $Gaia$ Data Release 2. Using hierarchical Bayesian modeling and Hamiltonian Monte Carlo, we incorporate planet radius uncertainties into an inhomogeneous Poisson point process model. We demonstrate that this model captures the general features of the outcome of the planet formation process around GK stars, and provides an infrastructure to use the $Kepler$ results to constrain analytic planet distribution models. We report an increased mean and variance in the marginal posterior distributions for the number of planets per $GK$ star when including planet radius measurement uncertainties. We estimate the number of planets-per-$GK$ star between 0.75 and 2.5 $R_{\oplus}$ and 50 to 300 day orbital periods to have a $68\%$ credible interval of $0.49$ to $0.77$ and a posterior mean of $0.63$. This posterior has a smaller mean and a larger variance than the occurrence rate calculated in this work and in Burke et al. (2015) for the same parameter space using the $Q1-Q16$ (previous $Kepler$ planet candidate and stellar catalog), and a larger mean and variance than when using the $DR25$ (latest $Kepler$ planet candidate and stellar catalog). We find that the accuracy and precision of our hierarchical Bayesian model posterior distributions are less sensitive to the total number of planets in the sample, and more so on the characteristics of the catalog completeness and reliability and the span of the planet parameter space.

The 2020s will be the most data-rich decade of astronomy in history. As the scale and complexity of our surveys increase, the problem of scheduling becomes more critical. We must develop high-quality scheduling approaches, implement them as open-source software, and begin linking the typically separate stages of observation and data analysis.

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

Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New methodologies derived from advances in statistics, computer science, and machine learning are beginning to be employed in sophisticated investigations that are not only bringing forth new discoveries, but are placing them on a solid footing. Progress in wide-field sky surveys, interferometric imaging, precision cosmology, exoplanet detection and characterization, and many subfields of stellar, Galactic and extragalactic astronomy, has resulted in complex data analysis challenges that must be solved to perform scientific inference. Research in astrostatistics and astroinformatics will be necessary to develop the state-of-the-art methodology needed in astronomy. Overcoming these challenges requires dedicated, interdisciplinary research. We recommend: (1) increasing funding for interdisciplinary projects in astrostatistics and astroinformatics; (2) dedicating space and time at conferences for interdisciplinary research and promotion; (3) developing sustainable funding for long-term astrostatisics appointments; and (4) funding infrastructure development for data archives and archive support, state-of-the-art algorithms, and efficient computing.

Protoplanetary disks are thought to be truncated at orbital periods of around 10 days. Therefore, origin of rocky short period planets with $P < 10$ days is a puzzle. We propose that many of these planets may form through the Type-I migration of planets locked into a chain of mutual mean motion resonances. We ran N-body simulations of planetary embryos embedded in a protoplanetary disk. The embryos experienced gravitational scatterings, collisions, disk torques, and dampening of orbital eccentricity and inclination. We then modelled Kepler observations of these planets using a forward model of both the transit probability and the detection efficiency of the Kepler pipeline. We found that planets become locked into long chains of mean motion resonances that migrate in unison. When the chain reaches the edge of the disk, the inner planets are pushed past the edge due to the disk torques acting on the planets farther out in the chain. Our simulated systems successfully reproduce the observed period distribution of short period Kepler planets between 1 and 2 $R_\oplus$. However, we obtain fewer closely packed short period planets than in the Kepler sample. Our results provide valuable insight into the planet formation process, and suggests that resonance locks, migration, and dynamical instabilities play important roles the the formation and evolution of close-in small exoplanets.

Transit timing variations (TTVs) are a valuable tool to determine the masses and orbits of transiting planets in multi-planet systems. TTVs can be readily modeled given knowledge of the interacting planets' orbital configurations and planet-star mass ratios, but such models are highly nonlinear and difficult to invert. Markov chain Monte Carlo (MCMC) methods are often used to explore the posterior distribution for model parameters, but, due to the high correlations between parameters, nonlinearity, and potential multi-modality in the posterior, many samplers perform very inefficiently. Therefore, we assess the performance of several MCMC samplers that use varying degrees of geometric information about the target distribution. We generate synthetic datasets from multiple models, including the TTVFaster model and a simple sinusoidal model, and test the efficiencies of various MCMC samplers. We find that sampling efficiency can be greatly improved for all models by sampling from a parameter space transformed using an estimate of the covariance and means of the target distribution. No one sampler performs the best for all datasets, but several samplers, such as Differential Evolution Monte Carlo and Geometric adaptive Monte Carlo, have consistently efficient performance. For datasets with near Gaussian posteriors, Hamiltonian Monte Carlo samplers with 2 or 3 leapfrog steps obtained the highest efficiencies. Based on differences in effective sample sizes per time, we show that the right choice of sampler can improve sampling efficiencies by several orders of magnitude.

We present empirical evidence, supported by a planet formation model, to show that the curve $R/R_\oplus = 1.05\,(F/F_\oplus)^{0.11}$ approximates the location of the so-called photo-evaporation valley. Planets below that curve are likely to have experienced complete photo-evaporation, and planets just above it appear to have inflated radii; thus we identify a new population of inflated super-Earths and mini-Neptunes. Our N-body simulations are set within an evolving protoplanetary disk and include prescriptions for orbital migration, gas accretion, and atmospheric loss due to giant impacts. Our simulated systems broadly match the sizes and periods of super-Earths in the Kepler catalog. They also reproduce the relative sizes of adjacent planets in the same system, with the exception of planet pairs that straddle the photo-evaporation valley. This latter group is populated by planet pairs with either very large or very small size ratios ($R_{\rm out} / R_{\rm in} \gg 1$ or $R_{\rm out} / R_{\rm in} \ll 1$) and a dearth of size ratios near unity. It appears that this feature could be reproduced if the planet outside the photo-evaporation valley (typically the outer planet, but some times not) has its atmosphere significantly expanded by stellar irradiation. This new population of planets may be ideal targets for future transit spectroscopy observations with the upcoming James Webb Space Telescope.

Many habitable zone exoplanets are expected to form with water mass fractions higher than that of the Earth. For rocky exoplanets with 10-1000x Earth's H2O but without H2, we model the multi-Gyr evolution of ocean temperature and chemistry, taking into account C partitioning, high-pressure ice phases, and atmosphere-lithosphere exchange. Within our model, for Sun-like stars, we find that: (1)~the duration of habitable surface water is strongly affected by ocean chemistry; (2)~possible ocean pH spans a wide range; (3)~surprisingly, many waterworlds retain habitable surface water for >1 Gyr, and (contrary to previous claims) this longevity does not necessarily involve geochemical cycling. The key to this cycle-independent planetary habitability is that C exchange between the convecting mantle and the water ocean is curtailed by seafloor pressure on waterworlds, so the planet is stuck with the ocean mass and ocean cations that it acquires during the first 1% of its history. In our model, the sum of positive charges leached from the planetary crust by early water-rock interactions is - coincidentally - often within an order of magnitude of the early-acquired atmosphere+ocean inorganic C inventory overlaps. As a result, pCO2 is frequently in the "sweet spot" (0.2-20 bar) for which the range of semimajor axis that permits surface liquid water is about as wide as it can be. Because the width of the HZ in semimajor axis defines (for Sun-like stars) the maximum possible time span of surface habitability, this effect allows for Gyr of habitability as the star brightens. We illustrate our findings by using the output of an ensemble of N-body simulations as input to our waterworld evolution code. Thus (for the first time in an end-to-end calculation) we show that chance variation of initial conditions, with no need for geochemical cycling, can yield multi-Gyr surface habitability on waterworlds.

We present a new framework to characterize the occurrence rates of planet candidates identified by Kepler based on hierarchical Bayesian modeling, Approximate Bayesian Computing (ABC), and sequential importance sampling. For this study we adopt a simple 2-D grid in planet radius and orbital period as our model and apply our algorithm to estimate occurrence rates for Q1-Q16 planet candidates orbiting around solar-type stars. We arrive at significantly increased planet occurrence rates for small planet candidates ($R_p<1.25 R_{\oplus}$) at larger orbital periods ($P>80$d) compared to the rates estimated by the more common inverse detection efficiency method. Our improved methodology estimates that the occurrence rate density of small planet candidates in the habitable zone of solar-type stars is $1.6^{+1.2}_{-0.5}$ per factor of 2 in planet radius and orbital period. Additionally, we observe a local minimum in the occurrence rate for strong planet candidates marginalized over orbital period between 1.5 and 2$R_{\oplus}$ that is consistent with previous studies. For future improvements, the forward modeling approach of ABC is ideally suited to incorporating multiple populations, such as planets, astrophysical false positives and pipeline false alarms, to provide accurate planet occurrence rates and uncertainties. Furthermore, ABC provides a practical statistical framework for answering complex questions (e.g., frequency of different planetary architectures) and providing sound uncertainties, even in the face of complex selection effects, observational biases, and follow-up strategies. In summary, ABC offers a powerful tool for accurately characterizing a wide variety of astrophysical populations.

We outline the important role that ground-based, Doppler monitoring of exoplanetary systems will play in advancing our theories of planet formation and dynamical evolution. A census of planetary systems requires a well designed survey to be executed over the course of a decade or longer. A coordinated survey to monitor several thousand targets each at ~1000 epochs (~3-5 million new observations) will require roughly 40 dedicated spectrographs. We advocate for improvements in data management, data sharing, analysis techniques, and software testing, as well as possible changes to the funding structures for exoplanet science.

This whitepaper discusses the diversity of exoplanets that could be detected by future observations, so that comparative exoplanetology can be performed in the upcoming era of large space-based flagship missions. The primary focus will be on characterizing Earth-like worlds around Sun-like stars. However, we will also be able to characterize companion planets in the system simultaneously. This will not only provide a contextual picture with regards to our Solar system, but also presents a unique opportunity to observe size dependent planetary atmospheres at different orbital distances. We propose a preliminary scheme based on chemical behavior of gases and condensates in a planet's atmosphere that classifies them with respect to planetary radius and incident stellar flux.

We present the Kepler Object of Interest (KOI) catalog of transiting exoplanets based on searching four years of Kepler time series photometry (Data Release 25, Q1-Q17). The catalog contains 8054 KOIs of which 4034 are planet candidates with periods between 0.25 and 632 days. Of these candidates, 219 are new and include two in multi-planet systems (KOI-82.06 and KOI-2926.05), and ten high-reliability, terrestrial-size, habitable zone candidates. This catalog was created using a tool called the Robovetter which automatically vets the DR25 Threshold Crossing Events (TCEs, Twicken et al. 2016). The Robovetter also vetted simulated data sets and measured how well it was able to separate TCEs caused by noise from those caused by low signal-to-noise transits. We discusses the Robovetter and the metrics it uses to sort TCEs. For orbital periods less than 100 days the Robovetter completeness (the fraction of simulated transits that are determined to be planet candidates) across all observed stars is greater than 85%. For the same period range, the catalog reliability (the fraction of candidates that are not due to instrumental or stellar noise) is greater than 98%. However, for low signal-to-noise candidates between 200 and 500 days around FGK dwarf stars, the Robovetter is 76.7% complete and the catalog is 50.5% reliable. The KOI catalog, the transit fits and all of the simulated data used to characterize this catalog are available at the NASA Exoplanet Archive.

The radial velocity method is one of the most successful techniques for detecting exoplanets. It works by detecting the velocity of a host star induced by the gravitational effect of an orbiting planet, specifically the velocity along our line of sight, which is called the radial velocity of the star. As astronomical instrumentation has improved, radial velocity surveys have become sensitive to low-mass planets that cause their host star to move with radial velocities of 1 m/s or less. While analysis of a time series of stellar spectra can in theory reveal such small radial velocities, in practice intrinsic stellar variability (e.g., star spots, convective motion, pulsations) affects the spectra and often mimics a radial velocity signal. This signal contamination makes it difficult to reliably detect low mass planets and planets orbiting magnetically active stars. A principled approach to recovering planet radial velocity signals in the presence of stellar activity was proposed by Rajpaul et al. (2015) and involves the use of a multivariate Gaussian process model to jointly capture time series of the apparent radial velocity and multiple indicators of stellar activity. We build on this work in two ways: (i) we propose using dimension reduction techniques to construct more informative stellar activity indicators that make use of a larger portion of the stellar spectrum; (ii) we extend the Rajpaul et al. (2015) model to a larger class of models and use a model comparison procedure to select the best model for the particular stellar activity indicators at hand. By combining our high-information stellar activity indicators, Gaussian process models, and model selection procedure, we achieve substantially improved planet detection power compared to previous state-of-the-art approaches.

The NASA Exoplanet Exploration Program's SAG15 group has solicited, collected, and organized community input on high-level science questions that could be addressed with future direct imaging exoplanet missions and the type and quality of data answering these questions will require. Input was solicited through a variety of forums and the report draft was shared with the exoplanet community continuously during the period of the report development (Nov 2015 -- May 2017). The report benefitted from the input of over 50 exoplanet scientists and from multiple open-forum discussions at exoplanet and astrobiology meetings. The SAG15 team has identified three group of high-level questions, those that focus on the properties of planetary systems (Questions A1--A2), those that focus on the properties of individual planets (Questions B1--B4), and questions that relate to planetary processes (Questions C1--C4). The questions in categories A, B, and C require different target samples and often different observational approaches. For each questions we summarize the current body of knowledge, the available and future observational approaches that can directly or indirectly contribute to answering the question, and provide examples and general considerations for the target sample required.

Photospheric velocities and stellar activity features such as spots and faculae produce measurable radial velocity signals that currently obscure the detection of sub-meter-per-second planetary signals. However, photospheric velocities are imprinted differently in a high-resolution spectrum than Keplerian Doppler shifts. Photospheric activity produces subtle differences in the shapes of absorption lines due to differences in how temperature or pressure affects the atomic transitions. In contrast, Keplerian Doppler shifts affect every spectral line in the same way. With high enough S/N and high enough resolution, statistical techniques can exploit differences in spectra to disentangle the photospheric velocities and detect lower-amplitude exoplanet signals. We use simulated disk-integrated time-series spectra and principal component analysis (PCA) to show that photospheric signals introduce spectral line variability that is distinct from Doppler shifts. We quantify the impact of instrumental resolution and S/N for this work.

Disk migration and high-eccentricity migration are two well-studied theories to explain the formation of hot Jupiters. The former predicts that these planets can migrate up until the planet-star Roche separation ($a_{Roche}$) and the latter predicts they will tidally circularize at a minimum distance of 2$a_{Roche}$. Considering long-running radial velocity and transit surveys have identified a couple hundred hot Jupiters to date, we can revisit the classic question of hot Jupiter formation in a data-driven manner. We approach this problem using data from several exoplanet surveys (radial velocity, Kepler, HAT, and WASP) allowing for either a single population or a mixture of populations associated with these formation channels, and applying a hierarchical Bayesian mixture model of truncated power laws of the form $x^{\gamma-1}$ to constrain the population-level parameters of interest (e.g., location of inner edges, $\gamma$, mixture fractions). Within the limitations of our chosen models, we find the current radial velocity and Kepler sample of hot Jupiters can be well explained with a single truncated power law distribution with a lower cutoff near 2$a_{Roche}$, a result that still holds after a decade, and $\gamma=-0.51\pm^{0.19}_{0.20}$. However, the HAT and WASP data show evidence for multiple populations (Bayes factor $\approx 10^{21}$). We find that $15\pm^{9}_{6}\%$ reside in a component consistent with disk migration ($\gamma=-0.04\pm^{0.53}_{1.27}$) and $85\pm^{6}_{9}\%$ in one consistent with high-eccentricity migration ($\gamma=-1.38\pm^{0.32}_{0.47}$). We find no immediately strong connections with some observed host star properties and speculate on how future exoplanet surveys could improve upon hot Jupiter population inference.

The outer architectures of Kepler's compact systems of multiple transiting planets remain poorly constrained and few of these systems have lower bounds on the orbital distance of any massive outer planets. We infer a minimum orbital distance and upper limits on the inclination of a hypothetical Jovian-mass planet orbiting exterior to the six transiting planets at Kepler-11. Our constraints are derived from dynamical models together with observations provided by the Kepler mission. Firstly, the lack of transit timing variations (TTV) in the outermost transiting planet Kepler-11 g imply that the system does not contain a Jovian mass perturber within 2 AU from the star. Secondly, we test under what initial conditions a Jovian-mass planet moderately inclined from the transiting planets would make their co-transiting configuration unlikely. The transiting planets are secularly coupled and exhibit small mutual inclinations over long time-scales, although the outermost transiting planet, Kepler-11 g is weakly coupled to the inner five. We rule out a Jovian-mass planet on a 3 degree inclination within 3.0 AU, and higher inclinations out to further orbital distances, unless an undetected planet exists orbiting in the dynamical gap between Kepler-11 f and Kepler-11 g. Our constraints depend little on whether we assume the six transiting planets of Kepler-11 were initially perfectly co-planar or whether a minimum initial mutual inclination between the transiting planets is adopted based on the measured impact parameters of the transiting planets.

Surveys have revealed many multi-planet systems containing super-Earths and Neptunes in orbits of a few days to a few months. There is debate whether in situ assembly or inward migration is the dominant mechanism of the formation of such planetary systems. Simulations suggest that migration creates tightly packed systems with planets whose orbital periods may be expressed as ratios of small integers (resonances), often in a many-planet series (chain). In the hundreds of multi-planet systems of sub-Neptunes, more planet pairs are observed near resonances than would generally be expected, but no individual system has hitherto been identified that must have been formed by migration. Proximity to resonance enables the detection of planets perturbing each other. Here we report transit timing variations of the four planets in the Kepler-223 system, model these variations as resonant-angle librations, and compute the long-term stability of the resonant chain. The architecture of Kepler-223 is too finely tuned to have been formed by scattering, and our numerical simulations demonstrate that its properties are natural outcomes of the migration hypothesis. Similar systems could be destabilized by any of several mechanisms, contributing to the observed orbital-period distribution, where many planets are not in resonances. Planetesimal interactions in particular are thought to be responsible for establishing the current orbits of the four giant planets in the Solar System by disrupting a theoretical initial resonant chain similar to that observed in Kepler-223.

We present ALMA Band 6 observations (1.3 mm/233 GHz) of Fomalhaut and its debris disc. The observations achieve a sensitivity of 17 $\mu$Jy and a resolution of 0.28 arcsec (2.1 au at a distance of 7.66 pc), which are the highest resolution observations to date of the millimetre grains in Fomalhaut's main debris ring. The ring is tightly constrained to $139^{+2}_{-3}$ au with a FWHM of $13\pm3$ au, following a Gaussian profile. The millimetre spectral index is constrained to $\alpha_{mm} = -2.62\pm0.12$. We explore fitting debris disc models in the image plane, as well as fitting models using visibility data directly. The results are compared and the potential advantages/disadvantages of each approach are discussed. The detected central emission is indistinguishable from a point source, with a most probable flux of $0.90\pm 0.12$ mJy (including calibration uncertainties). This implies that any inner debris structure, as was inferred from far-Infrared observations, must contribute little to the total central emission. Moreover, the stellar flux is less than 70\% of that predicted by extrapolating a black body from the constrained stellar photosphere temperature. This result emphasizes that unresolved inner debris components cannot be fully characterized until the behaviour of the host star's intrinsic stellar emission at millimetre wavelengths is properly understood.

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 present the discovery of three modestly-irradiated, roughly Neptune-mass planets orbiting three nearby Solar-type stars. HD 42618 b has a minimum mass of $15.4 \pm 2.4$ M$_{\oplus}$, a semi-major axis of 0.55 AU, an equilibrium temperature of 337 K, and is the first planet discovered to orbit the solar analogue host star, HD 42618. We also discover new planets orbiting the known exoplanet host stars HD 164922 and HD 143761 ($\rho$ CrB). The new planet orbiting HD 164922 has a minimum mass of $12.9 \pm 1.6$ M$_{\oplus}$ and orbits interior to the previously known Jovian mass planet orbiting at 2.1 AU. HD 164922 c has a semi-major axis of 0.34 AU and an equilibrium temperature of 418 K. HD 143761 c orbits with a semi-major axis of 0.44 AU, has a minimum mass of $25 \pm 2$ M$_{\oplus}$, and is the warmest of the three new planets with an equilibrium temperature of 445 K. It orbits exterior to the previously known warm Jupiter in the system. A transit search using space-based CoRoT data and ground-based photometry from the Automated Photometric Telescopes (APTs) at Fairborn Observatory failed to detect any transits, but the precise, high-cadence APT photometry helped to disentangle planetary-reflex motion from stellar activity. These planets were discovered as part of an ongoing radial velocity survey of bright, nearby, chromospherically-inactive stars using the Automated Planet Finder (APF) telescope at Lick Observatory. The high-cadence APF data combined with nearly two decades of radial velocity data from Keck Observatory and gives unprecedented sensitivity to both short-period low-mass, and long-period intermediate-mass planets.

We present a new transit timing catalog of 2599 Kepler Objects of Interest (=KOIs), using the PDC-MAP long-cadence light curves that include the full seventeen quarters of the mission (ftp://wise- ftp.tau.ac.il/pub/tauttv/TTV/ver_112). The goal is to produce an easy-to-use catalog that can stimulate further analyses of interesting systems. For 779 KOIs with high enough SNRs, we derived the timing, duration and depth of 69,914 transits. For 1820 KOIs with lower SNR, we derived only the timing of 225,273 transits. After removal of outlier timings, we derived various statistics for each KOI that were used to indicate significant variations. Including systems found by previous works, we have detected 260 KOIs which showed significant TTVs with long-term variations (>100 day), and another fourteen KOIs with periodic modulations shorter than 100 day and small amplitudes. For five of those, the periodicity is probably due to the crossing of rotating stellar spots by the transiting planets.

We present ALMA band 7 (345 GHz) continuum and $^{12}$CO(J = 3-2) observations of the circumstellar disk surrounding HD141569. At an age of about 5 Myr, the disk has a complex morphology that may be best interpreted as a nascent debris system with gas. Our $870\rm~\mu m$ ALMA continuum observations resolve a dust disk out to approximately $ 56 ~\rm au$ from the star (assuming a distance of 116 pc) with $0."38$ resolution and $0.07 ~ \rm mJy~beam^{-1}$ sensitivity. We measure a continuum flux density for this inner material of $3.8 \pm 0.4 ~ \rm mJy$ (including calibration uncertainties). The $^{12}$CO(3-2) gas is resolved kinematically and spatially from about 30 to 210 au. The integrated $^{12}$CO(3-2) line flux density is $15.7 \pm 1.6~\rm Jy~km~s^{-1}$. We estimate the mass of the millimeter debris and $^{12}$CO(3-2) gas to be $\gtrsim0.04~\rm M_{\oplus}$ and $\sim2\times 10^{-3}~\rm M_{\oplus}$, respectively. If the millimeter grains are part of a collisional cascade, then we infer that the inner disk ($<50$ au) has $\sim 160~\rm M_{\oplus}$ contained within objects less than 50 km in radius, depending on the planetesimal size distribution and density assumptions. MCMC modeling of the system reveals a disk morphology with an inclination of $53.4^{\circ}$ centered around a $\rm M=2.39~ M_{\odot}$ host star ($\rm Msin(i)=1.92~ M_{\odot}$). We discuss whether the gas in HD141569's disk may be second generation. If it is, the system can be used to study the clearing stages of planet formation.

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

The Kepler Mission has discovered thousands of planets with radii $<4\ R_\oplus$, paving the way for the first statistical studies of the dynamics, formation, and evolution of these sub-Neptunes and super-Earths. Planetary masses are an important physical property for these studies, and yet the vast majority of Kepler planet candidates do not have theirs measured. A key concern is therefore how to map the measured radii to mass estimates in this Earth-to-Neptune size range where there are no Solar System analogs. Previous works have derived deterministic, one-to-one relationships between radius and mass. However, if these planets span a range of compositions as expected, then an intrinsic scatter about this relationship must exist in the population. Here we present the first probabilistic mass-radius relationship (M-R relation) evaluated within a Bayesian framework, which both quantifies this intrinsic dispersion and the uncertainties on the M-R relation parameters. We analyze how the results depend on the radius range of the sample, and on how the masses were measured. Assuming that the M-R relation can be described as a power law with a dispersion that is constant and normally distributed, we find that $M/M_\oplus=2.7(R/R_\oplus)^{1.3}$, a scatter in mass of $1.9\ M_\oplus$, and a mass constraint to physically plausible densities, is the "best-fit" probabilistic M-R relation for the sample of RV-measured transiting sub-Neptunes ($R_{pl}<4\ R_\oplus$). More broadly, this work provides a framework for further analyses of the M-R relation and its probable dependencies on period and stellar properties.

We infer dynamical masses in eight multi-planet systems using transit times measured from Kepler's complete dataset, including short-cadence data where available. Of the eighteen dynamical masses that we infer, ten pass multiple tests for robustness. These are in systems; Kepler-26 (KOI-250), Kepler-29 (KOI-738), Kepler-60 (KOI-2086), Kepler-105 (KOI-115), and Kepler-307 (KOI-1576). Kepler-105 c has a density consistent with an Earth-like composition. Strong TTV signals were detected from additional planets, but their inferred masses were sensitive to outliers or consistent solutions could not be found with independently-measured transit times, including planets at; Kepler-49 (KOI-248), Kepler-57 (KOI-1270), Kepler-105 (KOI-115) and Kepler-177 (KOI-523). Strong upper limits on the mass of Kepler-177 c imply an extremely low density ~0.1 g cm$^{-3}$. In most cases, individual orbital eccentricities were poorly constrained due to degeneracies in TTV inversion. For five planet pairs in our sample, strong secular interactions imply a moderate-to-high likelihood of apsidal alignment over a wide range of possible eccentricities. We also find solutions for the three planets known to orbit Kepler-60 in a Laplace-like resonance chain. However, non-librating solutions also match the transit-timing data. For six systems, we calculate more precise stellar parameters than previously known, enabling useful constraints on planetary densities where we have robust mass measurements. Placing these exoplanets on the mass-radius diagram, we find a wide range of densities is observed among sub-Neptune mass planets and that the range in observed densities is anti-correlated with incident flux.

The Second Workshop on Extreme Precision Radial Velocities defined circa 2015 the state of the art Doppler precision and identified the critical path challenges for reaching 10 cm/s measurement precision. The presentations and discussion of key issues for instrumentation and data analysis and the workshop recommendations for achieving this precision are summarized here. Beginning with the HARPS spectrograph, technological advances for precision radial velocity measurements have focused on building extremely stable instruments. To reach still higher precision, future spectrometers will need to produce even higher fidelity spectra. This should be possible with improved environmental control, greater stability in the illumination of the spectrometer optics, better detectors, more precise wavelength calibration, and broader bandwidth spectra. Key data analysis challenges for the precision radial velocity community include distinguishing center of mass Keplerian motion from photospheric velocities, and the proper treatment of telluric contamination. Success here is coupled to the instrument design, but also requires the implementation of robust statistical and modeling techniques. Center of mass velocities produce Doppler shifts that affect every line identically, while photospheric velocities produce line profile asymmetries with wavelength and temporal dependencies that are different from Keplerian signals. Exoplanets are an important subfield of astronomy and there has been an impressive rate of discovery over the past two decades. Higher precision radial velocity measurements are required to serve as a discovery technique for potentially habitable worlds and to characterize detections from transit missions. The future of exoplanet science has very different trajectories depending on the precision that can ultimately be achieved with Doppler measurements.

Arnold (2005), Forgan (2013), and Korpela et al. (2015) noted that planet-sized artificial structures could be discovered with Kepler as they transit their host star. We present a general discussion of transiting megastructures, and enumerate ten potential ways their anomalous silhouettes, orbits, and transmission properties would distinguish them from exoplanets. We also enumerate the natural sources of such signatures. Several anomalous objects, such as KIC 12557548 and CoRoT-29, have variability in depth consistent with Arnold's prediction and/or an asymmetric shape consistent with Forgan's model. Since well motivated physical models have so far provided natural explanations for these signals, the ETI hypothesis is not warranted for these objects, but they still serve as useful examples of how nonstandard transit signatures might be identified and interpreted in a SETI context. Boyajian et al. 2015 recently announced KIC 8462852, an object with a bizarre light curve consistent with a "swarm" of megastructures. We suggest this is an outstanding SETI target. We develop the normalized information content statistic $M$ to quantify the information content in a signal embedded in a discrete series of bounded measurements, such as variable transit depths, and show that it can be used to distinguish among constant sources, interstellar beacons, and naturally stochastic or artificial, information-rich signals. We apply this formalism to KIC 12557548 and a specific form of beacon suggested by Arnold to illustrate its utility.

We characterize the eccentricity distribution of a sample of ~50 short-period planet candidates using transit and occultation measurements from NASA's Kepler Mission. First, we evaluate the sensitivity of our hierarchical Bayesian modeling and test its robustness to model misspecification using simulated data. When analyzing actual data assuming a Rayleigh distribution for eccentricity, we find that the posterior mode for the dispersion parameter is $\sigma=0.081 \pm^{0.014}_{0.003}$. We find that a two-component Gaussian mixture model for $e \cos \omega$ and $e \sin \omega$ provides a better model than either a Rayleigh or Beta distribution. Based on our favored model, we find that $\sim90\%$ of planet candidates in our sample come from a population with an eccentricity distribution characterized by a small dispersion ($\sim0.01$), and $\sim10\%$ come from a population with a larger dispersion ($\sim0.22$). Finally, we investigate how the eccentricity distribution correlates with selected planet and host star parameters. We find evidence that suggests systems around higher metallicity stars and planet candidates with smaller radii come from a more complex eccentricity distribution.

An intriguing trend among Kepler's multi-planet systems is an overabundance of planet pairs with period ratios just wide of mean motion resonances (MMR) and a dearth of systems just narrow of them. In a recently published paper Chatterjee & Ford (2015; henceforth CF15) has proposed that gas-disk migration traps planets in a MMR. After gas dispersal, orbits of these trapped planets are altered through interaction with a residual planetesimal disk. They found that for massive enough disks planet-planetesimal disk interactions can break resonances and naturally create moderate to large positive offsets from the initial period ratio for large ranges of planetesimal disk and planet properties. Divergence from resonance only happens if the mass of planetesimals that interact with the planets is at least a few percent of the total planet mass. This threshold, above which resonances are broken and the offset from resonances can grow, naturally explains why the asymmetric large offsets were not seen in more massive planet pairs found via past radial velocity surveys. In this article we will highlight some of the key findings of CF15. In addition, we report preliminary results from an extension of this study, that investigates the effects of planet-planetesimal disk interactions on initially non-resonant planet pairs. We find that planetesimal scattering typically increases period ratios of non-resonant planets. If the initial period ratios are below and in proximity of a resonance, under certain conditions, this increment in period ratios can create a deficit of systems with period ratios just below the exact integer corresponding to the MMR and an excess just above. From an initially uniform distribution of period ratios just below a 2:1 MMR, planetesimal interactions can create an asymmetric distribution across this MMR similar to what is observed for the Kepler planet pairs.

We report constraints on the three-dimensional orbital architecture for all four planets known to orbit the nearby M dwarf Gliese 876 based solely on Doppler measurements and demanding long-term orbital stability. Our dataset incorporates publicly available radial velocities taken with the ELODIE and CORALIE spectrographs, HARPS, and Keck HIRES as well as previously unpublished HIRES velocities. We first quantitatively assess the validity of the planets thought to orbit GJ 876 by computing the Bayes factors for a variety of different coplanar models using an importance sampling algorithm. We find that a four-planet model is preferred over a three-planet model. Next, we apply a Newtonian MCMC algorithm to perform a Bayesian analysis of the planet masses and orbits using an n-body model in three-dimensional space. Based on the radial velocities alone, we find that a 99% credible interval provides upper limits on the mutual inclinations for the three resonant planets ($\Phi_{cb}<6.20^\circ$ for the "c" and "b" pair and $\Phi_{be}<28.5^\circ$ for the "b" and "e" pair). Subsequent dynamical integrations of our posterior sample find that the GJ 876 planets must be roughly coplanar ($\Phi_{cb}<2.60^\circ$ and $\Phi_{be}<7.87^\circ$), suggesting the amount of planet-planet scattering in the system has been low. We investigate the distribution of the respective resonant arguments of each planet pair and find that at least one argument for each planet pair and the Laplace argument librate. The libration amplitudes in our three-dimensional orbital model supports the idea of the outer-three planets having undergone significant past disk migration.

Extrasolar planets that pass in front of their host star (transit) cause a temporary decrease in the apparent brightness of the star once per orbit, providing a direct measure of the planet's size and orbital period. In some systems with multiple transiting planets, the times of the transits are measurably affected by the gravitational interactions between neighbouring planets. In favorable cases, the departures from Keplerian orbits implied by the observed transit times permit planetary masses to be measured, which is key to determining bulk densities. Characterizing rocky planets is particularly difficult, since they are generally smaller and less massive than gaseous planets. Thus, few exoplanets near Earth's size have had their masses measured. Here we report the sizes and masses of three planets orbiting Kepler-138, a star much fainter and cooler than the Sun. We measure the mass of the Mars-sized inner planet based on on the transit times of its neighbour and thereby provide the first density measurement for an exoplanet smaller than Earth. The middle and outer planets are both slightly larger than Earth. The middle planet's density is similar to that of Earth, while the outer planet is less than half as dense, implying that it contains a greater portion of low density components such as H2O and/or H2.

We present new multi-band transit photometry of three small (R$_{p}$ < 6 R$_{\oplus}$), short-period (P < 6 days) Kepler planet candidates acquired with the Gran Telescopio Canarias. These observations supplement the results presented in Col\'on & Ford (2011) and Col\'on et al. (2012), where we used multicolor transit photometry of five Kepler planet candidates to search for wavelength-dependent transit depths and either validate planet candidates or identify eclipsing binary false positives within our sample. In those previous studies, we provided evidence that three targets were false positives and two targets were planets. Here, we present observations that provide evidence supporting a planetary nature for KOI 439.01 and KOI 732.01, and we find that KOI 531.01, a 6 R$_{\oplus}$ planet candidate around an M dwarf, is likely a false positive. We also present a discussion of the purported "sub-Jovian desert" in the orbital period-planet radius plane, which cannot be easily explained by observational bias. Both KOI 439.01 and KOI 732.01 are likely planets located within the so-called desert and should be investigated with further follow-up observations. As only ~30 of the ~3600 currently active Kepler planet candidates are located within the sub-Jovian desert, it will be interesting to see if these candidates also survive the vetting process and fill in the gap in the period-radius plane. Confirming planets in this regime will be important for understanding planetary migration and evolution processes, and we urge additional follow-up observations of these planet candidates to confirm their nature.

Mazeh, Holczer, and Shporer (2015) have presented an approach that can, in principle, use the derived transit timing variation (TTV) of some transiting planets observed by the $Kepler$ mission to distinguish between prograde and retrograde motion of their orbits with respect to their parent stars' rotation. The approach utilizes TTVs induced by spot-crossing events that occur when the planet moves across a spot on the stellar surface, looking for a correlation between the derived TTVs and the stellar brightness derivatives at the corresponding transits. This can work even in data that cannot temporally resolve the spot-crossing events themselves. Here we apply this approach to the $Kepler$ KOIs, identifying nine systems where the photometric spot modulation is large enough and the transit timing accurate enough to allow detection of a TTV-brightness-derivatives correlation. Of those systems five show highly significant prograde motion (Kepler-17b, Kepler-71b, KOI-883.01, KOI-895.01, and KOI-1074.01), while no system displays retrograde motion, consistent with the suggestion that planets orbiting cool stars have prograde motion. All five systems have impact parameter $0.2\lesssim b\lesssim0.5$, and all systems within that impact parameter range show significant correlation, except HAT-P-11b where the lack of a correlation follows its large stellar obliquity. Our search suffers from an observational bias against detection of high impact parameter cases, and the detected sample is extremely small. Nevertheless, our findings may suggest that stellar spots, or at least the larger ones, tend to be located at a low stellar latitude, but not along the stellar equator, similar to the Sun.

(Abridged) NASA's Kepler mission has provided several thousand transiting planet candidates, yet only a small subset have been confirmed as true planets. Therefore, the most fundamental question about these candidates is the fraction of bona fide planets. Estimating the rate of false positives of the overall Kepler sample is necessary to derive the planet occurrence rate. We present the results from two large observational campaigns that were conducted with the Spitzer telescope during the the Kepler mission. These observations are dedicated to estimating the false positive rate (FPR) amongst the Kepler candidates. We select a sub-sample of 51 candidates, spanning wide ranges in stellar, orbital and planetary parameter space, and we observe their transits with Spitzer at 4.5 microns. We use these observations to measures the candidate's transit depths and infrared magnitudes. A bandpass-dependent depth alerts us to the potential presence of a blending star that could be the source of the observed eclipse: a false-positive scenario. For most of the candidates (85%), the transit depths measured with Kepler are consistent with the depths measured with Spitzer as expected for planetary objects, while we find that the most discrepant measurements are due to the presence of unresolved stars that dilute the photometry. The Spitzer constraints on their own yield FPRs between 5-40%, depending on the KOIs. By considering the population of the Kepler field stars, and by combining follow-up observations (imaging) when available, we find that the overall FPR of our sample is low. The measured upper limit on the FPR of our sample is 8.8% at a confidence level of 3 sigma. This observational result, which uses the achromatic property of planetary transit signals that is not investigated by the Kepler observations, provides an independent indication that Kepler's false positive rate is low.

The Kepler mission has yielded a large number of planet candidates from among the Kepler Objects of Interest (KOIs), but spectroscopic follow-up of these relatively faint stars is a serious bottleneck in confirming and characterizing these systems. We present motivation and survey design for an ongoing project with the SDSS-III multiplexed APOGEE near-infrared spectrograph to monitor hundreds of KOI host stars. We report some of our first results using representative targets from our sample, which include current planet candidates that we find to be false positives, as well as candidates listed as false positives that we do not find to be spectroscopic binaries. With this survey, KOI hosts are observed over ~20 epochs at a radial velocity precision of 100-200 m/s. These observations can easily identify a majority of false positives caused by physically-associated stellar or substellar binaries, and in many cases, fully characterize their orbits. We demonstrate that APOGEE is capable of achieving RV precision at the 100-200 m/s level over long time baselines, and that APOGEE's multiplexing capability makes it substantially more efficient at identifying false positives due to binaries than other single-object spectrographs working to confirm KOIs as planets. These APOGEE RVs enable ancillary science projects, such as studies of fundamental stellar astrophysics or intrinsically rare substellar companions. The coadded APOGEE spectra can be used to derive stellar properties (T_eff, log(g)) and chemical abundances of over a dozen elements to probe correlations of planet properties with individual elemental abundances.

An intriguing trend among \kepler's multi-planet systems is an overabundance of planet pairs with period ratios just wide of a mean motion resonance (MMR) and a dearth of systems just narrow of them. Traditional planet formation models are at odds with these observations. They are also in contrast with the period ratios of radial-velocity-discovered multi-planet systems which tend to pile up at 2:1 MMR. We propose that gas-disk migration traps planets in a MMR. After gas dispersal, orbits of these trapped planets are altered through interaction with a residual planetesimal disk. We study the effects of planetesimal disk interactions on planet pairs trapped in 2:1 MMR using planets of mass typical of the Kepler planet candidates (KPC) and explore large ranges for the mass, and density profile of the planetesimal disk. We find that planet-planetesimal disk interactions naturally create the observed asymmetry in period-ratio distribution for large ranges of planetesimal disk and planet properties. If the planetesimal disk mass is above a threshold of ~0.2x the planet mass, these interactions typically disrupt MMR. Afterwards, the planets migrate in such a way that the final period-ratio is slightly higher than the integer ratio corresponding to the initial MMR. Below this threshold these interactions typically cannot disrupt the resonance and the period ratio stays close to the integer ratio. The threshold explains why the more massive planet pairs found by RV surveys are still in resonance. We encourage future research to explore how significantly the associated accretion would change the planets' atmospheric and surface properties.

We present an update to seven stars with long-period planets or planetary candidates using new and archival radial velocities from Keck-HIRES and literature velocities from other telescopes. Our updated analysis better constrains orbital parameters for these planets, four of which are known multi-planet systems. HD 24040 b and HD 183263 c are super-Jupiters with circular orbits and periods longer than 8 yr. We present a previously unseen linear trend in the residuals of HD 66428 indicative on an additional planetary companion. We confirm that GJ 849 is a multi-planet system and find a good orbital solution for the c component: it is a $1 M_{\rm Jup}$ planet in a 15 yr orbit (the longest known for a planet orbiting an M dwarf). We update the HD 74156 double-planet system. We also announce the detection of HD 145934 b, a $2 M_{\rm Jup}$ planet in a 7.5 yr orbit around a giant star. Two of our stars, HD 187123 and HD 217107, at present host the only known examples of systems comprising a hot Jupiter and a planet with a well constrained period $> 5$ yr, and with no evidence of giant planets in between. Our enlargement and improvement of long-period planet parameters will aid future analysis of origins, diversity, and evolution of planetary systems.

We present a solution to the long outstanding meter barrier problem in planet formation theory. As solids spiral inward due to aerodynamic drag, they will enter disk regions that are characterized by high temperatures, densities, and pressures. High partial pressures of rock vapor can suppress solid evaporation, and promote collisions between partially molten solids, allowing rapid growth. This process should be ubiquitous in planet-forming disks, which may be evidenced by the abundant class of Systems with Tightly-packed Inner Planets (STIPs) discovered by the NASA Kepler mission.

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.

We present an updated study of the planets known to orbit 55 Cancri A using 1,418 high-precision radial velocity observations from four observatories (Lick, Keck, Hobby-Eberly Telescope, Harlan J. Smith Telescope) and transit time/durations for the inner-most planet, 55 Cancri "e" (Winn et al. 2011). We provide the first posterior sample for the masses and orbital parameters based on self-consistent n-body orbital solutions for the 55 Cancri planets, all of which are dynamically stable (for at least $10^8$ years). We apply a GPU version of Radial velocity Using N-body Differential evolution Markov Chain Monte Carlo (RUN DMC; B. Nelson et al. 2014) to perform a Bayesian analysis of the radial velocity and transit observations. Each of the planets in this remarkable system has unique characteristics. Our investigation of high-cadence radial velocities and priors based on space-based photometry yields an updated mass estimate for planet "e" ($8.09\pm0.26$ M$_\oplus$), which affects its density ($5.51\pm^{1.32}_{1.00}$ g cm$^{-3}$) and inferred bulk composition. Dynamical stability dictates that the orbital plane of planet "e" must be aligned to within $60^o$ of the orbital plane of the outer planets (which we assume to be coplanar). The mutual interactions between the planets "b" and "c" may develop an apsidal lock about $180^o$. We find 36-45% of all our model systems librate about the anti-aligned configuration with an amplitude of $51^o\pm^{6^o}_{10^o}$. Other cases showed short-term perturbations in the libration of $\varpi_b-\varpi_c$, circulation, and nodding, but we find the planets are not in a 3:1 mean-motion resonance. A revised orbital period and eccentricity for planet "d" pushes it further toward the closest known Jupiter analog in the exoplanet population.

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

Doppler planet searches revealed that many giant planets orbit close to their host star or in highly eccentric orbits. These and subsequent observations inspired new theories of planet formation that invoke gravitation interactions in multiple planet systems to explain the excitation of orbital eccentricities and even short-period giant planets. Recently, NASA's Kepler mission has identified over 300 systems with multiple transiting planet candidates, including many potentially rocky planets. Most of these systems include multiple planets with closely-spaced orbits and sizes between that of Earth and Neptune. These systems represent yet another new and unexpected class of planetary systems and provide an opportunity to test the theories developed to explain the properties of giant exoplanets. Presently, we have limited knowledge about such planetary systems, mostly about their sizes and orbital periods. With the advent of long-term, nearly continuous monitoring by Kepler, the method of transit timing variations (TTVs) has blossomed as a new technique for characterizing the gravitational effects of mutual planetary perturbations for hundreds of planets. TTVs can provide precise (but complex) constraints on planetary masses, densities and orbits, even for planetary systems with faint host stars. In the coming years, astronomers will translate TTV observations into increasingly powerful constraints on the formation and orbital evolution of planetary systems with low-mass planets. Between TTVs, improved Doppler surveys, high-contrast imaging campaigns and microlensing surveys, astronomers can look forward to a much better understanding of planet formation in the coming decade.

The Kepler mission has discovered over 2500 exoplanet candidates in the first two years of spacecraft data, with approximately 40% of them in candidate multi-planet systems. The high rate of multiplicity combined with the low rate of identified false-positives indicates that the multiplanet systems contain very few false-positive signals due to other systems not gravitationally bound to the target star (Lissauer, J. J., et al., 2012, ApJ 750, 131). False positives in the multi- planet systems are identified and removed, leaving behind a residual population of candidate multi-planet transiting systems expected to have a false-positive rate less than 1%. We present a sample of 340 planetary systems that contain 851 planets that are validated to substantially better than the 99% confidence level; the vast majority of these have not been previously verified as planets. We expect ~2 unidentified false-positives making our sample of planet very reliable. We present fundamental planetary properties of our sample based on a comprehensive analysis of Kepler light curves and ground-based spectroscopy and high-resolution imaging. Since we do not require spectroscopy or high-resolution imaging for validation, some of our derived parameters for a planetary system may be systematically incorrect due to dilution from light due to additional stars in the photometric aperture. None the less, our result nearly doubles the number of verified exoplanets.

We extend the statistical analysis of Lissauer et al. (2012, ApJ 750, 112), which demonstrates that the overwhelming majority of Kepler candidate multiple transiting systems (multis) represent true transiting planets, and develop therefrom a procedure to validate large numbers of planet candidates in multis as bona fide exoplanets. We show that this statistical framework correctly estimates the abundance of false positives already identified around Kepler targets with multiple sets of transit-like signatures based on their abundance around targets with single sets of transit-like signatures. We estimate the number of multis that represent split systems of one or more planets orbiting each component of a binary star system. We use the high reliability rate for multis to validate more than one dozen particularly interesting multi-planet systems are validated in a companion paper by Rowe et al. (2014, ApJ, this issue). We note that few very short period (P < 1.6 days) planets orbit within multiple transiting planet systems and discuss possible reasons for their absence. There also appears to be a shortage of planets with periods exceeding a few months in multis.

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

We present narrow-band photometric measurements of the exoplanet GJ 1214b using the 10.4 m Gran Telescopio Canarias (GTC) and the OSIRIS instrument. Using tuneable filters we observed a total of five transits, three of which were observed at two wavelengths nearly simultaneously, producing a total of eight individual light curves, six of these probed the possible existence of a methane absorption feature in the 8770 - 8850 {\AA} region at high resolution. We detect no increase in the planet-to-star radius ratio across the methane feature with a change in radius ratio of $\Delta$R = -0.0007 $\pm$ 0.0017 corresponding to a scale height (H) change of -0.5 $\pm$ 1.2 H across the methane feature, assuming a hydrogen dominated atmosphere. We find a variety of water and cloudy atmospheric models fit the data well, but find that cloud-free models provide poor fits. These observations support a flat transmission spectrum resulting from the presence of a high-altitude haze or a water-rich atmosphere, in agreement with previous studies. In this study the observations are predominantly limited by the photometric quality and the limited number of data points (resulting from a long observing cadence), which make the determination of the systematic noise challenging. With tuneable filters capable of high resolution measurements (R ~ 600 - 750) of narrow absorption features, the interpretation of our results are also limited by the absence of high resolution methane models below 1 $\mu$m.

This chapter concerns the long-term dynamical evolution of planetary systems from both theoretical and observational perspectives. We begin by discussing the planet-planet interactions that take place within our own Solar System. We then describe such interactions in more tightly-packed planetary systems. As planet-planet interactions build up, some systems become dynamically unstable, leading to strong encounters and ultimately either ejections or collisions of planets. After discussing the basic physical processes involved, we consider how these interactions apply to extrasolar planetary systems and explore the constraints provided by observed systems. The presence of a residual planetesimal disc can lead to planetary migration and hence cause instabilities induced by resonance crossing; however, such discs can also stabilise planetary systems. The crowded birth environment of a planetary system can have a significant impact: close encounters and binary companions can act to destabilise systems, or sculpt their properties. In the case of binaries, the Kozai mechanism can place planets on extremely eccentric orbits which may later circularise to produce hot Jupiters.

In the 20+ years of Doppler observations of stars, scientists have uncovered a diverse population of extrasolar multi-planet systems. A common technique for characterizing the orbital elements of these planets is Markov chain Monte Carlo (MCMC), using a Keplerian model with random walk proposals and paired with the Metropolis-Hastings algorithm. For approximately a couple of dozen planetary systems with Doppler observations, there are strong planet-planet interactions due to the system being in or near a mean-motion resonance (MMR). An N-body model is often required to accurately describe these systems. Further computational difficulties arise from exploring a high-dimensional parameter space ($\sim$7 x number of planets) that can have complex parameter correlations. To surmount these challenges, we introduce a differential evolution MCMC (DEMCMC) applied to radial velocity data while incorporating self-consistent N-body integrations. Our Radial velocity Using N-body DEMCMC (RUN DMC) algorithm improves upon the random walk proposal distribution of the traditional MCMC by using an ensemble of Markov chains to adaptively improve the proposal distribution. We describe the methodology behind the algorithm, along with results of tests for accuracy and performance. We find that most algorithm parameters have a modest effect on the rate of convergence. However, the size of the ensemble can have a strong effect on performance. We show that the optimal choice depends on the number of planets in a system, as well as the computer architecture used and the resulting extent of parallelization. While the exact choices of optimal algorithm parameters will inevitably vary due to the details of individual planetary systems, we offer recommendations for choosing the DEMCMC algorithm's algorithmic parameters that result in excellent performance for a wide variety of planetary systems.

Stars hosting hot Jupiters are often observed to have high obliquities, whereas stars with multiple co-planar planets have been seen to have low obliquities. This has been interpreted as evidence that hot-Jupiter formation is linked to dynamical disruption, as opposed to planet migration through a protoplanetary disk. We used asteroseismology to measure a large obliquity for Kepler-56, a red giant star hosting two transiting co-planar planets. These observations show that spin-orbit misalignments are not confined to hot-Jupiter systems. Misalignments in a broader class of systems had been predicted as a consequence of torques from wide-orbiting companions, and indeed radial-velocity measurements revealed a third companion in a wide orbit in the Kepler-56 system.

We have analyzed new and previously published radial velocity observations of MARVELS-1, known to have an ostensibly substellar companion in a ~6- day orbit. We find significant (~100 m/s) residuals to the best-fit model for the companion, and these residuals are naively consistent with an interior giant planet with a P = 1.965d in a nearly perfect 3:1 period commensuribility (|Pb/Pc - 3| < 10^{-4}). We have performed several tests for the reality of such a companion, including a dynamical analysis, a search for photometric variability, and a hunt for contaminating stellar spectra. We find many reasons to be critical of a planetary interpretation, including the fact that most of the three-body dynamical solutions are unstable. We find no evidence for transits, and no evidence of stellar photometric variability. We have discovered two apparent companions to MARVELS-1 with adaptive optics imaging at Keck; both are M dwarfs, one is likely bound, and the other is likely a foreground object. We explore false-alarm scenarios inspired by various curiosities in the data. Ultimately, a line profile and bisector analysis lead us to conclude that the ~100 m/s residuals are an artifact of spectral contamination from a stellar companion contributing ~15-30% of the optical light in the system. We conclude that origin of this contamination is the previously detected radial velocity companion to MARVELS-1, which is not, as previously reported, a brown dwarf, but in fact a G dwarf in a face-on orbit.

We present an updated analysis of radial velocity data of the HD 82943 planetary system based on 10 years of measurements obtained with the Keck telescope. Previous studies have shown that the HD 82943 system has two planets that are likely in 2:1 mean-motion resonance (MMR), with the orbital periods about 220 and 440 days (Lee et al. 2006). However, alternative fits that are qualitatively different have also been suggested, with two planets in a 1:1 resonance (Gozdziewski & Konacki 2006) or three planets in a Laplace 4:2:1 resonance (Beauge et al. 2008). Here we use \c{hi}2 minimization combined with parameter grid search to investigate the orbital parameters and dynamical states of the qualitatively different types of fits, and we compare the results to those obtained with the differential evolution Markov chain Monte Carlo method. Our results support the coplanar 2:1 MMR configuration for the HD 82943 system, and show no evidence for either the 1:1 or 3-planet Laplace resonance fits. The inclination of the system with respect to the sky plane is well constrained at about 20(+4.9 -5.5) degree, and the system contains two planets with masses of about 4.78 MJ and 4.80 MJ (where MJ is the mass of Jupiter) and orbital periods of about 219 and 442 days for the inner and outer planet, respectively. The best fit is dynamically stable with both eccentricity-type resonant angles {\theta}1 and {\theta}2 librating around 0 degree.

We propose a scientific program to complete a census of planets, characterizing their masses, orbital properties, and dynamical histories using continued observations of the Kepler field of view with the Kepler spacecraft in a two reaction wheel mission. Even with a significantly reduced photometric precision, extending time-domain observations of this field is uniquely capable of pursuing several critical science goals: 1) measuring the architectures of planetary systems by identifying non-transiting planets interleaved among known transiting planets, 2) establishing the mass-radius relationship for planets in the important transition region between small, gas-rich sub-Neptune planets and large, rocky super-Earths, and 3) uncovering dynamical evidence of the formation and evolution of the inner regions of planetary systems. To meet these objectives, the unique multi-object observing capabilities of Kepler will be used in a set of concurrent campaigns with specific motivations. These campaigns focus largely on the ability to interpret Transit Timing Variations (TTVs) that result from dynamical interactions among planets in a system and include: 1) observations of systems that exhibit large TTVs and are particularly rich in dynamical information, 2) observations of systems where additional transit times will yield mass measurements of the constituent planets, 3) observations of systems where the TTV signal evolves over very long timescales, and 4) observations of systems with long-period planet candidates where additional transits will remove orbital period ambiguities caused by gaps in the original Kepler data.

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

Following Ford et al. (2011, 2012) and Steffen et al. (2012) we derived the transit timing of 1960 Kepler KOIs using the pre-search data conditioning (PDC) light curves of the first twelve quarters of the Kepler data. For 721 KOIs with large enough SNRs, we obtained also the duration and depth of each transit. The results are presented as a catalog for the community to use. We derived a few statistics of our results that could be used to indicate significant variations. Including systems found by previous works, we have found 130 KOIs that showed highly significant TTVs, and 13 that had short-period TTV modulations with small amplitudes. We consider two effects that could cause apparent periodic TTV - the finite sampling of the observations and the interference with the stellar activity, stellar spots in particular. We briefly discuss some statistical aspects of our detected TTVs. We show that the TTV period is correlated with the orbital period of the planet and with the TTV amplitude.

The NASA Kepler mission has revealed an abundant class of Systems with Tightly-packed Inner Planets (STIPs). The current paradigm for planet formation suggests that small planetesimals will quickly spiral into the host star due to aerodynamic drag, preventing rocky planet formation. In contrast, we find that aerodynamic drift, when acting on an ensemble of solids, can concentrate mass at short orbital periods in gaseous disks. Sublimation fronts may further aid this process. Kepler data suggest that the innermost known planets are found near the silicate sublimation zone. STIP planets should have a wide range of volatile fractions due to aerodynamic drift and H2 dissociation-driven gas accretion. We further propose that the low mass of Mars is evidence that the Solar System was once a proto-STIP.

Since the discovery of the first exoplanet we have known that other planetary systems can look quite unlike our own. However, until recently we have only been able to probe the upper range of the planet size distribution. The high precision of the Kepler space telescope has allowed us to detect planets that are the size of Earth and somewhat smaller, but no previous planets have been found that are smaller than those we see in our own Solar System. Here we report the discovery of a planet significantly smaller than Mercury. This tiny planet is the innermost of three planets that orbit the Sun-like host star, which we have designated Kepler-37. Owing to its extremely small size, similar to that of Earth's Moon, and highly irradiated surface, Kepler-37b is probably a rocky planet with no atmosphere or water, similar to Mercury.

We present the detection of five planets -- Kepler-62b, c, d, e, and f -- of size 1.31, 0.54, 1.95, 1.61 and 1.41 Earth radii, orbiting a K2V star at periods of 5.7, 12.4, 18.2, 122.4 and 267.3 days, respectively. The outermost planets (Kepler-62e & -62f) are super-Earth-size (1.25 < planet radius/earth radius < 2.0) planets in the habitable zone (HZ) of their host star, receiving 1.2 +- 0.2 and 0.41 +- 0.05 times the solar flux at Earth's orbit. Theoretical models of Kepler-62e and -62f for a stellar age of ~7 Gyr suggest that both planets could be solid: either with a rocky composition or composed of mostly solid water in their bulk.

The Kepler Mission was designed to measure the frequency of Earth-size planets in the habitable zone of Sun-like stars. A crucial component for recovering the underlying planet population from a sample of detected planets is understanding the completeness of that sample - what fraction of the planets that could have been discovered in a given data set were actually detected. Here we outline the information required to determine the sample completeness, and describe an experiment to address a specific aspect of that question, which is the issue of transit signal recovery. We investigate the extent to which the Kepler pipeline preserves individual transit signals by injecting simulated transits into the pixel-level data, processing the modified pixels through the pipeline, and comparing the measured transit signal-to-noise ratio (SNR) to that expected without perturbation by the pipeline. We inject simulated transit signals across the full focal plane for a set of observations of length 89 days. On average, we find that the SNR of the injected signal is recovered at MS = 0.9973(+/-0.0012)xBS-0.0151(+/-0.0049), where MS is the measured SNR and BS is the baseline, or expected, SNR. The 1{\sigma} width of the distribution around this correlation is +/-2.64%. We discuss the pipeline processes that cause the measured SNR to deviate significantly from the baseline SNR; these are primarily the handling of data adjacent to spacecraft re-pointings and the removal of harmonics prior to the measurement of the SNR. Finally we outline the further work required to characterise the completeness of the Kepler pipeline.

NASA's Kepler Mission has revealed two transiting planets orbiting Kepler-68. Follow-up Doppler measurements have established the mass of the innermost planet and revealed a third jovian-mass planet orbiting beyond the two transiting planets. Kepler-68b, in a 5.4 day orbit has mass 8.3 +/- 2.3 Earth, radius 2.31 +/- 0.07 Earth radii, and a density of 3.32 +/- 0.92 (cgs), giving Kepler-68b a density intermediate between that of the ice giants and Earth. Kepler-68c is Earth-sized with a radius of 0.953 Earth and transits on a 9.6 day orbit; validation of Kepler-68c posed unique challenges. Kepler-68d has an orbital period of 580 +/- 15 days and minimum mass of Msin(i) = 0.947 Jupiter. Power spectra of the Kepler photometry at 1-minute cadence exhibit a rich and strong set of asteroseismic pulsation modes enabling detailed analysis of the stellar interior. Spectroscopy of the star coupled with asteroseismic modeling of the multiple pulsation modes yield precise measurements of stellar properties, notably Teff = 5793 +/- 74 K, M = 1.079 +/- 0.051 Msun, R = 1.243 +/- 0.019 Rsun, and density 0.7903 +/- 0.0054 (cgs), all measured with fractional uncertainties of only a few percent. Models of Kepler-68b suggest it is likely composed of rock and water, or has a H and He envelope to yield its density of about 3 (cgs).

Planets embedded within dust disks may drive the formation of large scale clumpy dust structures by trapping dust into resonant orbits. Detection and subsequent modeling of the dust structures would help constrain the mass and orbit of the planet and the disk architecture, give clues to the history of the planetary system, and provide a statistical estimate of disk asymmetry for future exoEarth-imaging missions. Here we present the first search for these resonant structures in the inner regions of planetary systems by analyzing the light curves of hot Jupiter planetary candidates identified by the Kepler mission. We detect only one candidate disk structure associated with KOI 838.01 at the 3-sigma confidence level, but subsequent radial velocity measurements reveal that KOI 838.01 is a grazing eclipsing binary and the candidate disk structure is a false positive. Using our null result, we place an upper limit on the frequency of dense exozodi structures created by hot Jupiters. We find that at the 90% confidence level, less than 21% of Kepler hot Jupiters create resonant dust clumps that lead and trail the planet by ~90 degrees with optical depths >~5*10^-6, which corresponds to the resonant structure expected for a lone hot Jupiter perturbing a dynamically cold dust disk 50 times as dense as the zodiacal cloud.

We present a study of the relative sizes of planets within the multiple candidate systems discovered with the $Kepler$ mission. We have compared the size of each planet to the size of every other planet within a given planetary system after correcting the sample for detection and geometric biases. We find that for planet-pairs for which one or both objects is approximately Neptune-sized or larger, the larger planet is most often the planet with the longer period. No such size--location correlation is seen for pairs of planets when both planets are smaller than Neptune. Specifically, if at least one planet in a planet-pair has a radius of $\gtrsim 3R_\oplus$, $68\pm 6%$ of the planet pairs have the inner planet smaller than the outer planet, while no preferred sequential ordering of the planets is observed if both planets in a pair are smaller than $\lesssim3 R_\oplus$.

We report the radial-velocity discovery of a second planetary mass companion to the K0 V star HD 37605, which was already known to host an eccentric, P~55 days Jovian planet, HD 37605b. This second planet, HD 37605c, has a period of ~7.5 years with a low eccentricity and an Msini of ~3.4 MJup. Our discovery was made with the nearly 8 years of radial velocity follow-up at the Hobby-Eberly Telescope and Keck Observatory, including observations made as part of the Transit Ephemeris Refinement and Monitoring Survey (TERMS) effort to provide precise ephemerides to long-period planets for transit follow-up. With a total of 137 radial velocity observations covering almost eight years, we provide a good orbital solution of the HD 37605 system, and a precise transit ephemeris for HD 37605b. Our dynamic analysis reveals very minimal planet-planet interaction and an insignificant transit time variation. Using the predicted ephemeris, we performed a transit search for HD 37605b with the photometric data taken by the T12 0.8-m Automatic Photoelectric Telescope (APT) and the Microvariability and Oscillations of Stars (MOST) satellite. Though the APT photometry did not capture the transit window, it characterized the stellar activity of HD 37605, which is consistent of it being an old, inactive star, with a tentative rotation period of 57.67 days. The MOST photometry enabled us to report a dispositive null detection of a non-grazing transit for this planet. Within the predicted transit window, we exclude an edge-on predicted depth of 1.9% at >>10sigma, and exclude any transit with an impact parameter b>0.951 at greater than 5sigma. We present the BOOTTRAN package for calculating Keplerian orbital parameter uncertainties via bootstrapping. We found consistency between our orbital parameters calculated by the RVLIN package and error bars by BOOTTRAN with those produced by a Bayesian analysis using MCMC.

We present Swarm-NG, a C++ library for the efficient direct integration of many n-body systems using highly-parallel Graphics Processing Unit (GPU), such as NVIDIA's Tesla T10 and M2070 GPUs. While previous studies have demonstrated the benefit of GPUs for n-body simulations with thousands to millions of bodies, Swarm-NG focuses on many few-body systems, e.g., thousands of systems with 3...15 bodies each, as is typical for the study of planetary systems. Swarm-NG parallelizes the simulation, including both the numerical integration of the equations of motion and the evaluation of forces using NVIDIA's "Compute Unified Device Architecture" (CUDA) on the GPU. Swarm-NG includes optimized implementations of 4th order time-symmetrized Hermite integration and mixed variable symplectic integration, as well as several sample codes for other algorithms to illustrate how non-CUDA-savvy users may themselves introduce customized integrators into the Swarm-NG framework. To optimize performance, we analyze the effect of GPU-specific parameters on performance under double precision. Applications of Swarm-NG include studying the late stages of planet formation, testing the stability of planetary systems and evaluating the goodness-of-fit between many planetary system models and observations of extrasolar planet host stars (e.g., radial velocity, astrometry, transit timing). While Swarm-NG focuses on the parallel integration of many planetary systems,the underlying integrators could be applied to a wide variety of problems that require repeatedly integrating a set of ordinary differential equations many times using different initial conditions and/or parameter values.

While hundreds of planets have been discovered around field stars, only a few are known in star clusters. To explain the lack of short-period giant planets in globular clusters (GC), such as 47 Tucane and \omega\ Centauri, it has been suggested that their low metallicities may have prevented planet formation. Alternatively, the high rates of close stellar encounters in these clusters may have influenced the formation and subsequent evolution of planetary systems. How common are planets in clusters around normal main-sequence stars? Here we consider whether this question can be addressed using data from the Kepler mission. The Kepler field of view contains 4 low-density (relative to GCs) open clusters where the metallicities are about solar (or even higher) and stellar encounters are much less frequent than in typical GCs. We provide detailed $N$-body models and show that most planets in Kepler-detectable orbits are not significantly perturbed by stellar encounters in these open clusters. We focus on the most massive cluster, NGC 6791, which has super-solar metallicity, and find that if planets formed in this cluster at the same frequency as observed in the field, Kepler could detect 1 -- 20 transiting planets depending on the planet-size distribution and the duration of data collection. Due to the large distance to NGC 6791 Kepler will have to search relatively faint ($K_p<20$) stars for the full extended mission to achieve such a yield.

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

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

We confirm 27 planets in 13 planetary systems by showing the existence of statistically significant anti-correlated transit timing variations (TTVs), which demonstrates that the planet candidates are in the same system, and long-term dynamical stability, which places limits on the masses of the candidates---showing that they are planetary. %This overall method of planet confirmation was first applied to \kepler systems 23 through 32. All of these newly confirmed planetary systems have orbital periods that place them near first-order mean motion resonances (MMRs), including 6 systems near the 2:1 MMR, 5 near 3:2, and one each near 4:3, 5:4, and 6:5. In addition, several unconfirmed planet candidates exist in some systems (that cannot be confirmed with this method at this time). A few of these candidates would also be near first order MMRs with either the confirmed planets or with other candidates. One system of particular interest, Kepler-56 (KOI-1241), is a pair of planets orbiting a 12th magnitude, giant star with radius over three times that of the Sun and effective temperature of 4900 K---among the largest stars known to host a transiting exoplanetary system.

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.

At least two multi-planetary systems in a 4:3 mean motion resonance have been found by radial velocity surveys. These planets are gas giants and the systems are only stable when protected by a resonance. Additionally the Kepler mission has detected at least 4 strong candidate planetary systems with a period ratio close to 4:3. This paper investigates traditional dynamical scenarios for the formation of these systems. We systematically study migration scenarios with both N-body and hydro-dynamic simulations. We investigate scenarios involving the in-situ formation of two planets in resonance. We look at the results from finely tuned planet-planet scattering simulations with gas disk damping. Finally, we investigate a formation scenario involving isolation-mass embryos. Although the combined planet-planet scattering and damping scenario seems promising, none of the above scenarios is successful in forming enough systems in 4:3 resonance with planetary masses similar to the observed ones. This is a negative result but it has important implications for planet formation. Previous studies were successful in forming 2:1 and 3:2 resonances. This is generally believed to be evidence of planet migration. We highlight the main differences between those studies and our failure in forming a 4:3 resonance. We also speculate on more exotic and complicated ideas. These results will guide future investigators toward exploring the above scenarios and alternative mechanisms in a more general framework.

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.

Using the OSIRIS instrument installed on the 10.4-m Gran Telescopio Canarias (GTC) we acquired multi-color transit photometry of four small (Rp < 5 R_Earth) short-period (P < 6 days) planet candidates recently identified by the Kepler space mission. These observations are part of a program to constrain the false positive rate for small, short-period Kepler planet candidates. Since planetary transits should be largely achromatic when observed at different wavelengths (excluding the small color changes due to stellar limb darkening), we use the observed transit color to identify candidates as either false positives (e.g., a blend with a stellar eclipsing binary either in the background/foreground or bound to the target star) or validated planets. Our results include the identification of KOI 225.01 and KOI 1187.01 as false positives and the tentative validation of KOI 420.01 and KOI 526.01 as planets. The probability of identifying two false positives out of a sample of four targets is less than 1%, assuming an overall false positive rate for Kepler planet candidates of 10% (as estimated by Morton & Johnson 2011). Therefore, these results suggest a higher false positive rate for the small, short-period Kepler planet candidates than has been theoretically predicted by other studies which consider the Kepler planet candidate sample as a whole. Furthermore, our results are consistent with a recent Doppler study of short-period giant Kepler planet candidates (Santerne et al. 2012). We also investigate how the false positive rate for our sample varies with different planetary and stellar properties. Our results suggest that the false positive rate varies significantly with orbital period and is largest at the shortest orbital periods (P < 3 days), where there is a corresponding rise in the number of detached eclipsing binary stars... (truncated)

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

The chance that a planetary system will interact with another member of its host star's nascent cluster would be greatly increased if gas giant planets form in situ on wide orbits. In this paper, we explore the outcomes of planet-planet scattering for a distribution of multiplanet systems that all have one of the planets on an initial orbit of 100 AU. The scattering experiments are run with and without stellar flybys. We convolve the outcomes with distributions for protoplanetary disk and stellar cluster sizes to generalize the results where possible. We find that the frequencies of large mutual inclinations and high eccentricities are sensitive to the number of planets in a system, but not strongly to stellar flybys. However, flybys do play a role in changing the low and moderate portions of the mutual inclination distributions, and erase dynamically cold initial conditions on average. Wide-orbit planets can be mixed throughout the planetary system, and in some cases, can potentially become hot Jupiters, which we demonstrate using scattering experiments that include a tidal damping model. If planets form on wide orbits in situ, then there will be discernible differences in the proper motion distributions of a sample of wide-orbit planets compared with a pure scattering formation mechanism. Stellar flybys can enhance the frequency of ejections in planetary systems, but auto-ionization is likely to remain the dominant source of free-floating planets.

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.

The dynamical evolution of planetary systems leaves observable signatures in debris disks. Optical images trace micron-sized grains, which are strongly affected by stellar radiation and need not coincide with their parent body population. Observations of mm-size grains accurately trace parent bodies, but previous images lack the resolution and sensitivity needed to characterize the ring's morphology. Here we present ALMA 350 GHz observations of the Fomalhaut debris ring. These observations demonstrate that the parent body population is 13-19 AU wide with a sharp inner and outer boundary. We discuss three possible origins for the ring, and suggest that debris confined by shepherd planets is the most consistent with the ring's morphology.

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

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

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.

The Kepler mission has discovered a plethora of multiple transiting planet candidate exosystems, many of which feature putative pairs of planets near mean motion resonance commensurabilities. Identifying potentially resonant systems could help guide future observations and enhance our understanding of planetary formation scenarios. We develop and apply an algebraic method to determine which Kepler 2-planet systems cannot be in a 1st-4th order resonance, given the current, publicly available data. This method identifies when any potentially resonant angle of a system must circulate. We identify and list 70 near-resonant systems which cannot actually reside in resonance, assuming a widely-used formulation for deriving planetary masses from their observed radii and that these systems do not contain unseen bodies that affect the interactions of the observed planets. This work strengthens the argument that a high fraction of exoplanetary systems may be near resonance but not actually in resonance.

This paper reports the discovery and characterization of the transiting hot giant exoplanet Kepler-17b. The planet has an orbital period of 1.486 days, and radial velocity measurements from the Hobby-Eberly Telescope (HET) show a Doppler signal of 420+/-15 m.s-1. From a transit-based estimate of the host star's mean density, combined with an estimate of the stellar effective temperature T_eff=5630+/-100 K from high-resolution spectra, we infer a stellar host mass of 1.061+/-0.067 M_sun and a stellar radius of 1.019+/-0.033 R_jup. We estimate the planet mass and radius to be Mp=2.450+/-0.114 M_jup and Rp=1.312+/-0.018 R_jup and a planet density near 1.35 g.cm-3. The host star is active, with dark spots that are frequently occulted by the planet. The continuous monitoring of the star reveals a stellar rotation period of 11.89 days, 8 times the the planet's orbital period; this period ratio produces stroboscopic effects on the occulted starspots. The temporal pattern of these spot-crossing events shows that the planet's orbit is prograde and the star's obliquity is smaller than 15 deg. We detected planetary occultations of Kepler-17b with both the Kepler and Spitzer Space Telescopes. We use these observations to constrain the eccentricity, e, and find that it is consistent with a circular orbit (e<0.0011). The brightness temperatures of the planet the infrared bandpasses are T_3.6um=1880+/-100 K and T4.5um=1770+/-150 K. We measure the optical geometric albedo A_g in the Kepler bandpass and find A_g = 0.10+/-0.02. The observations are best described by atmospheric models for which most of the incident energy is re-radiated away from the day side.

We report the discovery of an eclipsing binary star (KIC 7025851) near KOI 565 (KIC 7025846) based on photometric observations of KOI 565 and several nearby stars acquired in two narrow bandpasses (790.2 and 794.3\pm2.0 nm) nearly-simultaneously with the GTC/OSIRIS. We use the individual photometry in each bandpass as well as the colors of KOI 565 and other nearby stars to determine that the source of the transit signal initially detected by Kepler is not due to a super-Earth-size planet around KOI 565. Instead, we find the source to be a background eclipsing binary star located \sim15 arcsec to the North of KOI 565. We discuss future prospects for using high-precision multi-color photometry from the GTC to determine whether additional Kepler planet candidates have a planetary nature or are instead false positives (e.g., foreground or background eclipsing binaries or hierarchical triple systems). Our approach is complementary to multi-color follow-up observations of Kepler planet candidates currently being conducted with the Spitzer space telescope in the infrared as well as to other follow-up techniques.

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 observations of HD 80606 using the 10.4-m Gran Telescopio Canarias (GTC) and the OSIRIS tunable filter imager. We acquired very-high-precision, narrow-band photometry in four bandpasses around the K I absorption feature during the January 2010 transit of HD 80606b and during out-of-transit observations conducted in January and April of 2010. We obtained differential photometric precisions of \sim 2.08e-4 for the in-transit flux ratio measured at 769.91-nm, which probes the K I line core. We find no significant difference in the in-transit flux ratio between observations at 768.76 and 769.91 nm. Yet, we find a difference of \sim 8.09 \pm 2.88e-4 between these observations and observations at a longer wavelength that probes the K I wing (777.36 nm). While the presence of red noise in the transit data has a non-negligible effect on the uncertainties in the flux ratio, the 777.36-769.91 nm colour during transit shows no effects from red noise and also indicates a significant colour change, with a mean value of \sim 8.99\pm0.62e-4. This large change in the colour is equivalent to a \sim 4.2% change in the apparent planetary radius with wavelength, which is much larger than the atmospheric scale height. This implies the observations probed the atmosphere at very low pressures as well as a dramatic change in the pressure at which the slant optical depth reaches unity between \sim770 and 777 nm. We hypothesize that the excess absorption may be due to K I in a high-speed wind being driven from the exoplanet's exosphere. We discuss the viability of this and alternative interpretations, including stellar limb darkening, starspots, and effects from Earth's atmosphere. We strongly encourage follow-up observations of HD 80606b to confirm the signal measured here. Finally, we discuss the future prospects for exoplanet characterization using tunable filter spectrophotometry.

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 report the discovery of planet Kepler-12b (KOI-20), which at 1.695\pm0.030 RJ is among the handful of planets with super-inflated radii above 1.65 RJ. Orbiting its slightly evolved G0 host with a 4.438-day period, this 0.431\pm0.041 MJ planet is the least-irradiated within this largest-planet-radius group, which has important implications for planetary physics. The planet's inflated radius and low mass lead to a very low density of 0.111\pm0.010 g cm-3. We detect the occultation of the planet at a significance of 3.7{\sigma} in the Kepler bandpass. This yields a geometric albedo of 0.14\pm0.04; the planetary flux is due to a combination of scattered light and emitted thermal flux. We use multiple observations with Warm Spitzer to detect the occultation at 7{\sigma} and 4{\sigma} in the 3.6 and 4.5 {\mu}m bandpasses, respectively. The occultation photometry timing is consistent with a circular orbit, at e < 0.01 (1{\sigma}), and e < 0.09 (3{\sigma}). The occultation detections across the three bands favor an atmospheric model with no dayside temperature inversion. The Kepler occultation detection provides significant leverage, but conclusions regarding temperature structure are preliminary, given our ignorance of opacity sources at optical wavelengths in hot Jupiter atmospheres. If Kepler-12b and HD 209458b, which intercept similar incident stellar fluxes, have the same heavy element masses, the interior energy source needed to explain the large radius of Kepler-12b is three times larger than that of HD 209458b. This may suggest that more than one radius-inflation mechanism is at work for Kepler-12b, or that it is less heavy-element rich than other transiting planets.

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

We apply standard Markov chain Monte Carlo (MCMC) analysis techniques for 50 short- period, single-planet systems discovered with radial velocity technique. We develop a new method for accessing the significance of a non-zero orbital eccentricity, namely {\Gamma} analysis, which combines frequentist bootstrap approach with Bayesian analysis of each simulated data set. We find the eccentricity estimations from {\Gamma} analysis are generally consistent with results from both standard MCMC analysis and previous references. The {\Gamma} method is particular useful for assessing the significance of small eccentricities. Our results suggest that the current sample size is insufficient to draw robust conclusions about the roles of tidal interaction and perturbations in shaping the eccentricity distribution of short-period single-planet systems. We use a Bayesian population analysis to show that a mixture of analytical distributions is a good approximation of the underlying eccentricity distribution. For short-period planets, we find the most probable values of parameters in the analytical functions given the observed eccentricities. These analytical functions can be used in theoretical investigations or as priors for the eccentricity distribution when analyzing short-period planets. As the measurement precision improves and sample size increases, the method can be applied to more complex parametrizations for the underlying distribution of eccentricity for extrasolar planetary systems.

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 present a Bayesian surrogate model for the analysis of periodic or quasi-periodic time series data. We describe a computationally efficient implementation that enables Bayesian model comparison. We apply this model to simulated and real exoplanet observations. We discuss the results and demonstrate some of the challenges for applying our surrogate model to realistic exoplanet data sets. In particular, we find that analyses of real world data should pay careful attention to the effects of uneven spacing of observations and the choice of prior for the "jitter" parameter.

We present results from the Radio Interferometric Planet (RIPL) search for compan- ions to the nearby star GJ 896A. We present 11 observations over 4.9 years. Fitting astrometric parameters to the data reveals a residual with peak-to-peak amplitude of ~ 3 mas in right ascension. This residual is well-fit by an acceleration term of 0.458 \pm 0.032 mas/y^2. The parallax is fit to an accuracy of 0.2 mas and the proper motion terms are fit to accuracies of 0.01 mas/y. After fitting astrometric and acceleration terms residuals are 0.26 mas in each coordinate, demonstrating that stellar jitter does not limit the ability to carry out radio astrometric planet detection and characterization. The acceleration term originates in part from the companion GJ 896B but the amplitude of the acceleration in declination is not accurately predicted by the orbital model. The acceleration sets a mass upper limit of 0.15 MJ at a semi-major axis of 2 AU for a planetary companion to GJ 896A. For semi-major axes between 0.3 and 2 AU upper limits are determined by the maximum angular separation; the upper limits scale from the minimum value in proportion to the inverse of the radius. Upper limits at larger radii are set by the acceleration and scale as the radius squared. An improved solution for the stellar binary system could improve the exoplanet mass sensitivity by an order of magnitude.

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.

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.

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.

If the two planets in the HAT-P-13 system are coplanar, the orbital states provide a probe of the internal planetary structure. Previous analyses of radial velocity and transit timing data of the system suggested that the observational constraints on the orbital states were rather small. We reanalyze the available data, treating the jitter as an unknown MCMC parameter, and find that a wide range of jitter values are plausible, hence the system parameters are less well constrained than previously suggested. For slightly increased levels of jitter ($\sim 4.5\,m\,s^{-1}$) the eccentricity of the inner planet can be in the range $0<e_{inner}<0.07$, the period and eccentricity of the outer planet can be $440<P_{outer}<470$ days and $0.55<e_{outer}<0.85$ respectively, while the relative pericenter alignment, $\eta$, of the planets can take essentially any value $-180^{\circ}<\eta<+180^{\circ}$. It is therefore difficult to determine whether $e_{inner}$ and $\eta$ have evolved to a fixed-point state or a limit cycle, or to use $e_{inner}$ to probe the internal planetary structure. We perform various transit timing variation (TTV) analyses, demonstrating that current constraints merely restrict $e_{outer}<0.85$, and rule out relative planetary inclinations within $\sim 2^{\circ}$ of $i_{rel}=90^{\circ}$, but that future observations could significantly tighten the restriction on both these parameters. We demonstrate that TTV profiles can readily distinguish the theoretically favored inclinations of $i_{rel}=0^{\circ}\,&\,45^{\circ}$, provided that sufficiently precise and frequent transit timing observations of HAT-P-13b can be made close to the pericenter passage of HAT-P-13c. We note the relatively high probability that HAT-P-13c transits and suggest observational dates and strategies.

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.

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.

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.

This paper reports the detection and the measurements of occultations of the two transiting hot giant exoplanets Kepler-5b and Kepler-6b by their parent stars. The observations are obtained in the near infrared with Spitzer Space Telescope and at optical wavelengths by combining more than a year of Kepler photometry. The investigation consists of constraining the eccentricities of these systems and of obtaining broad band emergent spectra for individual planets. For both targets, the occultations are detected at 3 sigma level at each wavelength with mid-occultation times consistent with circular orbits. The brightness temperatures of these planets are deduced from the infrared observations and reach T=1930+/-100K and T=1660+/-120K for Kepler-5b and Kepler-6b respectively. We measure optical geometric albedos A_g in the Kepler bandpass and find A_g=0.12+/-0.04 for Kepler-5b and A_g=0.11+/-0.04 for Kepler-6b leading to an upper limit for the Bond albedo of A_B < 0.17 in both cases. The observations for both planets are best described by models for which most of the incident energy is redistributed on the dayside, with only less than 10% of the absorbed stellar flux redistributed to the night side of these planets. The data for Kepler-5b favor a model without a temperature inversion, whereas for Kepler-6b they do not allow distinguishing between models with and without inversion.

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 new radial velocities from Keck Observatory and both Newtonian and Keplerian solutions for the triple-planet system orbiting HD 37124. The orbital solution for this system has improved dramatically since the third planet was first reported in Vogt et al. 2005 with an ambiguous orbital period. We have resolved this ambiguity, and the outer two planets have an apparent period commensurability of 2:1. A dynamical analysis finds both resonant and non-resonant configurations consistent with the radial velocity data, and constrains the mutual inclinations of the planets to be less than about 30 degrees. We discuss HD 37124 in the context of the other 19 exoplanetary systems with apparent period commenserabilities, which we summarize in a table. We show that roughly one in three well-characterized multiplanet systems has a apparent low-order period commensuribility, which is more than would naively be expected if the periods of exoplanets in known multiplanet systems were drawn randomly from the observed distribution of planetary orbital periods.

Recent observations have revealed two new classes of planetary orbits. Rossiter- Mclaughlin (RM) measurements have revealed hot Jupiters in high-obliquity orbits. In addition, direct-imaging has discovered giant planets at large (~ 100 AU) separations via direct-imaging technique. Simple-minded disk-migration scenarios are inconsistent with the high-inclination (and even retrograde) orbits as seen in recent RM measurements. Furthermore, forming giant planets at large semi-major axis (a) may be challenging in the core-accretion paradigm. We perform many N-body simulations to explore the two above-mentioned orbital architectures. Planet-planet scattering in a multi-planet system can naturally excite orbital inclinations. Planets can also get scattered to large distances. Large-a planetary orbits created from planet-planet scattering are expected to have high eccentricities (e). Theoretical models predict that the observed long-period planets, such as Fomalhaut-b have moderate e \approx 0.3. Interestingly, these are also in systems with disks. We find that if a massive-enough outer disk is present, a scattered planet may be circularized at large a via dynamical friction from the disk and repeated scattering of the disk particles.

We present a new short-period brown dwarf candidate around the star TYC 1240-00945-1. This candidate was discovered in the first year of the Multi-object APO Radial Velocity Exoplanets Large-area Survey (MARVELS), which is part of the third phase of the Sloan Digital Sky Survey (SDSS-III), and we designate the brown dwarf as MARVELS-1b. MARVELS uses the technique of dispersed fixed-delay interferometery to simultaneously obtain radial velocity measurements for 60 objects per field using a single, custom-built instrument that is fiber fed from the SDSS 2.5-m telescope. From our 20 radial velocity measurements spread over a ~370 d time baseline, we derive a Keplerian orbital fit with semi-amplitude K=2.533+/-0.025 km/s, period P=5.8953+/-0.0004 d, and eccentricity consistent with circular. Independent follow-up radial velocity data confirm the orbit. Adopting a mass of 1.37+/-0.11 M_Sun for the slightly evolved F9 host star, we infer that the companion has a minimum mass of 28.0+/-1.5 M_Jup, a semimajor axis 0.071+/-0.002 AU assuming an edge-on orbit, and is probably tidally synchronized. We find no evidence for coherent instrinsic variability of the host star at the period of the companion at levels greater than a few millimagnitudes. The companion has an a priori transit probability of ~14%. Although we find no evidence for transits, we cannot definitively rule them out for companion radii ~<1 R_Jup.

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.

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.

Both ground and space-based transit observatories are poised to significantly increase the number of known transiting planets and the number of precisely measured transit times. The variation in a planet's transit times may be used to infer the presence of additional planets. Deducing the masses and orbital parameters of such planets from transit time variations (TTVs) alone is a rich and increasingly relevant dynamical problem. In this work, we evaluate the extent of the degeneracies in this process, systematically explore the dependence of TTV signals on several parameters and provide phase space plots that could aid observers in planning future observations. Our explorations are focused on a likely-to-be prevalent situation: a known transiting short-period Neptune or Jupiter-sized planet and a suspected external low-mass perturber on a nearly-coplanar orbit. Through approximately 10^7 N-body simulations, we demonstrate how TTV signal amplitudes may vary by orders of magnitude due to slight variations in any one orbital parameter (0.001 AU in semimajor axis, 0.005 in eccentricity, or a few degrees in orbital angles), and quantify the number of consecutive transit observations necessary in order to obtain a reasonable opportunity to characterize the unseen planet (approximately greater or equal to 50 observations). Planets in or near period commensurabilities of the form p:q, where p < 21 and q < 4, produce distinct TTV signatures, regardless of whether the planets are actually locked in a mean motion resonance. We distinguish these systems from the secular systems in our explorations. Additionally, we find that computing the autocorrelation function of a TTV signal can provide a useful diagnostic for identifying possible orbits for additional planets and suggest that this method could aid integration of TTV signals in future studies of particular exosystems.

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.

We investigate potential biases in the measurements of exoplanet orbital parameters obtained from radial velocity observations for single-planet systems. We create a mock catalog of radial velocity data, choosing input planet masses, periods, and observing patterns from actual radial velocity surveys and varying input eccentricities. We apply Markov Chain Monte Carlo (MCMC) simulations and compare the resulting orbital parameters to the input values. We find that a combination of the effective signal-to-noise ratio of the data, the maximal gap in phase coverage, and the total number of periods covered by observations is a good predictor of the quality of derived orbit parameters. As eccentricity is positive definite, we find that eccentricities of planets on nearly circular orbits are preferentially overestimated, with typical bias of 1-2 times the median eccentricity uncertainty in a survey (e.g., 0.04 in the Butler et al. 2006 catalog). When performing population analysis, we recommend using the mode of the marginalized posterior eccentricity distribution to minimize potential biases. While the Butler et al. (2006) catalog reports eccentricities below 0.05 for just 17% of single-planet systems, we estimate that the true fraction of e < 0.05 orbits is about f(0.05)=38\pm 9%. For planets with P > 10 days, we find f(0.05)=28\pm 8% versus 10% from Butler et al. (2006). These planets either never acquired a large eccentricity or were circularized following any significant eccentricity excitation.

We report radial velocity measurements of the G-type subgiants 24 Sextanis (=HD90043) and HD200964. Both are massive, evolved stars that exhibit periodic variations due to the presence of a pair of Jovian planets. Photometric monitoring with the T12 0.80m APT at Fairborn Observatory demonstrates both stars to be constant in brightness to <= 0.002 mag, thus strengthening the planetary interpretation of the radial velocity variations. 24 Sex b,c have orbital periods of 453.8 days and 883~days, corresponding to semimajor axes 1.333 AU and 2.08 AU, and minimum masses (Msini) 1.99 Mjup and 0.86 Mjup, assuming a stellar mass 1.54 Msun. HD200964 b,c have orbital periods of 613.8 days and 825 days, corresponding to semimajor axes 1.601 AU and 1.95 AU, and minimum masses 1.85 Mjup and 0.90 Mjup, assuming M* = 1.44 Msun. We also carry out dynamical simulations to properly account for gravitational interactions between the planets. Most, if not all, of the dynamically stable solutions include crossing orbits, suggesting that each system is locked in a mean motion resonance that prevents close encounters and provides long-term stability. The planets in the 24 Sex system likely have a period ratio near 2:1, while the HD200964 system is even more tightly packed with a period ratio close to 4:3. However, we caution that further radial velocity observations and more detailed dynamical modelling will be required to provide definitive and unique orbital solutions for both cases, and to determine whether the two systems are truly resonant.

We report the first extrasolar planet observations from the 10.4-m Gran Telescopio Canarias (GTC), currently the world's largest, fully steerable, single-aperture optical telescope. We used the OSIRIS tunable filter imager on the GTC to acquire high-precision, narrow-band photometry of the transits of the giant exoplanets, TrES-2b and TrES-3b. We obtained near-simultaneous observations in two near-infrared (NIR) wavebands (790.2 and 794.4 +/- 2.0 nm) specifically chosen to avoid water vapor absorption and skyglow so as to minimize the atmospheric effects that often limit the precision of ground-based photometry. Our results demonstrate a very-high photometric precision with minimal atmospheric contamination despite relatively poor atmospheric conditions and some technical problems with the telescope. We find the photometric precision for the TrES-2 observations to be 0.343 and 0.412 mmag for the 790.2 and 794.4 nm light curves, and the precision of the TrES-3 observations was found to be 0.470 and 0.424 mmag for the 790.2 and 794.4 nm light curves. We also discuss how future follow-up observations of transiting planets with this novel technique can contribute to the characterization of Neptune- and super-Earth-size planets to be discovered by space-based missions like CoRoT and Kepler, as well as measure atmospheric properties of giant planets, such as the strength of atmospheric absorption features.

The discovery of multi-planet extrasolar systems has kindled interest in using their orbital evolution as a probe of planet formation. Accurate descriptions of planetary orbits identify systems which could hide additional planets or be in a special dynamical state, and inform targeted follow-up observations. We combine published radial velocity data with Markov Chain Monte Carlo analyses in order to obtain an ensemble of masses, semimajor axes, eccentricities and orbital angles for each of 5 dynamically active multi-planet systems: HD 11964, HD 38529, HD 108874, HD 168443, and HD 190360. We dynamically evolve these systems using 52,000 long-term N-body integrations that sample the full range of possible line-of-sight and relative inclinations, and we report on the system stability, secular evolution and the extent of the resonant interactions. We find that planetary orbits in hierarchical systems exhibit complex dynamics and can become highly eccentric and maybe significantly inclined. Additionally we incorporate the effects of general relativity in the long-term simulations and demonstrate that can qualitatively affect the dynamics of some systems with high relative inclinations. The simulations quantify the likelihood of different dynamical regimes for each system and highlight the dangers of restricting simulation phase space to a single set of initial conditions or coplanar orbits.

We perform numerical calculations of the expected transit timing variations (TTVs) induced on a Hot-Jupiter by an Earth-mass perturber. Motivated by the recent discoveries of retrograde transiting planets, we concentrate on an investigation of the effect of varying relative planetary inclinations, up to and including completely retrograde systems. We find that planets in low order (E.g. 2:1) mean-motion resonances (MMRs) retain approximately constant TTV amplitudes for $0<\,^{\circ}i<170\,^{\circ}$, only reducing in amplitude for $i > 170\,^{\circ}$. Systems in higher order MMRs (E.g. 5:1) increase in TTV amplitude as inclinations increase towards $45\,^{\circ}$, becoming approximately constant for $45 < i < 135$, and then declining for $i > 135\,^{\circ}$. Planets away from resonance slowly decrease in TTV amplitude as inclinations increase from 0 to 180, where-as planets adjacent to resonances can exhibit a huge range of variability in TTV amplitude as a function of both eccentricity and inclination. For highly retrograde systems ($135\,^{\circ} < i \leq 180\,^{\circ}$), TTV signals will be undetectable across almost the entirety of parameter space, with the exceptions occurring when the perturber has high eccentricity or is very close to a MMR. This high inclination decrease in TTV amplitude (on and away from resonance) is important for the analysis of the known retrograde and multi-planet transiting systems, as inclination effects need to be considered if TTVs are to be used to exclude the presence of any putative planetary companions: absence of evidence is not evidence of absence.

The recent discoveries of massive planets on ultra-wide orbits of HR 8799 (Marois et al. 2008) and Fomalhaut (Kalas et al. 2008) present a new challenge for planet formation theorists. Our goal is to figure out which of three giant planet formation mechanisms--core accretion (with or without migration), scattering from the inner disk, or gravitational instability--could be responsible for Fomalhaut b, HR 8799 b, c and d, and similar planets discovered in the future. This paper presents the results of numerical experiments comparing the long-period planet formation efficiency of each possible mechanism in model A star, G star and M star disks. First, a simple core accretion simulation shows that planet cores forming beyond 35 AU cannot reach critical mass, even under the most favorable conditions one can construct. Second, a set of N-body simulations demonstrates that planet-planet scattering does not create stable, wide-orbit systems such as HR 8799. Finally, a linear stability analysis verifies previous work showing that global spiral instabilities naturally arise in high-mass disks. We conclude that massive gas giants on stable orbits with semimajor axes greater than 35 AU form by gravitational instability in the disk. We recommend that observers examine the planet detection rate as a function of stellar age, controlling for planet dimming with time. If planet detection rate is found to be independent of stellar age, it would confirm our prediction that gravitational instability is the dominant mode of producing detectable planets on wide orbits. We also predict that the occurrence ratio of long-period to short-period gas giants should be highest for M dwarfs due to the inefficiency of core accretion and the expected small fragment mass in their disks.

We present the results of a transcontinental campaign to observe the 2009 June 5 transit of the exoplanet HD 80606b. We report the first detection of the transit ingress, revealing the transit duration to be 11.64 +/- 0.25 hr and allowing more robust determinations of the system parameters. Keck spectra obtained at midtransit exhibit an anomalous blueshift, giving definitive evidence that the stellar spin axis and planetary orbital axis are misaligned. The Keck data show that the projected spin-orbit angle is between 32-87 deg with 68.3% confidence and between 14-142 deg with 99.73% confidence. Thus the orbit of this planet is not only highly eccentric (e=0.93), but is also tilted away from the equatorial plane of its parent star. A large tilt had been predicted, based on the idea that the planet's eccentric orbit was caused by the Kozai mechanism. Independently of the theory, it is noteworthy that all 3 exoplanetary systems with known spin-orbit misalignments have massive planets on eccentric orbits, suggesting that those systems migrate differently than lower-mass planets on circular orbits.

Currently, over forty transiting planets have been discovered by ground-based photometric surveys, and space-based missions like Kepler and CoRoT are expected to detect hundreds more. Follow-up photometric observations from the ground will play an important role in constraining both orbital and physical parameters for newly discovered planets, especially those with small radii (R_p less than approximately 4 Earth radii) and/or intermediate to long orbital periods (P greater than approximately 30 days). Here, we simulate transit light curves from Kepler-like photometry and ground-based observations in the near-infrared (NIR) to determine how jointly modeling space-based and ground-based light curves can improve measurements of the transit duration and planet-star radius ratio. We find that adding observations of at least one ground-based transit to space-based observations can significantly improve the accuracy for measuring the transit duration and planet-star radius ratio of small planets (R_p less than approximately 4 Earth radii) in long-period (~1 year) orbits, largely thanks to the reduced effect of limb darkening in the NIR. We also demonstrate that multiple ground-based observations are needed to gain a substantial improvement in the measurement accuracy for small planets with short orbital periods (~3 days). Finally, we consider the role that higher ground-based precisions will play in constraining parameter measurements for typical Kepler targets. Our results can help inform the priorities of transit follow-up programs (including both primary and secondary transit of planets discovered with Kepler and CoRoT), leading to improved constraints for transit durations, planet sizes, and orbital eccentricities.

Radio astrometry of nearby, low-mass stars has the potential to be a powerful tool for the discovery and characterization of planetary companions. We present a Very Large Array survey of 172 active M dwarfs at distances of less than 10 pc. Twenty nine stars were detected with flux densities greater than 100 microJy. We observed 7 of these stars with the Very Long Baseline Array at milliarcsecond resolution in three separate epochs. With a detection threshold of 500 microJy in images of sensitivity 1 sigma ~ 100 microJy, we detected three stars three times (GJ 65B, GJ896A, GJ 4247), one star twice (GJ 285), and one star once (GJ 803). Two stars were undetected (GJ 412B and GJ 1224). For the four stars detected in multiple epochs, residuals from the optically-determined proper motions have an rms deviation of ~0.2 milliarcseconds, consistent with statistical noise limits. Combined with previous optical astrometry, these residuals provide acceleration upper limits that allow us to exclude planetary companions more massive than 3-6 M_Jup at a distance of ~1 AU with a 99% confidence level.

We present new precision radial velocities and a three-planet Keplerian orbit fit for the V = 8.5, G5 V star HIP 14810. We began observing this star at Keck Observatory as part of the N2K Planet Search Project. Wright et al. (2007) announced the inner two planets to this system, and subsequent observations have revealed the outer planet planet and the proper orbital solution for the middle planet. The planets have minimum masses of 3.9, 1.3, and 0.6 M_Jup and orbital periods of 6.67, 147.7, and 952 d, respectively. We have numerically integrated the family of orbital solutions consistent with the data and find that they are stable for at least 10^6 yr. Our photometric search shows that the inner planet does not transit.

The observed extrasolar planets possess both large masses (with a median M sin i of 1.65 MJ) and a wide range in orbital eccentricity (0 < e < 0.94). As planets are thought to form in circumstellar disks, one important question in planet formation is determining whether, and to what degree, a gaseous disk affects an eccentric planet's orbit. Recent studies have probed the interaction between a disk and a terrestrial planet on an eccentric orbit, and the interaction between a disk and a gas giant on a nearly circular orbit, but little is known about the interaction between a disk and an eccentric gas giant. Such a scenario could arise due to scattering while the disk is still present, or perhaps through planet formation via gravitational instability. We fill this gap with simulations of eccentric, massive (gap-forming) planets in disks using the hydrodynamical code FARGO. Although the long-term orbital evolution of the planet depends on disk properties, including the boundary conditions used, the disk always acts to damp eccentricity when the planet is released into the disk. This eccentricity damping takes place on a timescale of 40 years, 15 times faster than the migration timescale.

Astrometry can detect rocky planets in a broad range of masses and orbital distances and measure their masses and three-dimensional orbital parameters, including eccentricity and inclination, to provide the properties of terrestrial planets. The masses of both the new planets and the known gas giants can be measured unambiguously, allowing a direct calculation of the gravitational interactions, both past and future. Such dynamical interactions inform theories of the formation and evolution of planetary systems, including Earth-like planets. Astrometry is the only technique technologically ready to detect planets of Earth mass in the habitable zone (HZ) around solar-type stars within 20 pc. These Earth analogs are close enough for follow-up observations to characterize the planets by infrared imaging and spectroscopy with planned future missions such as the James Webb Space Telescope (JWST) and the Terrestrial Planet Finder/Darwin. Employing a demonstrated astrometric precision of 1 microarcsecond and a noise floor under 0.1 micro-arcseconds, SIM Lite can make multiple astrometric measurements of the nearest 60 F-, G-, and K-type stars during a five-year mission. SIM Lite directly tests theories of rocky planet formation and evolution around Sun-like stars and identifies the nearest potentially habitable planets for later spaceborne imaging, e.g., with Terrestrial Planet Finder and Darwin. SIM was endorsed by the two recent Decadal Surveys and it meets the highest-priority goal of the 2008 AAAC Exoplanet Task Force.

We ask if Earth-like planets (terrestrial mass and habitable-zone orbit) can be detected in multi-planet systems, using astrometric and radial velocity observations. We report here the preliminary results of double-blind calculations designed to answer this question.

We present the latest velocities for 10 multi-planet systems, including a re-analysis of archival Keck and Lick data, resulting in improved velocities that supersede our previously published measurements. We derive updated orbital fits for ten Lick and Keck systems, including two systems (HD 11964, HD 183263) for which we provide confirmation of second planets only tentatively identified elsewhere, and two others (HD 187123, and HD 217107) for which we provide a major revision of the outer planet's orbit. We compile orbital elements from the literature to generate a catalog of the 28 published multiple-planet systems around stars within 200 pc. From this catalog we find several intriguing patterns emerging: - Including those systems with long-term radial velocity trends, at least 28% of known planetary systems appear to contain multiple planets. - Planets in multiple-planet systems have somewhat smaller eccentricities than single planets. - The distribution of orbital distances of planets in multi-planet systems and single planets are inconsistent: single-planet systems show a pile-up at P ~ 3 days and a jump near 1 AU, while multi-planet systems show a more uniform distribution in log-period. In addition, among all planetary systems we find: - There may be an emerging, positive correlation between stellar mass and giant-planet semi-major axis. - Exoplanets more massive than Jupiter have eccentricities broadly distributed across 0 < e < 0.5, while lower-mass exoplanets exhibit a distribution peaked near e = 0.

Direct imaging searches have begun to detect planetary and brown dwarf companions and to place constraints on the presence of giant planets at large separations from their host star. This work helps to motivate such planet searches by predicting a population of young giant planets that could be detectable by direct imaging campaigns. Both the classical core accretion and the gravitational instability model for planet formation are hard-pressed to form long-period planets in situ. Here, we show that dynamical instabilities among planetary systems that originally formed multiple giant planets much closer to the host star could produce a population of giant planets at large (~100 AU - 100000 AU) separations. We estimate the limits within which these planets may survive, quantify the efficiency of gravitational scattering into both stable and unstable wide orbits, and demonstrate that population analyses must take into account the age of the system. We predict that planet scattering creates a population of detectable giant planets on wide orbits that decreases in number on timescales of ~10 Myr. We demonstrate that several members of such populations should be detectable with current technology, quantify the prospects for future instruments, and suggest how they could place interesting constraints on planet formation models.

[Abridged] We present the results of a highly parallel Kepler equation solver using the Graphics Processing Unit (GPU) on a commercial nVidia GeForce 280GTX and the "Compute Unified Device Architecture" programming environment. We apply this to evaluate a goodness-of-fit statistic (e.g., chi^2) for Doppler observations of stars potentially harboring multiple planetary companions (assuming negligible planet-planet interactions). We tested multiple implementations using single precision, double precision, pairs of single precision, and mixed precision arithmetic. We find that the vast majority of computations can be performed using single precision arithmetic, with selective use of compensated summation for increased precision. However, standard single precision is not adequate for calculating the mean anomaly from the time of observation and orbital period when evaluating the goodness-of-fit for real planetary systems and observational data sets. Using all double precision, our GPU code outperforms a similar code using a modern CPU by a factor of over 60. Using mixed-precision, our GPU code provides a speed-up factor of over 600, when evaluating N_sys > 1024 models planetary systems each containing N_pl = 4 planets and assuming N_obs = 256 observations of each system. We conclude that modern GPUs also offer a powerful tool for repeatedly evaluating Kepler's equation and a goodness-of-fit statistic for orbital models when presented with a large parameter space.

HD 69830 exhibits radial velocity variations attributed to three planets as well as infrared emission attributed to a warm debris disk. Previous studies have developed models for the planet migration and mass growth (Alibert et al. 2006) and the replenishment of warm grains (Wyatt et al. 2007). We perform n-body integrations in order to explore the implications of these models for: 1) the excitation of planetary eccentricity, 2) the accretion and clearing of a putative planetesimal disk, 3) the distribution of planetesimal orbits following migration, and 4) the implications for the origin of the IR emission. We find that: i) It is not possible to explain the observed planetary eccentricities (e~0.1) purely as the result of planetary perturbations during migration unless the planetary system is nearly face-on. ii) The rate of accretion of planetesimals onto planets in our n-body simulations is significantly different to that assumed in the semi-analytic models, suggesting that one cannot successfully treat planetesimal accretion in the simplified manner of Alibert et al. (2006). iii) Eccentricity damping of planetesimals does not act as an insurmountable obstacle to the existence of an excited eccentric disk: All simulations result in at least 25 Earth-masses of material remaining bound in the region ~1-9 AU, even after all three planets have migrated through the region. iv) Gas drag works to size-sort the planetesimals, with the largest bodies preferentially occupying the highest eccentricity and longest-lived orbits. Further work will be required to understand whether these eccentricity distributions are high enough to explain the level of dust emission observed despite mass loss via steady state collisional evolution. [abridged]

The eccentricity evolution of multiple planet systems can provide valuable constraints on planet formation models. Unfortunately, the inevitable uncertainties in the current orbital elements can lead to significant ambiguities in the nature of the secular evolution. Integrating any single set of orbital elements inadequately describes the full range of secular evolutions consistent with current observations. Thus, we combine radial velocity observations of HD 12661 with Markov Chain Monte Carlo sampling to generate ensembles of initial conditions for direct n-body integrations. We find that any mean motion resonances are quite weak and do not significantly impact the secular evolution, and that current observations indicate circulation or large amplitude libration of the periapses. The eccentricity of the outer planet undergoes large oscillations for nearly all of the allowed two-planet orbital solutions. This type of secular evolution would arise if planet c had been impulsively perturbed, perhaps due to strong scattering of an additional planet that was subsequently accreted onto the star. Finally, we note that the secular evolution implied by the current orbital configuration implies that planet c spends ~96% of the time following an orbit more eccentric than that presently observed. Either this system is being observed during a relatively rare state, or additional planets are affecting the observed radial velocities and/or the system's secular eccentricity evolution.

(Abridged) In planetary systems with two or more giant planets, dynamical instabilities can lead to collisions or ejections through strong planet--planet scattering. Previous studies for simple initial configurations with two equal-mass planets revealed some discrepancies between the results of numerical simulations and the observed orbital elements of extrasolar planets. Here, we show that simulations with two unequal mass planets starting on nearly circular orbits predict a reduced frequency of collisions and a broader range of final eccentricities. The two-planet scattering model can easily reproduce the observed eccentricities with a plausible distribution of planet mass ratios. Further, the two-planet scattering model predicts a maximum eccentricity of about 0.8, independent of the distribution of planet mass ratios, provided that both planets are initially place on nearly circular orbits. This compares favorably with current observations and will be tested by future planet discoveries. The combination of planet--planet scattering and tidal circularization may be able to explain the existence of some giant planets with very short period orbits. Orbital migration due to planet scattering could play an important role in explaining the increased rate of giant planets with very short period orbits. We also re-examine and discuss various possible correlations between eccentricities and other properties of observed extrasolar planets. We demonstrate that the observed distribution of planet masses, orbital periods, and eccentricities can provide constraints for models of planet formation and evolution.

Observations in the past decade have revealed extrasolar planets with a wide range of orbital semimajor axes and eccentricities. Based on the present understanding of planet formation via core accretion and oligarchic growth, we expect that giant planets often form in closely packed configurations. While the protoplanets are embedded in a protoplanetary gas disk, dissipation can prevent eccentricity growth and suppress instabilities from becoming manifest. However, once the disk dissipates, eccentricities can grow rapidly, leading to close encounters between planets. Strong planet--planet gravitational scattering could produce both high eccentricities and, after tidal circularization, very short-period planets, as observed in the exoplanet population. We present new results for this scenario based on extensive dynamical integrations of systems containing three giant planets, both with and without residual gas disks. We assign the initial planetary masses and orbits in a realistic manner following the core accretion model of planet formation. We show that, with realistic initial conditions, planet--planet scattering can reproduce quite well the observed eccentricity distribution. Our results also make testable predictions for the orbital inclinations of short-period giant planets formed via strong planet scattering followed by tidal circularization.

With the recent discoveries of hundreds of extrasolar planets, the search for planets like Earth and life in the universe, is quickly gaining momentum. In the future, large space observatories could directly detect the light scattered from rocky planets, but they would not be able to spatially resolve a planet's surface. Using reflectance models and real cloud data from satellite observations, here we show that, despite Earth's dynamic weather patterns, the light scattered by the Earth to a hypothetical distant observer as a function of time contains sufficient information to accurately measure Earth's rotation period. This is because ocean currents and continents result in relatively stable averaged global cloud patterns. The accuracy of these measurements will vary with the viewing geometry and other observational constraints. If the rotation period can be measured with accuracy, data spanning several months could be coherently combined to obtain spectroscopic information about individual regions of the planetary surface. Moreover, deviations from a periodic signal can be used to infer the presence of relatively short-live structures in its atmosphere (i.e., clouds). This could provide a useful technique for recognizing exoplanets that have active weather systems, changing on a timescale comparable to their rotation. Such variability is likely to be related to the atmospheric temperature and pressure being near a phase transition and could support the possibility of liquid water on the planet's surface.

The discovery of over 200 extrasolar planets with the radial velocity (RV) technique has revealed that many giant planets have large eccentricities, in striking contrast with most of the planets in the solar system and prior theories of planet formation. The realization that many giant planets have large eccentricities raises a fundamental question: ``Do terrestrial-size planets of other stars typically have significantly eccentric orbits or nearly circular orbits like the Earth?'' Here, we demonstrate that photometric observations of transiting planets could be used to characterize the orbital eccentricities for individual transiting planets, as well the eccentricity distribution for various populations of transiting planets (e.g., those with a certain range of orbital periods or physical sizes). Such characterizations can provide valuable constraints on theories for the excitation of eccentricities and tidal dissipation. We outline the future prospects of the technique given the exciting prospects for future transit searches, such as those to be carried out by the CoRoT and Kepler missions.

To achieve maximum planet yield for a given radial velocity survey, the observing strategy must be carefully considered. In particular, the adopted cadence can greatly affect the sensitivity to exoplanetary parameters such as period and eccentricity. Here we describe simulations which aim to maximise detections based upon the target parameter space of the survey.

High-precision radial velocity planet searches have surveyed over ~2000 nearby stars and detected over ~200 planets. While these same stars likely harbor many additional planets, they will become increasingly challenging to detect, as they tend to have relatively small masses and/or relatively long orbital periods. Therefore, observers are increasing the precision of their observations, continuing to monitor stars over decade timescales, and also preparing to survey thousands more stars. Given the considerable amounts of telescope time required for such observing programs, it is important use the available resources as efficiently as possible. Previous studies have found that a wide range of predetermined scheduling algorithms result in planet searches with similar sensitivities. We have developed adaptive scheduling algorithms which have a solid basis in Bayesian inference and information theory and also are computationally feasible for modern planet searches. We have performed Monte Carlo simulations of plausible planet searches to test the power of adaptive scheduling algorithms. Our simulations demonstrate that planet searches performed with adaptive scheduling algorithms can simultaneously detect more planets, detect less massive planets, and measure orbital parameters more accurately than comparable surveys using a non-adaptive scheduling algorithm. We expect that these techniques will be particularly valuable for the N2K radial velocity planet search for short-period planets as well as future astrometric planet searches with the Space Interferometry Mission which aim to detect terrestrial mass planets.

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

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.

The extraordinary astrometric accuracy of radio interferometry creates an important and unique opportunity for the discovery and characterization of exo-planets. Currently, the Very Long Baseline Array can routinely achieve better than 100 microarcsecond accuracy, and can approach 10 microarcsecond with careful calibration. We describe here RIPL, the Radio Interferometric PLanet search, a new program with the VLBA and the Green Bank 100 m telescope that will survey 29 low-mass, active stars over 3 years with sub-Jovian planet mass sensitivity at 1 AU. An upgrade of the VLBA bandwidth will increase astrometric accuracy by an order of magnitude. Ultimately, the colossal collecting area of the Square Kilometer Array could push astrometric accuracy to 1 microarcsecond, making detection and characterizaiton of Earth mass planets possible. RIPL and other future radio astrometric planet searches occupy a unique volume in planet discovery and characterization parameter space. The parameter space of astrometric searches gives greater sensitivity to planets at large radii than radial velocity searches. For the VLBA and the expanded VLBA, the targets of radio astrometric surveys are by necessity nearby, low-mass, active stars, which cannot be studied efficiently through the radial velocity method, coronagraphy, or optical interferometry. For the SKA, detection sensitivity will extend to solar-type stars. Planets discovered through radio astrometric methods will be suitable for characterization through extreme adaptive optics. The complementarity of radio astrometric techniques with other methods demonstrates that radio astrometry can play an important role in the roadmap for exoplanet discovery and characterization.

As many as 5 ice giants--Neptune-mass planets composed of 90% ice and rock and 10% hydrogen--are thought to form at heliocentric distances of 10-25 AU on closely packed orbits spaced ~5 Hill radii apart. Such oligarchies are ultimately unstable. Once the parent disk of planetesimals is sufficiently depleted, oligarchs perturb one another onto crossing orbits. We explore both the onset and the outcome of the instability through numerical integrations, including dynamical friction cooling of planets by a planetesimal disk whose properties are held fixed. To trigger instability and the ejection of the first ice giant in systems having an original surface density in oligarchs of Sigma ~ 1 g/cm^2, the disk surface density s must fall below 0.1 g/cm^2. Ejections are predominantly by Jupiter and occur within 10 Myr. To eject more than 1 oligarch requires s < 0.03 g/cm^2. Systems starting with up to 4 oligarchs in addition to Jupiter and Saturn can readily yield solar-system-like outcomes in which 2 surviving ice giants lie inside 30 AU and have their orbits circularized by dynamical friction. Our numerical simulations support the idea that planetary systems begin in more crowded and compact configurations, like those of shear-dominated oligarchies. In contrast to previous studies, we identify s < 0.1 Sigma as the regime relevant for understanding the evolution of the outer solar system, and we encourage future studies to concentrate on this regime while relaxing our assumption of a fixed planetesimal disk.

We determine stellar parameters for the M dwarf GJ 436 that hosts a Neptune-mass planet. We employ primarily spectral modeling at low and high resolution, examining the agreement between model and observed optical spectra of five comparison stars of type, M0-M3. Modeling high resolution optical spectra suffers from uncertainties in TiO transitions, affecting the predicted strengths of both atomic and molecular lines in M dwarfs. The determination of Teff, gravity, and metallicity from optical spectra remains at ~10%. As molecules provide opacity both in lines and as an effective continuum, determing molecular transition parameters remains a challenge facing models such as the PHOENIX series, best verified with high resolution and spectrophotometric spectra. Our analysis of GJ 436 yields an effective temperature of Teff = 3350 +/- 300 K and a mass of 0.44 Msun. New Doppler measurements for GJ 436 with a precision of 3 m/s taken during 6 years improve the Keplerian model of the planet, giving a minimum mass, M sin i = 0.0713 Mjup = 22.6 Mearth, period, P = 2.6439 d, and e = 0.16 +/- 0.02. The noncircular orbit contrasts with the tidally circularized orbits of all close-in exoplanets, implying either ongoing pumping of eccentricity by a more distant companion, or a higher Q value for this low-mass planet. The velocities indeed reveal a long term trend, indicating a possible distant companion.

Theoretical studies predict that Trojans are likely a frequent byproduct of planet formation and evolution. We present a novel method of detecting Trojan companions to transiting extrasolar planets which involves comparing the time of central eclipse with the time of the stellar reflex velocity null. We demonstrate that this method offers the potential to detect terrestrial-mass Trojans using existing ground-based observatories. This method rules out Trojan companions to HD 209458b and HD 149026b more massive than ~13 Earth masses and \~25 Earth masses at a 99.9% confidence level. Such a Trojan would be dynamically stable, would not yet have been detected by photometric or spectroscopic monitoring, and would be unrecognizable from radial velocity observations alone. We outline the future prospects for this method, and show that the detection of a "Hot Trojan" of any mass would place a significant constraint on theories of orbital migration.

We derive detailed theoretical models for 1074 nearby stars from the SPOCS (Spectroscopic Properties of Cool Stars) Catalog. The California and Carnegie Planet Search has obtained high-quality echelle spectra of over 1000 nearby stars taken with the Hamilton spectrograph at Lick Observatory, the HIRES spectrograph at Keck, and UCLES at the Anglo Australian Observatory. A uniform analysis of the high-resolution spectra has yielded precise stellar parameters, enabling systematic error analyses and accurate theoretical stellar modeling. We have created a large database of theoretical stellar evolution tracks using the Yale Stellar Evolution Code (YREC) to match the observed parameters of the SPOCS stars. Our very dense grids of evolutionary tracks eliminate the need for interpolation between stellar evolutionary tracks and allow precise determinations of physical stellar parameters (mass, age, radius, size and mass of the convective zone, etc.). Combining our stellar models with the observed stellar atmospheric parameters and uncertainties, we compute the likelihood for each set of stellar model parameters separated by uniform time steps along the stellar evolutionary tracks. The computed likelihoods are used for a Bayesian analysis to derive posterior probability distribution functions for the physical stellar parameters of interest. We provide a catalog of physical parameters for 1074 stars that are based on a uniform set of high quality spectral observations, a uniform spectral reduction procedure, and a uniform set of stellar evolutionary models. We explore this catalog for various possible correlations between stellar and planetary properties, which may help constrain the formation and dynamical histories of other planetary systems.

Radial velocity (RV) planet searches are increasingly finding planets with small velocity amplitudes, with long orbital periods, or in multiple planet systems. Bayesian inference has the potential to improve the interpretation of existing observations, the planning of future observations and ultimately inferences concerning the overall population of planets. In recent years, the refinement of Markov chain Monte Carlo (MCMC) algorithms has made it practical to accurately characterize orbital parameters and their uncertainties from RV observations of single-planet and weakly interacting multiple-planet systems. Unfortunately, MCMC is not sufficient for Bayesian model selection, i.e., comparing the marginal posterior probability of models, as is necessary to determine how strongly the observational data favor a model with n+1 planets over a model with just n planets. Many of the obvious estimators for the marginal posterior probability suffer from poor convergence properties. We compare several estimators of the marginal likelihood and feature those that display desirable convergence properties based on the analysis of a sample data set for HD 88133b. We find that methods based on importance sampling are most efficient, provided that a good analytic approximation of the posterior probability distribution is available. We present a simple algorithm for using a sample from the posterior to construct a mixture distribution that approximates the posterior and can be used for importance sampling and Bayesian model selection. We conclude with some suggestions for the development and refinement of computationally efficient and robust estimators of marginal posterior probabilities.

We present a modification to the iodine cell Doppler technique that eliminates the need for an observed stellar template spectrum. For a given target star, we iterate toward a synthetic template spectrum beginning with an existing spectrum of a similar star. We then perturb the shape of this first-guess template to match the program observation of the target star taken through an iodine cell. The elimination of a separate template observation saves valuable telescope time, a feature that is ideally suited for the quick-look strategy employed by the ``Next 2000 Stars'' (N2K) planet search program. Tests using Keck/HIRES spectra indicate that synthetic templates yield a short-term precision of 3 m/s and a long-term, run-to-run precision of 5 m/s. We used this new Doppler technique to discover three new planets: a 1.5 Mjup planet in a 2.1375 d orbit around HD 86081; a 0.71 Mjup planet in circular, 26.73 d orbit around HD 224693; and a Saturn-mass planet in an 18.179 d orbit around HD 33283. The remarkably short period of HD 86081b bridges the gap between the extremely short-period planets detected in the OGLE survey and the 16 Doppler-detected hot jupiters (P < 15 d), which have an orbital period distribution that piles up at about three days. We have acquired photometric observations of two of the planetary host stars with the automated photometric telescopes at Fairborn Observatory. HD 86081 and HD 224693 both lack detectable brightness variability on their radial velocity periods, supporting planetary-reflex motion as the cause of the radial velocity variability. HD 86081 shows no evidence of planetary transits in spite of a 17.6% transit probability.

We report the detection of the first extrasolar planet, ET-1 (HD 102195b), using the Exoplanet Tracker (ET), a new generation Doppler instrument. The planet orbits HD 102195, a young star with solar metallicity that may be part of the local association. The planet imparts radial velocity variability to the star with a semiamplitude of $63.4\pm2.0$ m s$^{-1}$ and a period of 4.11 days. The planetary minimum mass ($m \sin i$) is $0.488\pm0.015$ $M_J$.

For centuries, our knowledge of planetary systems and ideas about planet formation were based on a single example, our solar system. During the last thirteen years, the discovery of ~170 planetary systems has ushered in a new era for astronomy. I review the surprising properties of extrasolar planetary systems and discuss how they are reshaping theories of planet formation. I focus on how multiple planet systems constrain the mechanisms proposed to explain the large eccentricities typical of extrasolar planets. I suggest that strong planet-planet scattering is common and most planetary systems underwent a phase of large eccentricities. I propose that a planetary system's final eccentricities may be strongly influenced by how much mass remains in a planetesimal disk after the last strong planet-planet scattering event.

Many of the known extrasolar planets are ``hot Jupiters,'' giant planets with orbital periods of just a few days. We use the observed distribution of hot Jupiters to constrain the location of its inner edge in the mass--period diagram. If we assume a slope corresponding to the classical Roche limit, then we find that the edge corresponds to a separation close to_twice_ the Roche limit, as expected if the planets started on highly eccentric orbits that were later circularized. In contrast, any migration scenario would predict an inner edge right at the Roche limit, which applies to planets approaching on nearly circular orbits. However, the current sample of hot Jupiters is not sufficient to provide a precise constraint simultaneously on both the location and slope of the inner edge.

Precise radial velocity measurements have led to the discovery of ~170 extrasolar planetary systems. Understanding the uncertainties in the orbital solutions will become increasingly important as the discovery space for extrasolar planets shifts to planets with smaller masses and longer orbital periods. The method of Markov chain Monte Carlo (MCMC) provides a rigorous method for quantifying the uncertainties in orbital parameters in a Bayesian framework (Ford 2005a). The main practical challenge for the general application of MCMC is the need to construct Markov chains which quickly converge. The rate of convergence is very sensitive to the choice of the candidate transition probability distribution function (CTPDF). Here we explain one simple method for generating alternative CTPDFs which can significantly speed convergence by one to three orders of magnitude. We have numerically tested dozens of CTPDFs with simulated radial velocity data sets to identify those which perform well for different types of orbits and suggest a set of CTPDFs for general application. Additionally, we introduce other refinements to the MCMC algorithm for radial velocity planets, including an improved treatment of the uncertainties in the radial velocity observations, an algorithm for automatically choosing step sizes, an algorithm for automatically determining reasonable stopping times, and the use of importance sampling for including the dynamical evolution of multiple planet systems. Together, these improvements make it practical to apply MCMC to multiple planet systems. We demonstrate the improvements in efficiency by analyzing a variety of extrasolar planetary systems.

The Space Interferometry Mission (SIM) is expected to make precise astrometric measurements that can be used to detect low mass planets around nearby stars. Since most nearby stars are members of multiple star systems, many stars will have a measurable acceleration due to their companion, which must be included when solving for astrometric parameters and searching for planetary perturbations. Additionally, many of the stars with one radial velocity planet show indications of additional planets. Therefore, astrometric surveys like SIM must be capable of detecting planets and measuring orbital parameters in systems with multiple stellar and/or planetary companions. We have conducted Monte Carlo simulations to investigate how the presence of multiple companions affects the sensitivity of an astrometric survey such as SIM. We find that the detection efficiency for planets in wide binary systems is relatively unaffected by the presence of a binary companion, if the planetary orbital period is less than half the duration of the astrometric survey. For longer orbital periods, there are significant reductions in the sensitivity of an astrometric survey. Additionally, we find that the signal required to detect a planet can be increased significantly due to the presence of an additional planet orbiting the same star. Fortunately, adding a modest number of precision radial velocity observations significantly improves the sensitivity for many multiple planet systems. Thus, the combination of radial velocity observations and astrometric observations by SIM will be a particularly valuable for studying multiple planet systems.

Many of the known extrasolar planets are ``hot Jupiters,'' giant planets with orbital periods of just a few days. We use the observed distribution of hot Jupiters to constrain the location of the ``inner edge'' and planet migration theory. If we assume the location of the inner edge is proportional to the Roche limit, then we find that this edge is located near twice the Roche limit, as expected if the planets were circularized from a highly eccentric orbit. If confirmed, this result would place significant limits on migration via slow inspiral. However, if we relax our assumption for the slope of the inner edge, then the current sample of hot Jupiters is not sufficient to provide a precise constraint on both the location and power law index of the inner edge.

Earth's deciduous plants have a sharp order-of-magnitude increase in leaf reflectance between approximately 700 and 750 nm wavelength. This strong reflectance of Earth's vegetation suggests that surface biosignatures with sharp spectral features might be detectable in the spectrum of scattered light from a spatially unresolved extrasolar terrestrial planet. We assess the potential of Earth's step-function-like spectroscopic feature, referred to as the "red edge", as a tool for astrobiology. We review the basic characteristics and physical origin of the red edge and summarize its use in astronomy: early spectroscopic efforts to search for vegetation on Mars and recent reports of detection of the red edge in the spectrum of Earthshine (i.e., the spatially integrated scattered light spectrum of Earth). We present Earthshine observations from Apache Point Observatory to emphasize that time variability is key to detecting weak surface biosignatures such as the vegetation red edge. We briefly discuss the evolutionary advantages of vegetation's red edge reflectance, and speculate that while extraterrestrial "light harvesting organisms" have no compelling reason to display the exact same red edge feature as terrestrial vegetation, they might have similar spectroscopic features at different wavelengths than terrestrial vegetation. This implies that future terrestrial-planet-characterizing space missions should obtain data that allow time-varying, sharp spectral features at unknown wavelengths to be identified. We caution that some mineral reflectance edges are similar in slope and strength to vegetation's red edge (albeit at different wavelengths); if an extrasolar planet reflectance edge is detected care must be taken with its interpretation.

Doppler spectroscopy has detected 136 planets around nearby stars. A major puzzle is why their orbits are highly eccentric, while all planets in our Solar System are on nearly circular orbits, as expected if they formed by accretion processes in a protostellar disk. Several mechanisms have been proposed to generate large eccentricities after planet formation, but so far there has been little observational evidence to support any particular one. Here we report that the current orbital configuration of the three giant planets around Upsilon Andromedae (ups And) provides evidence for a close dynamical interaction with another planet, now lost from the system. The planets started on nearly circular orbits, but chaotic evolution caused the outer planet (ups And d) to be perturbed suddenly into a higher-eccentricity orbit. The coupled evolution of the system then causes slow periodic variations in the eccentricity of the middle planet (ups And c). Indeed, we show that ups And c periodically returns to a very nearly circular state every 9000 years. Our analysis shows that strong planet-planet scattering, one of several mechanisms previously discussed for increasing orbital eccentricities, must have operated in this system.

Precise radial velocity measurements have led to the discovery of ~100 extrasolar planetary systems. We investigate the uncertainty in the orbital solutions that have been fit to these observations. Understanding these uncertainties will become more and more important as the discovery space for extrasolar planets shifts to longer and longer periods. While detections of short period planets can be rapidly refined, planets with long orbital periods will require observations spanning decades to constrain the orbital parameters precisely. Already in some cases, multiple distinct orbital solutions provide similarly good fits, particularly in multiple planet systems. We present a method for quantifying the uncertainties in orbital fits and addressing specific questions directly from the observational data rather than relying on best fit orbital solutions. This Markov chain Monte Carlo (MCMC) technique has the advantage that it is well suited for the high dimensional parameter spaces necessary for the multiple planet systems. We apply the MCMC technique to several extrasolar planetary systems, assessing the uncertainties in orbital elements for several systems. Our MCMC simulations demonstrate that for some systems there are strong correlations between orbital parameters and/or significant non-Gaussianities in parameter distributions, even though the measurement errors are nearly Gaussian. Once these effects are considered the actual uncertainties in orbital elements can be significantly larger or smaller than the published uncertainties. We also present simple applications of our methods such as predicting the times of possible transits for GJ 876.

The Space Interferometry Mission (SIM) will make precise astrometric measurements that can be used to detect planets around nearby stars. Since observational time will be extremely valuable, it is important to consider how the choice of the observing schedule influences the efficiency of SIM planet searches. We have conducted Monte Carlo simulations of astrometric observations to understand the effects of different scheduling algorithms. We find that the efficiency of planet searches is relatively insensitive to the observing schedule for most reasonable observing schedules.

The Space Interferometry Mission (SIM) will make precise astrometric measurements that can be used to detect planets around nearby stars. We have simulated SIM observations and estimated the ability of SIM to detect planets with given masses and orbital periods and measure their orbital elements. We combine these findings with an estimate of the mass and period distribution of planets determined from radial velocity surveys to predict the number and characteristics of planets SIM would likely find. Our predictions are based on extrapolating the mass distribution of known extrasolar planets by up to a factor of 100. This extrapolation provides the best prediction we can make of the actual number of planets that SIM will detect and characterize, but may substantially over- or underestimate the frequency of Earth-mass planets, especially if these form by a different mechanism than giant planets. We find that SIM is likely to detect 1-5 planets with masses less than 3 M_Earth (depending on mission parameters). SIM would measure masses and orbits with 30% accuracy for 0-2 of these planets, but is unlikely to measure orbits with 10% accuracy for more than one of them. SIM is likely to detect 5-25 planets with mass less than 20 M_Earth, measure masses and orbits with 30% accuracy for 2-12 of these, and measure masses and orbits with 10% accuracy for 2-8 such planets. SIM is likely to find 25-160 planets of all masses, depending on the observing strategy and mission lifetime.

The search for Earth-like extrasolar planets is in part motivated by the potential detection of spectroscopic biomarkers. Spectroscopic biomarkers are spectral features that are either consistent with life, indicative of habitability, or provide clues to a planet's habitability. Most attention so far has been given to atmospheric biomarkers, gases such as O2, O3, H2O, CO, and CH4. Here we discuss surface biomarkers. Surface biomarkers that have large, distinct, abrupt changes in their spectra may be detectable in an extrasolar planet's spectrum at wavelengths that penetrate to the planetary surface. Earth has such a surface biomarker: the vegetation "red edge" spectroscopic feature. Recent interest in Earth's surface biomarker has motivated Earthshine observations of the spatially unresolved Earth and two recent studies may have detected the vegetation red edge feature in Earth's hemispherically integrated spectrum. A photometric time series in different colors should help in detecting unusual surface features in extrasolar Earth-like planets.

For the first time in human history the possibility of detecting and studying Earth-like planets is on the horizon. Terrestrial Planet Finder (TPF), with a launch date in the 2015 timeframe, is being planned by NASA to find and characterize planets in the habitable zones of nearby stars. The mission Darwin from ESA has similar goals. The motivation for both of these space missions is the detection and spectroscopic characterization of extrasolar terrestrial planet atmospheres. Of special interest are atmospheric biomarkers--such as O2, O3, H2O, CO and CH4--which are either indicative of life as we know it, essential to life, or can provide clues to a planet's habitability. A mission capable of measuring these spectral features would also obtain sufficient signal-to-noise to characterize other terrestrial planet properties. For example, physical characteristics such as temperature and planetary radius can be constrained from low- resolution spectra. In addition, planet characteristics such as weather, rotation rate, presence of large oceans or surface ice, and existence of seasons could be derived from photometric measurements of the planet's variability. We will review the potential to characterize terrestrial planets beyond their spectral signatures. We will also discuss the possibility to detect strong surface biomarkers--such as Earth's vegetation red edge near 700 nm--that are different from any known atomic, molecular, or mineralogical signature.

Instabilities and strong dynamical interactions between multiple giant planets have been proposed as a possible explanation for the surprising orbital properties of extrasolar planetary systems. In particular, dynamical instabilities seem to provide a natural mechanism for producing the highly eccentric orbits seen in many systems. Previously, we performed numerical integrations for the dynamical evolution of planetary systems containing two giant planets of equal masses initially in nearly circular orbits very close to the dynamical stability limit. We found the ratio of collisions to ejections in these simulations was greater than the ratio of circular orbits to eccentric orbits among the known extrasolar planets. Further, the mean eccentricity of the planets remaining after an ejection was larger than the mean eccentricity of the known extrasolar planets. Recently, we have performed additional integrations, generalizing to consider two planets of unequal masses. Our new simulations reveal that the two-planet scattering model can produce a distribution of eccentricities consistent with the observed eccentricity distribution for plausible mass distributions. Additionally, this model predicts a maximum eccentricity of about 0.8, in agreement with observations. Early results from simulations of three equal-mass planets also reveal a reduced frequency of collisions and a broad range of final eccentricities for the retained inner planet.

The light scattered by an extrasolar Earth-like planet's surface and atmosphere will vary in intensity and color as the planet rotates; the resulting light curve will contain information about the planet's properties. Since most of the light comes from a small fraction of the planet's surface, the temporal flux variability can be quite significant, $\sim$ 10-100%. In addition, for cloudless Earth-like extrasolar planet models, qualitative changes to the surface (such as ocean fraction, ice cover) significantly affect the light curve. Clouds dominate the temporal variability of the Earth but can be coherent over several days. In contrast to Earth's temporal variability, a uniformly, heavily clouded planet (e.g. Venus), would show almost no flux variability. We present light curves for an unresolved Earth and for Earth-like model planets calculated by changing the surface features. This work suggests that meteorological variability and the rotation period of an Earth-like planet could be derived from photometric observations. The inverse problem of deriving surface properties from a given light curve is complex and will require much more investigation.

The detection of massive planets orbiting nearby stars has become almost routine, but current techniques are as yet unable to detect terrestrial planets with masses comparable to the Earth's. Future space-based observatories to detect Earth-like planets are being planned. Terrestrial planets orbiting in the habitable zones of stars-where planetary surface conditions are compatible with the presence of liquid water-are of enormous interest because they might have global environments similar to Earth's and even harbor life. The light scattered by such a planet will vary in intensity and colour as the planet rotates; the resulting light curve will contain information about the planet's properties. Here we report a model that predicts features that should be discernible in light curves obtained by low-precision photometry. For extrasolar planets similar to Earth we expect daily flux variations up to hundreds of percent, depending sensitively on ice and cloud cover. Qualitative changes in surface or climate generate significant changes in the predicted light curves. This work suggests that the meteorological variability and the rotation period of an Earth-like planet could be derived from photometric observations. Other properties such as the composition of the surface (e.g., ocean versus land fraction), climate indicators (for example ice and cloud cover), and perhaps even signatures of Earth-like plant life could be constrained or possibly, with further study, even uniquely determined.

We study the kilohertz quasi-periodic oscillations (kHz QPOs) in the low-mass X-ray binaries 4U 1728-34, 4U 1608-52, and Aql X-1. Each source traces out a set of nearly parallel lines in a frequency vs. X-ray count rate diagram. We find that between two of these tracks, for similar QPO frequency, the source count rate can differ by up to a factor of ~ 4, whereas at the same time the rms amplitude of the kHz QPOs is only a factor of ~ 1.1 different. We also find that, for 4U 1608-52 and Aql X-1, the rms spectrum of the kHz QPOs does not depend upon which track the source occupies in the frequency vs. X-ray count rate diagram. Our results for 4U 1728-34, 4U 1608-52, and Aql X-1 are inconsistent with scenarios in which the properties of the kHz QPOs are only determined by the mass accretion rate through the disk, whereas X-ray count rate also depends upon other sources of energy that do not affect the QPOs.

High-resolution N-body simulations of four popular Cold Dark Matter cosmologies (LCDM, OCDM, QCDM, and tilted SCDM), each containing 10^5 clusters of galaxies in a cubic gigaparsec volume, are used to determine the evolution of the cluster mass function from z=3 to z=0. The large volume and high resolution of these simulations allow an accurate measure of the evolution of cosmologically important (but rare) massive clusters at high redshift. The simulated mass function is presented for cluster masses within several radii typically used observationally (R=0.5, 1.0, and 1.5 Mpc/h, both comoving and physical) in order to enable direct comparison with current and future observations. The simulated evolution is compared with current observations of massive clusters at redshifts 0.3<z<0.8. The Omega_m=1 tilted SCDM model, which exhibits very rapid evolution of the cluster abundance, produces too few clusters at z>0.3 and no massive clusters at z>0.5, in stark contradiction with observations. The Omega_m=0.3 models- LCDM, OCDM, and QCDM- all exhibit considerably weaker evolution and are consistent with current data. Among these low density models, OCDM evolves the least. These trends are enhanced at high redshift and can be used to discriminate between flat and open low density models. The simulated mass functions are compared with the Press-Schechter approximation. Standard Press-Schechter predicts too many low mass clusters at z=0, and too few clusters at higher redshift. We modify the approximation by a simple parameterization of the density contrast threshold for collapse, which has a redshift dependence. This modified Press-Schechter approximation provides a good fit to the simulated mass functions.

We report the discovery of a kilohertz quasi-periodic oscillation (kHz QPO) in the low-mass X-ray binary (LMXB) 2S 0918-549. The kHz QPO has a frequency of 1156+-9 Hz, a FWHM of 70+-20 Hz, and a fractional rms amplitude of 18%+-2%. We also detected for the first time a type I X-ray burst in this source. The compact object must therefore be a neutron star. Finally, we were able to classify the source as an atoll source exhibiting the full complement of spectral/timing states of this class (island and banana). From the peak burst flux an upper limit on the persistent luminosity can be derived of 0.5% of the Eddington luminosity, making 2S 0918-549 one of the least luminous LMXBs showing a kHz QPO. We compare the fractional rms amplitudes of the upper kHz QPO across the ensemble of LMXBs. We find a strong anticorrelation with luminosity. In LMXBs with luminosities ~100 times lower than those of Z-sources, the fractional rms amplitude is a factor ~10 larger.

Instabilities and strong dynamical interactions between several giant planets have been proposed as a possible explanation for the surprising orbital properties of extrasolar planetary systems. In particular, dynamical instabilities would seem to provide a natural mechanism for producing the highly eccentric orbits seen in many systems. Here we present results from numerical integrations for the dynamical evolution of planetary systems containing two identical giant planets in nearly circular orbits very close to the dynamical stability limit. We determine the statistical properties of the three main types of systems resulting from the development of an instability: systems containing one planet, following either a collision between the two initial planets, or the ejection of one of them to infinity, and systems containing two planets in a new, quasi-stable configuration. We discuss the implications of our results for the formation and evolution of observed extrasolar planetary systems. We conclude that the distributions of eccentricities and semimajor axes for observed systems cannot be explained easily by invoking dynamical interactions between two planets initially on circular orbits. While highly eccentric orbits can be produced naturally by these interactions, collisions between the two planets, which occur frequently in the range of observed semimajor axes, would result in many more nearly circular orbits than in the observed sample.

We have obtained spectral and timing data on GX 3+1 and Ser X-1. Both sources were observed simultaneously with BeppoSAX and RXTE. The RXTE data is used to provide power spectra and colour-colour diagrams in order to constrain the state (and thus track $\dot M$) the sources are in. The BeppoSAX data provide the broad-band spectra. The spectra of both sources are reasonably well-fit using a model consisting of a disk-blackbody, a comptonized component and a Fe line, absorbed by interstellar absorption. The electron temperature (kT$_{\rm e}$) of the Comptonizing plasma is in both cases $\sim$2.5 keV. This implies that no strong high-energy tail from the Comptonized component is present in either of the sources. We discuss the similarities between these burst sources and the luminous X-ray sources located in globular clusters. We find that the spectral parameters of the comptonized component provide information about the mass-accretion rate, which agrees well with estimates from the timing and spectral variations.

We study the timing properties of the bursting atoll source 4U 1728-34 as a function of its position in the X-ray color-color diagram. In the island part of the color-color diagram (corresponding to the hardest energy spectra) the power spectrum of 4U 1728-34 shows several features such as a band-limited noise component present up to a few tens of Hz, a low frequency quasi-periodic oscillation (LFQPO) at frequencies between 20 and 40 Hz, a peaked noise component around 100 Hz, and one or two QPOs at kHz frequencies. In addition to these, in the lower banana (corresponding to softer energy spectra) we also find a very low frequency noise (VLFN) component below ~1 Hz. In the upper banana (corresponding to the softest energy spectra) the power spectra are dominated by the VLFN, with a peaked noise component around 20 Hz. We find that the frequencies of the kHz QPOs are well correlated with the position in the X-ray color-color diagram. For the frequency of the LFQPO and the break frequency of the broad-band noise component the relation appears more complex. These frequencies both increase when the frequency of the upper kHz QPO increases from 400 to 900 Hz, but at this frequency a jump in the values of the parameters occurs. We interpret this jump in terms of the gradual appearance of a QPO at the position of the break at high inferred mass accretion rate, while the previous LFQPO disappears. Simultaneously, another kind of noise appears with a break frequency of ~7 Hz, similar to the NBO of Z sources. The 100 Hz peaked noise does not seem to correlate with the position of the source in the color-color diagram, but remains relatively constant in frequency. This component may be similar to several 100 Hz QPOs observed in black hole binaries.

We have used Rossi X-ray Timing Explorer data to measure the lags between soft (2-5 keV) and hard (5-13 keV) photons and to study the aperiodic variability of the superluminal black hole candidate GRS 1915+105 during low-flux states. The power density spectra exhibit quasi-periodic oscillations (QPO) whose frequency increases with increasing count rate and varies in the frequency range 0.6-8 Hz. A correlation between the QPO frequency and the phase lag spectra is reported for the first time. This correlation is found for both the phase lag continuum and the phase lag at the QPO frequency. We find that as the QPO frequency moves to higher values the phase lags reverse sign from positive to negative. The absolute value of the lag always increases with photon energy. The negative (soft) lags are associated with a softer energy spectrum, whereas the positive (hard) lags are seen when the source is harder. We describe a possible scenario that may account for the change in the sign of the lags.

We measure simultaneously the properties of the energy spectra and the frequencies of the kilohertz quasi-periodic oscillations (QPOs) in fifteen low mass X-ray binaries covering a wide range of X-ray luminosities. In each source the QPO frequencies cover the same range of approximately 300 Hz to 1300 Hz, though the sources differ by two orders of magnitude in their X-ray luminosities (as measured from the unabsorbed 2-50 keV flux). So the X-ray luminosity does not uniquely determine the QPO frequency. This is difficult to understand since the evidence from individual sources indicates that the frequency and luminosity are very well correlated at least over short timescales. Perhaps beaming effects or bolometric corrections change the observed luminosities, or perhaps part of the energy in mass accretion is used to power outflows reducing the energy emitted in X-rays. It is also possible that the parameters of a QPO model are tuned in such a way that the same range of frequencies appears in all sources. Different modes of accretion may be involved for example (disk and radial) or multiple parameters may conspire to yield the same frequencies.

We analyze the properties of the pulsed emission from the accreting millisecond pulsar SAX J1808.4-3658 in observations of its April 1998 outburst by the Rossi X-Ray Timing Explorer. Pulse phase spectroscopy shows that the emission evolves from a hard spectrum (power law with photon index 2.39+-0.06) to a soft spectrum (index 3.39+-0.24). This softening is also observable as a phase lag in the fundamental of low-energy photons with respect to high-energy photons. We show that this lag is roughly constant over ten days of the outburst. We fit these data with a model where the pulse emission is from a hot spot on the rotating neutron star and the flux as a function of phase is calculated including the effects of general relativity. The energy-dependent lags are very well described by this model. The harder spectra at earlier phases (as the spot approaches) are the result of larger Doppler boosting factors which are important for this fast pulsar.

We present an analysis of ~390 ksec of data of the Z source GX 340+0 taken during 24 observations with the Rossi X-ray Timing Explorer satellite. We report the discovery of a new broad component in the power spectra. The frequency of this component varied between 9 and 14 Hz, and remained close to half that of the horizontal branch quasi-periodic oscillations (HBO). Its rms amplitude was consistent with being constant around ~5%, while its FWHM increased with frequency from 7 to 18 Hz. If this sub-HBO component is the fundamental frequency, then the HBO and its second harmonic are the second and fourth harmonic component, while the third harmonic was not detected. This is similar to what was recently found for the black hole candidate XTE J1550-564. The profiles of both the horizontal- and the normal branch quasi-periodic oscillation peaks were asymmetric when they were strongest. We describe this in terms of a shoulder component at the high frequency side of the quasi-periodic oscillation peak, whose rms amplitudes were approximately constant at ~4% and \~3%, respectively. The peak separation between the twin kHz quasi-periodic oscillations was consistent with being constant at 339+-8 Hz but a trend similar to that seen in, e.g. Sco X-1 could not be excluded. We discuss our results within the framework of the various models which have been proposed for the kHz QPOs and low frequency peaks.

We study the correlations between timing and X-ray spectral properties in the low mass X-ray binary 4U 0614+09 using a large (265-ks) data set obtained with the Rossi X-ray Timing Explorer. We find strong quasi-periodic oscillations (QPOs) of the X-ray flux, like the kilohertz QPOs in many other X-ray binaries with accreting neutron stars, with frequencies ranging from 1329 Hz down to 418 Hz and, perhaps, as low as 153 Hz. We report the highest frequency QPO yet from any low mass X-ray binary at 1329+-4 Hz, which has implications for neutron star structure. This QPO has a 3.5-sigma single-trial significance, for an estimated 40 trials the significance is 2.4-sigma. Besides the kilohertz QPOs, the Fourier power spectra show four additional components: high frequency noise (HFN), described by a broken power-law with a break frequency between 0.7 and 45 Hz, very low frequency noise (VLFN), which is fitted as a power-law below 1 Hz, and two broad Lorentzians with centroid frequencies varying from 6 to 38 Hz and 97 to 158 Hz, respectively. We find strong correlations between the frequencies of the kilohertz QPOs, the frequency of the 6 to 38 Hz broad Lorentzian, the break frequency of the HFN, the strength of both the HFN and the VLFN and the position of the source in the hard X-ray color vs. intensity diagram. The frequency of the 97 to 158 Hz Lorentzian does not correlate with these parameters. We also find that the relation between power density and break frequency of the HFN is similar to that established for black hole candidates in the low state. We suggest that the changing mass accretion rate is responsible for the correlated changes in all these parameters.

We derive octupole-level secular perturbation equations for hierarchical triple systems, using classical Hamiltonian perturbation techniques. By extending previous work done to leading (quadrupole) order to octupole level (i.e., including terms of order alpha^3, where alpha = a_1/a_2 < 1 is the ratio of semimajor axes) we obtain expressions that are applicable to a much wider range of parameters. For triple systems containing a close inner binary, we also discuss the possible interaction between the classical Newtonian perturbations and the general relativistic precession of the inner orbit. In some cases we show that this interaction can lead to resonances and a significant increase in the maximum amplitude of eccentricity perturbations. We establish the validity of our analytic expressions by providing detailed comparisons with the results of direct numerical integrations of the three-body problem obtained for a large number of representative cases. We also discuss applications of the theory in the context of several observed triple systems of current interest, including the millisecond pulsar PSR B1620-26 in M4, the giant planet in 16 Cygni, and the protostellar binary TMR-1. [short version]

We present a new theoretical analysis of the PSR B1620-26 triple system in the globular cluster M4, based on the latest radio pulsar timing data, which now include measurements of five time derivatives of the pulse frequency. These data allow us to determine the mass and orbital parameters of the second companion completely (up to the usual unknown orbital inclination angle $i_2$). The current best-fit parameters correspond to a second companion of planetary mass, $m_2 \sin i_2 \simeq 7\times10^{-3} M_\odot$, in an orbit of eccentricity $e_2\simeq 0.45$ and semimajor axis $a_2\simeq 60 $AU. Using numerical scattering experiments, we study a possible formation scenario for the triple system, which involves a dynamical exchange interaction between the binary pulsar and a primordial star-planet system. The current orbital parameters of the triple are consistent with such a dynamical origin, and suggest that the separation of the parent star-planet system was very large, $\go 50 $AU. We also examine the possible origin of the anomalously high eccentricity of the inner binary pulsar. While this eccentricity could have been induced during the same dynamical interaction that created the triple, we find that it could equally well arise from long-term secular perturbation effects in the triple, combining the general relativistic precession of the inner orbit with the Newtonian gravitational perturbation of the planet. The detection of a planet in this system may be taken as evidence that large numbers of extrasolar planetary systems, not unlike those discovered recently in the solar neighborhood, also exist in old star clusters.

We have discovered a ~1 Hz quasi-periodic oscillation (QPO) in the persistent X-ray emission and during type I X-ray bursts of the globular cluster source, dipper and low-mass X-ray binary (LMXB) 4U 1746-37. The QPO properties resemble those of QPOs found recently in the LMXB dippers 4U 1323-62, and EXO 0748-676, which makes 4U 1746-37 the third source known to exhibit this type of QPOs. We present evidence for X-ray spectral changes in this source similar to those observed in LMXBs referred to as atoll sources. We detect two states, a low intensity and spectrally hard state, and a higher intensity and spectrally soft state. This may explain the different spectral characteristics reported for 4U 1746-37 earlier. The high intensity state resembles the banana branch state of atoll sources. The QPOs are only seen in the low intensity state, and are absent when the source is in the banana branch. This strongly suggests that either the accretion disk or an extended central source change shape between the low intensity state and the banana branch. Twelve bursts were detected, of which 5 took place while the source was on the banana branch and 7 when the source was in the low intensity state. The bursts occurring on the banana branch had an e-folding time ~3 times longer than those which occurred in the low intensity state. Whereas previously detected dips showed only a decrease in count rate of ~15%, we found in one observation a dip in which the count rate dropped from ~200 counts per second to ~20 counts per second. This dip lasted only ~250 seconds, during which clear spectral hardening occured. This is the first time strong evidence for spectral changes during a dip are reported for this source.

Based on Rossi X-ray Timing Explorer data, we discuss the classification of the soft X-ray transient Aquila X-1 in the Z/atoll scheme, and the relation of its kilohertz quasi-periodic oscillations (kHz QPO) properties to the X-ray colors. The color-color diagram shows one elongated ("banana") structure and several "islands" of data points. The power spectra of the island are best represented by a broken power-law, whereas those of the banana by a power-law below ~ 1 Hz plus an exponentially cut-off component at intermediate frequencies (30-60 Hz). The parameters of these two components change in correlation with the position of the source in the color-color diagram. Based on the pattern that the source shows in the color-color diagram and its aperiodic variability we conclude that Aquila X-1 is an atoll source. We have also investigated the possible correlation between the frequency of the kHz QPO and the position of the source in the color-color diagram. The complexity seen in the frequency versus count rate diagram is reduced to a single track when the frequency is plotted against hard or soft color.

We report the first simultaneous measurement of the broad band X-ray (0.3-150 keV) spectrum of the neutron star x-ray binary 4U0614+091. Our data confirm the presence of a hard x-ray tail that can be modeled as thermal Comptonization of low-energy photons on electrons having a very high temperature, greater than 220 keV, or as a non-thermal powerlaw. We detected a spectral feature that can be interpreted as reprocessing, via Compton reflection, of the direct emission by an optically-thick disk and found a correlation between the photon index of the power-law tail and the fraction of radiation reflected which is similar to the correlation found for black hole candidate x-ray binaries and Seyfert galaxies.

The behavior of quasi-periodic oscillations (QPOs) at frequencies near 1 kHz in the x-ray emission from the neutron star x-ray binary 4U1820-30 has been interpreted as evidence for the existence of the marginally stable orbit, a key prediction of strong-field general relativity. The signature of the marginally stable orbit is a saturation in QPO frequency, assumed to track inner disk radius, versus mass accretion rate. Previous studies of 4U1820-30 have used x-ray count rate as an indicator of mass accretion rate. However, x-ray count rate is known to not correlate robustly with mass accretion rate or QPO frequency in other sources. Here, we examine the QPO frequency dependence on two other indicators of mass accretion rate: energy flux and x-ray spectral shape. Using either of these indicators, we find that the QPO frequency saturates at high mass accretion rates. We interpret this as strong evidence for the existence of the marginally stable orbit.

Numerous X-ray bursts show strong oscillations in their flux at several hundred Hz as revealed by RXTE. Analyzing one such oscillation from the X-ray binary Aql X-1, I find that low energy photons (3.5-5.7 keV) lag high energy photons (>5.7 keV) by approximately 1 radian. The oscillations are thought to be produced by hot spots on the spinning neutron star. The lags can then be explained by a Doppler shifting of emission from the hot spots; higher energy photons being emitted earlier in the spin phase as the spot approaches the observer. A quantitative test of this simple model shows a remarkable agreement with the data. Similar low energy lags have been measured in kilohertz quasi-periodic oscillations and in the accreting millisecond pulsar SAX J1808.4-3658. A Doppler delay mechanism may be at work there as well.

Exploiting the presence of kilohertz quasi-periodic oscillations (QPOs) in the timing power spectrum, we find that the soft x-ray emission of the neutron-star X-ray binary and atoll source 4U1636-536 modulated at the QPO frequency lags behind that of the hard x-ray emission. Emission in the 3.8-6.4 keV band is delayed by 25.0 +/- 3.3 microseconds relative to the 9.3-69 keV band. The delay increases in magnitude with increasing energy. Our results are consistent with those of Vaughan et al. (1997), when the sign is corrected (Vaughan et al. 1998), for the atoll source 4U1608-52. The soft lag could be produced by Comptonization of hard photons injected into a cooler electron cloud or by intrinsic spectral softening of the emission during each oscillation cycle.

Using the Rossi X-ray Timing Explorer we have measured lags of the 9 to 33 keV photons relative to the 2 to 9 keV photons in the timing noise between 0.01 and 100 Hz in the accreting neutron stars 4U 0614+091 and 4U 1705-44. We performed similar measurements on the accreting black hole candidates Cyg X-1 and GX 339-4 as a comparison. During the observations these sources were all in low (hard) states. We find phase lags of between 0.03 and 0.2 radians in all these sources, with a variation in frequency much less than expected for a lag constant in time. We also measure a coherence consistent with unity in all sources. As already noted for the black hole candidates, these data are inconsistent with simple Comptonization models invoking a constant time delay. Comptonization in a non-uniform medium can perhaps explain the lags. However, the magnitudes of the lags imply that the hot electron gas extends to more than 1000 Schwarzschild radii. This may constitute an energy problem. We argue that while a large hot cloud is possible for black holes which may hide some of their accretion energy in advection, such a distribution may not be possible neutron stars where all the accretion energy is eventually released at the neutron star surface. This casts doubt on the Comptonization model, though the energy problem may be resolved for example by a wind from the inner disk.

We present new results based on observations carried out with the Rossi X-ray Timing Explorer during the decay of an outburst of the low-mass X-ray binary (LMXB) and atoll source 4U 1608-52. Our results appear to resolve, at least in 4U 1608-52, one of the long-standing issues about the phenomenology of the kilohertz quasi-periodic oscillations (kHz QPOs), namely, the lack of a unique relation between the frequency of the kHz QPOs and the X-ray flux. We show that despite its complex dependence on the X-ray flux, the frequency of the kHz QPOs is monotonically related to the position of the source in the color-color diagram. Our findings strengthen the idea that, as in the case of Z sources, in the atoll sources the X-ray flux is not a good indicator of $\dot M$, and that the observed changes in the frequency of the kHz QPOs in LMXBs are driven by changes in $\dot M$. These results raise some concern about the recently reported detection of the orbital frequency at the innermost stable orbit in 4U 1820-30.

Using the Yale stellar evolution code, we have calculated theoretical models for nearby stars with planetary-mass companions in short-period nearly circular orbits: 51 Pegasi, Tau Bootis, Upsilon Andromedae, Rho Cancri, and Rho Coronae Borealis. We present tables listing key stellar parameters such as mass, radius, age, and size of the convective envelope as a function of the observable parameters (luminosity, effective temperature, and metallicity), as well as the unknown helium fraction. For each star we construct best models based on recently published spectroscopic data and the present understanding of galactic chemical evolution. We discuss our results in the context of planet formation theory, and, in particular, tidal dissipation effects and stellar metallicity enhancements.

In recent observations with the Rossi X-Ray Timing Explorer we have detected two simultaneous quasi-periodic oscillation (QPO) peaks in the low mass X-ray binary and atoll source 4U 1735-44. The lower and higher frequency QPOs have frequencies varying between 632 and 729 Hz, and 982 and 1026 Hz, respectively. The fractional rms amplitudes are 3.7 to 8.1% and 5.0 to 5.8%. The frequency separation between the two QPOs changes from 341+-7 Hz to 296+-12 Hz. The inferred mass accretion rate during our observations is relatively low compared to that during the previous observations, where only a single QPO was present. There is weak evidence that the frequency of the QPOs correlates with the mass accretion rate, as observed in other binaries. Five X-ray bursts were recorded with no detectable oscillations with upper limits for the rms fraction of 4% to 13%.

We study the timing properties of the low mass X-ray binary 4U 1728-34 using recently released data from the Rossi X-Ray Timing Explorer. This binary, like many others with accreting neutron stars, is known to exhibit strong quasi-periodic oscillations (QPOs) of its X-ray flux near 1 kHz. In addition to the kilohertz QPOs, the Fourier power spectra show a broken power law noise component, with a break frequency between 1 and 50 Hz, and a Lorentzian between 10 and 50 Hz. We find that the frequencies of the break and the low-frequency Lorentzian are well correlated with the frequencies of the kilohertz QPOs. The slope of the correlation is similar to that expected if the oscillations are due to relativistic frame dragging (Lense-Thirring precession) in the inner accretion disk (Stella & Vietri 1998). The correlation is also nearly identical to the one found in Z-sources between the the well known QPOs on the horizontal branch and the kilohertz QPOs, suggesting that the low frequency oscillations are a similar phenomenon in these sources. The frequency of the break in the power spectra is also correlated with the frequencies of the kilohertz QPOs. As previously noted for the similar binaries 4U 1608-50 and 4U 1705-44, this broken power law component closely resembles that of black hole candidates in the low state, where the break frequency is taken as an indicator of mass accretion rate. The relation between break frequency and kilohertz QPO frequency thus provides additional proof that the frequency of the kilohertz QPOs increases with mass accretion rate.

We present new Rossi X-ray Timing Explorer observations of the low-mass X-ray binary 4U 1608-52 during the decay of its 1998 outburst. We detect by a direct FFT method the existence of a second kilohertz quasi-periodic oscillation (kHz QPO) in its power density spectrum, previously only seen by means of the sensitivity-enhancing `shift and add' technique. This result confirms that 4U 1608-52 is a twin kHz QPO source. The frequency separation between these two QPO decreased significantly, from 325.5 +/- 3.4 Hz to 225.3 +/- 12.0 Hz, as the frequency of the lower kHz QPO increased from 470 Hz to 865 Hz, in contradiction with a simple beat-frequency interpretation. This change in the peak separation of the kHz QPOs is closely similar to that previously seen in Sco X-1, but takes place at a ten times lower average luminosity. We discuss this result within the framework of models that have been proposed for kHz QPO. Beat frequency models where the peak separation is identified with the neutron star spin rate, as well as the explanations previously proposed to account for the similar behavior of the QPOs in Sco X-1, are strongly challenged by this result.

In observations with the Rossi X-ray Timing Explorer we have discovered quasi-periodic oscillations (QPOs) near 1 kHz from 4U 1705-44, a low--mass X-ray binary with a neutron star classified as an atoll source. In six separate observations, we detect one QPO with a frequency ranging between 770 and 870 Hz and a 4% rms fraction in the full detector energy band. There is evidence for a second QPO at 1073 Hz in one interval. The separation in frequency of the two QPOs is 298+-11 Hz. The QPOs are present only in observations where the mass accretion rate is inferred to be at an intermediate level, based on the atoll source phenomenology. At the highest accretion rates, the QPOs are not detected with upper limits to the rms fraction of about 2%. At the lowest accretion rates the upper limits are about 4%. The QPO frequency increases with inferred mass accretion rate. This is expected in models where the QPO frequency is generated by motion at an inner edge of the accretion disk. An increased mass accretion rate causes the disk edge to move in, increasing the orbital frequency. Five Type-I X-ray bursts are observed with no detectable oscillations.

We identify the periodic transient hard X-ray source GRO J1849-03 with the transient x-ray pulsar GS 1843-02 = X1845-024 based on the detection of x-ray outbursts from X1845-024 coincident with hard x-ray outbursts of GRO J1849--03. Based on its spin period of 94.8 s and its orbital period of 241 days, we classify the system as a Be/X-ray binary.

We present the results from RXTE/PCA observations of 4U 1608-52 in its island state on March 15, 18 and 22 of 1996. Three type I X-ray bursts were detected in one RXTE orbit on March 22. We observed QPO features peaking at 567-800 Hz on March 15 and 22, with source fractional rms amplitude of 13%-17% and widths of 78-180 Hz in the power density spectra averaged over each spacecraft orbit. The rms amplitudes of these QPOs are positively correlated with the photon energy. The three X-ray bursts, with burst intervals of 16 and 8 minutes, have a duration of 16s. The blackbody emission region of the smallest X-ray burst among the three suggest it was a local nuclear burning. We also discuss a type I X-ray burst candidate in the observation.

We report RXTE/PCA observations of 4U 1608-52 on March 15, 18 and 22 immediately after the outburst in early 1996. The persistent count rates ranged from 190 to 450 cps (1-60 keV). During this period of time, 4U 1608-52 was in the island state. We detected QPO features in the power density spectra (PDS) at 567-800 Hz on March 15 and 22, with source fractional root-mean-square (rms) amplitude of 13%-17% and widths of 78-180 Hz. The average rms amplitude of these QPO features is positively correlated with the energy. Our results imply that the neutron star spin frequency is possibly between 300 Hz and 365 Hz.

We investigate the behavior of the high frequency quasi-periodic oscillations (QPOs) in 4U 0614+091, combining timing and spectral analysis of RXTE observations. The energy spectra of the source can be described by a power law (alpha ~ 2.8) and a blackbody (kT ~ 1.5), with the blackbody accounting for 10 - 20% of the total energy flux. We find a robust correlation of the frequency, nu, of the higher frequency QPO near 1 kHz with the flux of the blackbody, F_BB. The slope of this correlation, dlog(nu)/dlog(F_BB), is 0.27 to 0.37. The source follows the same relation even in observations separated by several months. The QPO frequency does not have a similarly unique correlation with the total flux or the flux of the power law component. The RMS fraction of the higher frequency QPO rises with energy from 6.8 +- 1.5 % (3-5 keV) to 21.3 +- 4.0 % (10-12 keV). For the lower frequency QPO, however, it is consistent with a constant value of 5.4 +- 0.9 %. The results may be interpreted in terms of a beat frequency model for the production of the high frequency QPOs.

Quasi-periodic oscillations (QPOs) at frequencies near 1000 Hz were recently discovered in several x-ray binaries containing neutron stars. Two sources show no correlation between QPO frequency and source count rate (Berger et al. 1996, Zhang et al. 1996). We suggest that the QPO frequency is determined by the Keplerian orbital frequency near the marginally stable orbit predicted by general relativity in strong gravitational fields (Muchotrzeb-Czerny 1986, Paczynski 1987, Kluzniak et al. 1990). The QPO frequencies observed from 4U 1636-536 imply that the mass of the neutron star is 2.02 +/- 0.12 solar masses. Interpretation of the 4.1 keV absorption line observed from 4U 1636-536 (Waki et al. 1984) as due to Fe XXV ions then implies a neutron star radius of 9.6 +/-0.6 km.

We report the results of a CGRO 3-week observation of the binary system containing the 47 ms pulsar PSR B1259-63 orbiting around a Be star companion in a very eccentric orbit. The PSR B1259-63 binary is a unique system for the study of the interaction of a rapidly rotating pulsar with time-variable nebular surroundings. CGRO observed the PSR B1259-63 system in coincidence with its most recent periastron passage (January 3-23, 1994). Unpulsed and non-thermal hard X-ray emission was detected up to 200 keV, with a photon index $1.8 \pm 0.2$ and a flux of ~4 mCrab, corresponding to a luminosity of a few 10^{34} erg/s at the distance of 2 kpc. The hard X-ray flux and spectrum detected by CGRO agrees with the X-ray emission measured by simultaneous ASCA observations. EGRET upper limits are significant, and exclude strong inverse Compton cooling in the PSR B1259-63 system. We interpret the observed non-thermal emission as synchrotron radiation of shocked electron/positron pairs of the relativistic pulsar wind interacting with the mass outflow from the Be star. Our results clearly indicate, for the first time in a binary pulsar, that high energy emission can be shock-powered rather than caused by accretion. The lack of X-ray/gamma-ray pulsations constrains models of high-energy emission from rapidly rotating pulsars.

We have detected quasi-periodic oscillations (QPOs) near 1 kHz from the low mass X-ray binary 4U 0614+091 in observations with XTE. The observations span several months and sample the source over a large range of X-ray luminosity. In every interval QPOs are present above 400 Hz with fractional RMS amplitudes from 3 to 12%. At high count rates, two high frequency QPOs are detected simultaneously. The difference of their frequency centroids is consistent with 323 Hz in all observations. During one interval a third signal is detected at 328 +/- 2 Hz. This suggests the system has a stable `clock' which is most likely the neutron star with spin period 3.1 msec. Thus, our observations and those of another neutron star system by Strohmayer et al. (1996) provide the first evidence for millisecond pulsars within low-mass X-ray binary systems and reveal the `missing-link' between millisecond radiopulsars and the late stages of binary evolution in low mass X-ray binaries (Alpar et al. 1982). We suggest that the kinematics of the magnetospheric beat-frequency model (Alpar and Shaham 1985) applies to these QPOs. In this interpretation the high frequency signal is associated with the Keplerian frequency of the inner accretion disk and the lower frequency `beat' signal arises from the differential rotation frequency of the inner disk and the spinning neutron star. Assuming the high frequency QPO is a Keplerian orbital frequency for the accretion disk, we find a maximum mass of 1.9 solar masses and a maximum radius of 17 km for the neutron star.

We have detected transient X-ray activity from the X-ray burster 4U~0614+091 simultaneously with BATSE/CGRO (20-100 keV) and ASM/RXTE (1-12 keV). The peak fluxes reach approximately 40 mCrab in both instruments over a period of about 20 days. The variable emission shows a clear anticorrelation of the hard X-ray flux with the soft X-ray count rate. The observed anticorrelation is another clear counterexample to the notion that only black hole binaries exhibit such correlations. The individual spectra during this period can be fit by power laws with photon indices 2.2+-0.3 (ASM) and 2.7+-0.4 (BATSE), while the combined spectra can be described by a single power law with index 2.09+-0.08. BATSE and the ASM/RXTE are a good combination for monitoring X-ray sources over a wide energy band.

We have monitored a sample of 27 nearby globular clusters in the hard X-ray band (20-120 keV) for approximately 1400 days using the BATSE instrument on board the Compton Gamma-Ray Observatory. Globular clusters may contain a large number of compact objects (e.g., pulsars or X-ray binaries containing neutron stars) which can produce hard X-ray emission. Our search provides a sensitive (~50 mCrab) monitor for hard X-ray transient events on time scales of >1 day and a means for observing persistent hard X-ray emission. We have discovered no transient events from any of the clusters and no persistent emission. Our observations include a sensitive search of four nearby clusters containing dim X-ray sources: 47 Tucanae, NGC 5139, NGC 6397, and NGC 6752. The non-detection in these clusters implies a lower limit for the recurrence time of transients of 2 to 6 years for events with luminosities >10^36 erg s^-1 (20-120 keV) and ~20 years if the sources in these clusters are taken collectively. This suggests that the dim X-ray sources in these clusters are not transients similar to Aql~X-1. We also place upper limits on the persistent emission in the range 2-10*10^34 erg s^-1 (2 sigma, 20-120 keV) for these four clusters. For 47 Tuc the upper limit is more sensitive than previous measurements by a factor of 3. We find a model dependent upper limit of 19 isolated millisecond pulsars (MSPs) producing gamma-rays in 47 Tuc, compared to the 11 observed radio MSPs in this cluster.