Source author record

Matthew J. Holman

Matthew J. Holman appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

95works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

95 published item(s)

preprint2026arXiv

You Only Stack Once (YOSO): A Motion-Filtered, Deep-Learning Framework for Detecting Faint Moving Sources

We present You Only Stack Once (YOSO), an automated pipeline designed to detect faint, slow-moving Solar System objects in wide-field astronomical surveys. The pipeline integrates a novel Gaussian Motion Filter (GMoF) that operates at the pixel level to enhance signal-to-noise for objects exhibiting a range of apparent rates of motion. Unlike conventional shift-and-stack methods, which rely on discrete velocity trials, GMoF amplifies trails while suppressing random noise and static background features. Applied to a subset of DEEP observations from the Dark Energy Camera, YOSO recovered 45 out of 73 previously detected objects, as well as 11 new TNOs. It also discovered 216 objects in the near Solar System. Although alternative shift-and-stack methods are sensitive to objects about 0.88 magnitudes fainter, YOSO's false positive rate is extremely low, since it detects only sources that exhibit a trail and are consistent with a point source when shifted at the right rate. We show how this method can be deployed on large surveys like LSST, and adapted for other domains that require motion-based signal enhancement, including exoplanet imaging through Angular Differential Imaging (ADI), and near-Earth object (NEO) detection for missions like NEO Surveyor. YOSO thus provides a versatile, scalable approach for extracting faint, motion-dependent signals in the era of data-intensive astronomy.

preprint2022arXiv

From Data to Software to Science with the Rubin Observatory LSST

The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) dataset will dramatically alter our understanding of the Universe, from the origins of the Solar System to the nature of dark matter and dark energy. Much of this research will depend on the existence of robust, tested, and scalable algorithms, software, and services. Identifying and developing such tools ahead of time has the potential to significantly accelerate the delivery of early science from LSST. Developing these collaboratively, and making them broadly available, can enable more inclusive and equitable collaboration on LSST science. To facilitate such opportunities, a community workshop entitled "From Data to Software to Science with the Rubin Observatory LSST" was organized by the LSST Interdisciplinary Network for Collaboration and Computing (LINCC) and partners, and held at the Flatiron Institute in New York, March 28-30th 2022. The workshop included over 50 in-person attendees invited from over 300 applications. It identified seven key software areas of need: (i) scalable cross-matching and distributed joining of catalogs, (ii) robust photometric redshift determination, (iii) software for determination of selection functions, (iv) frameworks for scalable time-series analyses, (v) services for image access and reprocessing at scale, (vi) object image access (cutouts) and analysis at scale, and (vii) scalable job execution systems. This white paper summarizes the discussions of this workshop. It considers the motivating science use cases, identified cross-cutting algorithms, software, and services, their high-level technical specifications, and the principles of inclusive collaborations needed to develop them. We provide it as a useful roadmap of needs, as well as to spur action and collaboration between groups and individuals looking to develop reusable software for early LSST science.

preprint2022arXiv

The Effects of Disk Induced Apsidal Precession on Planets Captured into Mean Motion Resonance

The process of migration into resonance capture has been well studied for planetary systems where the gravitational potential is generated exclusively by the star and planets. However, massive protoplanetary disks add a significant perturbation to these models. In this paper we consider two limiting cases of disk-induced precession on migrating planets and find that small amounts of precession significantly affect the equilibrium reached by migrating planets. We investigate these effects with a combination of semi-analytic models of the resonance and numerical integrations. We also consider the case of the disk's dispersal, which can excite significant libration amplitude and can cause ejection from resonance for large enough precession rates. Both of these effects have implications for interpreting the known exoplanet population and may prove to be important considerations as the population of well-characterized exoplanet systems continues to grow.

preprint2020arXiv

EnckeHH: an integrator for gravitational dynamics with a dominant mass that achieves optimal error behaviour

We present EnckeHH, a new, highly accurate code for orbital dynamics of perturbed Keplerian systems such as planetary systems or galactic centre systems. It solves Encke's equations of motion, which assume perturbed Keplerian orbits. By incorporating numerical techniques, we have made the code follow optimal roundoff error growth. In a $10^ {12}$ day integration of the outer Solar System, EnckeHH was $3.5$ orders of magnitude more accurate than IAS15 in a fixed time step test. Adaptive steps are recommended for IAS15. Through study of efficiency plots, we show that EnckeHH reaches significantly higher accuracy than the Rebound integrators IAS15 and WHCKL for fixed step size.

preprint2020arXiv

The continuing search for evidence of tidal orbital decay of hot Jupiters

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

preprint2020arXiv

The Occurrence of Rocky Habitable Zone Planets Around Solar-Like Stars from Kepler Data

We present occurrence rates for rocky planets in the habitable zones (HZ) of main-sequence dwarf stars based on the Kepler DR25 planet candidate catalog and Gaia-based stellar properties. We provide the first analysis in terms of star-dependent instellation flux, which allows us to track HZ planets. We define $η_\oplus$ as the HZ occurrence of planets with radius between 0.5 and 1.5 $R_\oplus$ orbiting stars with effective temperatures between 4800 K and 6300 K. We find that $η_\oplus$ for the conservative HZ is between $0.37^{+0.48}_{-0.21}$ (errors reflect 68\% credible intervals) and $0.60^{+0.90}_{-0.36}$ planets per star, while the optimistic HZ occurrence is between $0.58^{+0.73}_{-0.33}$ and $0.88^{+1.28}_{-0.51}$ planets per star. These bounds reflect two extreme assumptions about the extrapolation of completeness beyond orbital periods where DR25 completeness data are available. The large uncertainties are due to the small number of detected small HZ planets. We find similar occurrence rates using both a Poisson likelihood Bayesian analysis and Approximate Bayesian Computation. Our results are corrected for catalog completeness and reliability. Both completeness and the planet occurrence rate are dependent on stellar effective temperature. We also present occurrence rates for various stellar populations and planet size ranges. We estimate with $95\%$ confidence that, on average, the nearest HZ planet around G and K dwarfs is about 6 pc away, and there are about 4 HZ rocky planets around G and K dwarfs within 10 pc of the Sun.

preprint2019arXiv

The Orbit of WASP-12b is Decaying

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

preprint2016arXiv

A Dynamical Analysis of the Kepler-80 System of Five Transiting Planets

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.

preprint2016arXiv

Discovery of A New Retrograde Trans-Neptunian Object: Hint of A Common Orbital Plane for Low Semi-Major Axis, High Inclination TNOs and Centaurs

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

preprint2016arXiv

Liberating exomoons in white dwarf planetary systems

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

preprint2016arXiv

Observational Constraints on Planet Nine: Astrometry of Pluto and Other Trans-Neptunian Objects

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

preprint2016arXiv

The Fate of Exomoons in White Dwarf Planetary Systems

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

preprint2016arXiv

The Pan-STARRS 1 Discoveries of five new Neptune Trojans

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

preprint2016arXiv

Uncovering Circumbinary Planetary Architectural Properties from Selection Biases

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

preprint2015arXiv

Dynamical Evolution of Multi-Resonant Systems: the Case of GJ876

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

preprint2015arXiv

Tests of the planetary hypothesis for PTFO 8-8695b

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

preprint2014arXiv

Architecture of Kepler's Multi-transiting Systems: II. New investigations with twice as many candidates

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.

preprint2014arXiv

Masses, Radii, and Orbits of Small Kepler Planets: The Transition from Gaseous to Rocky Planets

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

preprint2014arXiv

Planet Hunters VII. Discovery of a New Low-Mass, Low-Density Planet (PH3 c) Orbiting Kepler-289 with Mass Measurements of Two Additional Planets (PH3 b and d)

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

preprint2014arXiv

spotrod: a semi-analytic model for transits of spotted stars

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

preprint2014arXiv

Stellar rotation--planetary orbit period commensurability in the HAT-P-11 system

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

preprint2014arXiv

The Transiting Exoplanet Survey Satellite

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

preprint2014arXiv

TTVFast: An efficient and accurate code for transit timing inversion problems

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

preprint2013arXiv

2011 HM102: Discovery of a High-Inclination L5 Neptune Trojan in the Search for a post-Pluto New Horizons Target

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

preprint2013arXiv

First order resonance overlap and the stability of close two planet systems

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

preprint2013arXiv

Relative photometry of HAT-P-1b occultations

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

preprint2013arXiv

Stability of Satellites in Closely Packed Planetary Systems

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

preprint2012arXiv

Alignment of the stellar spin with the orbits of a three-planet system

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.

preprint2012arXiv

Almost All of Kepler's Multiple Planet Candidates are Planets

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.

preprint2012arXiv

CfA4: Light Curves for 94 Type Ia Supernovae

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

preprint2012arXiv

Improved Asteroid Astrometry and Photometry with Trail Fitting

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

preprint2012arXiv

Improved spectroscopic parameters for transiting planet hosts

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

preprint2012arXiv

Kepler constraints on planets near hot Jupiters

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.

preprint2012arXiv

Kepler-20: A Sun-like Star with Three Sub-Neptune Exoplanets and Two Earth-size Candidates

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

preprint2012arXiv

Kepler-36: A Pair of Planets with Neighboring Orbits and Dissimilar Densities

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.

preprint2012arXiv

Kepler-47: A Transiting Circumbinary Multi-Planet System

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.

preprint2012arXiv

Long-Term Stability of Horseshoe Orbits

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

preprint2012arXiv

Long-Term Transit Timing Monitoring and Refined Light Curve Parameters of HAT-P-13b

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

preprint2012arXiv

Planetary Candidates Observed by Kepler, III: Analysis of the First 16 Months of Data

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.

preprint2012arXiv

Rapid dynamical chaos in an exoplanetary system

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

preprint2012arXiv

The Neptune-Sized Circumbinary Planet Kepler-38b

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.

preprint2012arXiv

The Transiting Circumbinary Planets Kepler-34 and Kepler-35

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.

preprint2012arXiv

Transit Timing Observations from Kepler: II. Confirmation of Two Multiplanet Systems via a Non-parametric Correlation Analysis

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.

preprint2012arXiv

Transit Timing Observations from Kepler: III. Confirmation of 4 Multiple Planet Systems by a Fourier-Domain Study of Anti-correlated Transit Timing Variations

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.

preprint2012arXiv

Transit Timing Observations from Kepler: IV. Confirmation of 4 Multiple Planet Systems by Simple Physical Models

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.

preprint2012arXiv

Transit Timing Observations from Kepler: VI. Potentially interesting candidate systems from Fourier-based statistical tests

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.

preprint2012arXiv

Transit Timing Observations from Kepler: VI. Transit Timing Variation Candidates in the First Seventeen Months from Polynomial Models

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.

preprint2011arXiv

A Closely-Packed System of Low-Mass, Low-Density Planets Transiting Kepler-11

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.

preprint2011arXiv

A First Comparison of Kepler Planet Candidates in Single and Multiple Systems

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.

preprint2011arXiv

A Search for Additional Planets in Five of the Exoplanetary Systems Studied by the NASA EPOXI Mission

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

preprint2011arXiv

A Super-Earth Transiting a Naked-Eye Star

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

preprint2011arXiv

Architecture and Dynamics of Kepler's Candidate Multiple Transiting Planet 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.

preprint2011arXiv

Characteristics of planetary candidates observed by Kepler, II: Analysis of the first four months of data

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.

preprint2011arXiv

Dynamically Excited Outer Solar System Objects in the Hubble Space Telescope Archive

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

preprint2011arXiv

Kepler 18-b, c, and d: A System Of Three Planets Confirmed by Transit Timing Variations, Lightcurve Validation, Spitzer Photometry and Radial Velocity Measurements

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.

preprint2011arXiv

Kepler-10c, a 2.2-Earth radius transiting planet in a multiple system

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

preprint2011arXiv

Kepler-14b: A massive hot Jupiter transiting an F star in a close visual binary

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.

preprint2011arXiv

Kepler-16: A Transiting Circumbinary Planet

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.

preprint2011arXiv

Kepler-22b: A 2.4 Earth-radius Planet in the Habitable Zone of a Sun-like Star

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

preprint2011arXiv

KOI-126: A Triply-Eclipsing Hierarchical Triple with Two Low-Mass Stars

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.

preprint2011arXiv

Planet Occurrence within 0.25 AU of Solar-type Stars from Kepler

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^α. 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]

preprint2011arXiv

Starspots and spin-orbit alignment in the WASP-4 exoplanetary system

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

preprint2011arXiv

The Distribution of Transit Durations for Kepler Planet Candidates and Implications for their Orbital Eccentricities

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.

preprint2011arXiv

The First Kepler Mission Planet Confirmed With The Hobby-Eberly Telescope: Kepler-15b, a Hot Jupiter Enriched In Heavy Elements

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.

preprint2011arXiv

The Kepler-19 System: A Transiting 2.2 R_Earth Planet and a Second Planet Detected via Transit Timing Variations

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 μ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]

preprint2011arXiv

The Transit Light Curve Project. XIII. Sixteen Transits of the Super-Earth GJ 1214b

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

preprint2011arXiv

The Transit Light Curve project. XIV. Confirmation of Anomalous Radii for the Exoplanets TrES-4b, HAT-P-3b, and WASP-12b

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

preprint2011arXiv

Transit Timing Observations from Kepler: I. Statistical Analysis of the First Four Months

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.

preprint2011arXiv

Two Earth-sized planets orbiting Kepler-20

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.

preprint2010arXiv

A Prograde, Low-Inclination Orbit for the Very Hot Jupiter WASP-3b

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

preprint2010arXiv

A Search for a Sub-Earth Sized Companion to GJ 436 and a Novel Method to Calibrate Warm Spitzer IRAC Observations

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

preprint2010arXiv

A Search for Additional Planets in the NASA EPOXI Observations of the Exoplanet System GJ 436

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

preprint2010arXiv

Five Kepler target stars that show multiple transiting exoplanet candidates

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.

preprint2010arXiv

Modeling Kepler transit light curves as false positives: Rejection of blend scenarios for Kepler-9, and validation of Kepler-9d, a super-Earth-size planet in a multiple system

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.

preprint2010arXiv

NLTT 41135: a field M-dwarf + brown dwarf eclipsing binary in a triple system, discovered by the MEarth observatory

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

preprint2010arXiv

Orbital Orientations of Exoplanets: HAT-P-4b is Prograde and HAT-P-14b is Retrograde

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

preprint2010arXiv

Precise Estimates of the Physical Parameters for the Exoplanet System HD-17156 Enabled by HST FGS Transit and Asteroseismic Observations

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

preprint2010arXiv

System parameters, transit times and secondary eclipse constraints of the exoplanet systems HAT-P-4, TrES-2, TrES-3 and WASP-3 from the NASA EPOXI Mission of Opportunity

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

preprint2010arXiv

The inclinations of faint TNOs

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

preprint2010arXiv

The Rossiter-McLaughlin Effect of the Transiting Exoplanet XO-4b

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

preprint2010arXiv

The Value of Systems with Multiple Transiting Planets

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

preprint2010arXiv

Trans-Neptunian Objects with Hubble Space Telescope ACS/WFC

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

preprint2009arXiv

A Search for Occultations of Bright Stars by Small Kuiper Belt Objects using Megacam on the MMT

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

preprint2009arXiv

A super-Earth transiting a nearby low-mass star

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

preprint2009arXiv

HAT-P-7: A Retrograde or Polar Orbit, and a Third Body

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

preprint2009arXiv

Near-infrared transit photometry of the exoplanet HD 149026b

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

preprint2009arXiv

Studying the atmosphere of the exoplanet HAT-P-7b via secondary eclipse measurements with EPOXI, Spitzer and Kepler

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

preprint2009arXiv

The Transit Light Curve Project. XII. Six Transits of the Exoplanet XO-2b

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

preprint2008arXiv

A Subaru Pencil-beam Search for m_R~27 Trans-neptunian bodies

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

preprint2008arXiv

Photometric Follow-up Observations of the Transiting Neptune-Mass Planet GJ 436b

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

preprint2008arXiv

The Transit Light Curve Project. X. A Christmas Transit of HD 17156b

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

preprint2006arXiv

Precise Radius Estimates for the Exoplanets WASP-1b and WASP-2b

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

preprint2005arXiv

Measurement of Spin-Orbit Alignment in an Extrasolar Planetary System

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

preprint2004arXiv

The Use of Transit Timing to Detect Extrasolar Planets with Masses as Small as Earth

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

preprint2003arXiv

Photometric Survey of the Irregular Satellites

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