The Kepler mission revealed a class of planets known as ''super-puffs,'' with masses only a few times larger than Earth's but radii larger than Neptune, giving them very low mean densities. All three of the known planets orbiting the young solar-type star Kepler 51 are super-puffs. The Kepler 51 system thereby provides an opportunity for a comparative study of the structures and atmospheres of this mysterious class of planets, which may provide clues about their formation and evolution. We observed two transits each of Kepler 51b and 51d with the Wide Field Camera 3 (WFC3) on the Hubble Space Telescope. Combining new WFC3 transit times with re-analyzed Kepler data and updated stellar parameters, we confirmed that all three planets have densities lower than 0.1 g/cm$^{3}$. We measured the WFC3 transmission spectra to be featureless between 1.15 and 1.63 $\mu$m, ruling out any variations greater than 0.6 scale heights (assuming a H/He dominated atmosphere), thus showing no significant water absorption features. We interpreted the flat spectra as the result of a high-altitude aerosol layer (pressure $<$3 mbar) on each planet. Adding this new result to the collection of flat spectra that have been observed for other sub-Neptune planets, we find support for one of the two hypotheses introduced by Crossfield and Kreidberg (2017), that planets with cooler equilibrium temperatures have more high-altitude aerosols. We strongly disfavor their other hypothesis that the H/He mass fraction drives the appearance of large amplitude transmission features.

While planets between the size of Uranus and Saturn are absent within the Solar System, the star K2-24 hosts two such planets, K2-24b and c, with radii equal to $5.4~R_E$ and $7.5~R_E$, respectively. The two planets have orbital periods of 20.9 days and 42.4 days, residing only 1% outside the nominal 2:1 mean-motion resonance. In this work, we present results from a coordinated observing campaign to measure planet masses and eccentricities that combines radial velocity (RV) measurements from Keck/HIRES and transit-timing measurements from K2 and Spitzer. K2-24b and c have low, but non-zero, eccentricities of $e_1 \sim e_2 \sim 0.08$. The low observed eccentricities provide clues regarding the formation and dynamical evolution of K2-24b and K2-24c, suggesting that they could be the result of stochastic gravitational interactions with a turbulent protoplanetary disk, among other mechanisms. K2-24b and c are $19\pm2~M_E$ and $15\pm2~M_E$, respectively; K2-24c is 20% less massive than K2-24b, despite being 40% larger. Their large sizes and low masses imply large envelope fractions, which we estimate at $26^{+3}_{-3}\%$ and $52^{+5}_{-3}\%$. In particular, K2-24c's large envelope presents an intriguing challenge to the standard model of core nucleated accretion that predicts the onset of runaway accretion when $f_{env} \approx 50\%$.

Measuring precise planet masses, densities, and orbital dynamics in individual planetary systems is an important pathway toward understanding planet formation. The WASP-47 system has an unusual architecture that motivates a complex formation theory. The system includes a hot Jupiter ("b") neighbored by interior ("e") and exterior ("d") sub-Neptunes, and a long-period eccentric giant planet ("c"). We simultaneously modeled transit times from the Kepler K2 Mission and 118 radial velocities to determine precise masses, densities, and Keplerian orbital elements of the WASP-47 planets. Combining RVs and TTVs provides a better estimate of the mass of planet d (13.6\pm2.0~M_\odot) than obtained with only RVs ($12.75\pm2.70~M_\odot$) or TTVs ($16.1\pm3.8~M_\odot$). Planets e and d have high densities for their size, consistent with a history of photo-evaporation and/or formation in a volatile-poor environment. Through our RV and TTV analysis, we find that the planetary orbits have eccentricities similar to the solar system planets. The WASP-47 system has three similarities to our own solar system: (1) the planetary orbits are nearly circular and coplanar, (2) the planets are not trapped in mean motion resonances, and (3) the planets have diverse compositions. None of the current single-process exoplanet formation theories adequately reproduce these three characteristics of the WASP-47 system (or our solar system). We propose that WASP-47, like the solar system, formed in two stages: first, the giant planets formed in a gas-rich disk and migrated to their present locations, and second, the high-density sub-Neptunes formed in situ in a gas-poor environment.

One focus of modern astronomy is to detect temperate terrestrial exoplanets well-suited for atmospheric characterisation. A milestone was recently achieved with the detection of three Earth-sized planets transiting (i.e. passing in front of) a star just 8% the mass of the Sun 12 parsecs away. Indeed, the transiting configuration of these planets with the Jupiter-like size of their host star - named TRAPPIST-1 - makes possible in-depth studies of their atmospheric properties with current and future astronomical facilities. Here we report the results of an intensive photometric monitoring campaign of that star from the ground and with the Spitzer Space Telescope. Our observations reveal that at least seven planets with sizes and masses similar to the Earth revolve around TRAPPIST-1. The six inner planets form a near-resonant chain such that their orbital periods (1.51, 2.42, 4.04, 6.06, 9.21, 12.35 days) are near ratios of small integers. This architecture suggests that the planets formed farther from the star and migrated inward. The seven planets have equilibrium temperatures low enough to make possible liquid water on their surfaces.

We report precise radial velocity (RV) measurements of WASP-47, a G star that hosts three transiting planets in close proximity (a hot Jupiter, a super-Earth and a Neptune-sized planet) and a non-transiting planet at 1.4 AU. Through a joint analysis of previously published RVs and our own Keck-HIRES RVs, we significantly improve the planet mass and bulk density measurements. For the super-Earth WASP-47e ($P$ = 0.79 days), we measure a mass of 9.11 $\pm$ 1.17 $M_\oplus$, and a bulk density of 7.63 $\pm$ 1.90 g cm$^{-3}$, consistent with a rocky composition. For the hot Jupiter WASP-47b ($P$ = 4.2 days), we measure a mass of 356 $\pm$ 12 $M_\oplus$ (1.12 $\pm$ 0.04 $M_\rm{Jup}$) and constrain its eccentricity to $<0.021$ at 3-$\sigma$ confidence. For the Neptune-size planet WASP-47d ($P$ = 9.0 days), we measure a mass of 12.75 $\pm$ 2.70 $M_\oplus$, and a bulk density of 1.36 $\pm$ 0.42 g cm$^{-3}$, suggesting it has a thick H/He envelope. For the outer non-transiting planet, we measure a minimum mass of 411 $\pm$ 18 $M_\oplus$ (1.29 $\pm$ 0.06 $M_\rm{Jup}$), an orbital period of 595.7 $\pm$ 5.0 days, and an orbital eccentricity of 0.27 $\pm$ 0.04. Our new measurements are consistent with but 2$-$4$\times$ more precise than previous mass measurements.

We report the discovery and confirmation of two sub-Saturn planets orbiting a bright (V = 11.3), metal-rich ([Fe/H] = 0.42 $\pm$ 0.04 dex) G3 dwarf in the K2 Campaign 2 field. The planets are 5.68 $\pm$ 0.56 Earth-radii and 7.82 $\pm$ 0.72 Earth-radii and have orbital periods of 20.8851 $\pm$ 0.0003 d and 42.3633$\pm$0.0006 d, near to the 2:1 mean-motion resonance. We obtained 32 radial velocities (RVs) with Keck/HIRES and detected the reflex motion due to EPIC-203771098b and c. These planets have masses of 21.0 $\pm$ 5.4 Earth-masses and 27.0 $\pm$ 6.9 Earth-masses, respectively. With low densities of 0.63 $\pm$ 0.25 g/cc and 0.31 $\pm$ 0.12 g/cc, respectively, the planets require thick envelopes of H/He to explain their large sizes and low masses. Interior structure models predict that the planets have fairly massive cores of 17.6 $\pm$ 4.3 Earth-masses and 16.1 $\pm$ 4.2 Earth-masses, respectively. They may have formed exterior to their present locations, accreted their H/He envelopes at large orbital distances, and migrated in as a resonant pair. The proximity to resonance, large transit depths, and host star brightness offer rich opportunities for TTV follow-up. Finally, the low surface gravities of the EPIC-203771098 planets make them favorable targets for transmission spectroscopy by HST, Spitzer, and JWST.

We derive the transit timing variations (TTVs) of two planets near a second order mean motion resonance on nearly circular orbits. We show that the TTVs of each planet are given by sinusoids with a frequency of $j n_2-(j-2)n_1$, where $j \ge 3$ is an integer characterizing the resonance and $n_2$ and $n_1$ are the mean motions of the outer and inner planets, respectively. The amplitude of the TTV depends on the mass of the perturbing planet, relative to the mass of the star, and on both the eccentricities and longitudes of pericenter of each planet. The TTVs of the two planets are approximated anti-correlated, with phases of $\phi$ and $\approx \phi+\pi$, where the phase $\phi$ also depends on the eccentricities and longitudes of pericenter. Therefore, the TTVs caused by proximity to a second order mean motion resonance do not in general uniquely determine both planet masses, eccentricities, and pericenters. This is completely analogous to the case of TTVs induced by two planets near a first order mean motion resonance. We explore how other TTV signals, such as the short-period synodic TTV or a first order resonant TTV, in combination with the second order resonant TTV, can break degeneracies. Lastly, we derive approximate formulae for the TTVs of planets near any order eccentricity-type mean motion resonance; this shows that the same basic sinusoidal TTV structure holds for all eccentricity-type resonances. Our general formula reduces to previously derived results near first order mean motion resonances.

Characterization of transiting planets with transit timing variations (TTVs) requires understanding how to translate the observed TTVs into masses and orbital elements of the planets. This can be challenging in multi-planet transiting systems, but fortunately these systems tend to be nearly plane-parallel and low eccentricity. Here we present a novel derivation of analytic formulae for TTVs that are accurate to first order in the planet-star mass ratios and in the orbital eccentricities. These formulae are accurate in proximity to first order resonances, as well as away from resonance, and compare well with more computationally expensive N-body integrations in the low eccentricity, low mass-ratio regime when applied to simulated and to actual multi-transiting Kepler planet systems. We make code available for implementing these formulae.

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

We have used the {\it Spitzer Space Telescope} to observe two transiting planetary systems orbiting low mass stars discovered in the \Kepler \Ktwo mission. The system K2-3 (EPIC 201367065) hosts three planets while EPIC 202083828 (K2-26) hosts a single planet. Observations of all four objects in these two systems confirm and refine the orbital and physical parameters of the planets. The refined orbital information and more precise planet radii possible with \Spitzer will be critical for future observations of these and other \Ktwo targets. For K2-3b we find marginally significant evidence for a Transit Timing Variation between the \Ktwo and \Spitzer\ epochs.

We combine new high resolution imaging and spectroscopy from Keck/NIRC2, Discovery Channel Telescope/DSSI, and Keck/HIRES with published astrometry and radial velocities to measure individual masses and orbital elements of the GJ 3305 AB system, a young (~20 Myr) M+M binary (unresolved spectral type M0) member of the beta Pictoris moving group comoving with the imaged exoplanet host 51 Eri. We measure a total system mass of 1.11 \pm 0.04 M_sun, a period of 29.03 \pm 0.50$ yr, a semimajor axis of 9.78 \pm 0.14 AU, and an eccentricity of 0.19 \pm 0.02. The primary component has a dynamical mass of 0.67 \pm 0.05 M_sun and the secondary has a mass of 0.44 \pm 0.05 M_sun. The recently updated BHAC15 models are consistent with the masses of both stars to within 1.5 sigma. Given the observed masses the models predict an age of the GJ 3305 AB system of 37 \pm 9 Myr. Based on the the observed system architecture and our dynamical mass measurement, it is unlikely that the orbit of 51 Eri b has been significantly altered by the Kozai-Lidov mechanism.

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

Discoveries from the prime Kepler mission demonstrated that small planets (< 3 Earth-radii) are common outcomes of planet formation. While Kepler detected many such planets, all but a handful orbit faint, distant stars and are not amenable to precise follow up measurements. Here, we report the discovery of two small planets transiting K2-21, a bright (K = 9.4) M0 dwarf located 65$\pm$6 pc from Earth. We detected the transiting planets in photometry collected during Campaign 3 of NASA's K2 mission. Analysis of transit light curves reveals that the planets have small radii compared to their host star, 2.60 $\pm$ 0.14% and 3.15 $\pm$ 0.20%, respectively. We obtained follow up NIR spectroscopy of K2-21 to constrain host star properties, which imply planet sizes of 1.59 $\pm$ 0.43 Earth-radii and 1.92 $\pm$ 0.53 Earth-radii, respectively, straddling the boundary between high-density, rocky planets and low-density planets with thick gaseous envelopes. The planets have orbital periods of 9.32414 days and 15.50120 days, respectively, and have a period ratio of 1.6624, very near to the 5:3 mean motion resonance, which may be a record of the system's formation history. Transit timing variations (TTVs) due to gravitational interactions between the planets may be detectable using ground-based telescopes. Finally, this system offers a convenient laboratory for studying the bulk composition and atmospheric properties of small planets with low equilibrium temperatures.

We consider the dynamical evolution of two planets orbiting in the vicinity of a first order mean motion reso- nance while simultaneously undergoing eccentricity damping and convergent migration. Following Goldreich & Schlichting (2014), we include a coupling between the dissipative semimajor axis evolution and the damping of the eccentricities. In agreement with past studies, we find that this coupling can lead to overstability of the resonance and that for a certain range of parameters capture into resonance is only temporary. Using a more general model, we show that whether overstable motion can occur depends in a characteristic way on the mass ratio between the two planets as well as their relative eccentricity damping timescales. Moreover, we show that even when escape from resonance does occur, the timescale for escape is long enough such at any given time a pair of planets is more likely to be found in a resonance rather than migrating between them. Thus, we argue that overstability of resonances cannot singlehandedly reconcile convergent migration with the observed lack of Kepler planet pairs found near resonances. However, it is possible that overstable motion in combination with other effects such as large scale orbital instability could produce the observed period ratio distribution.

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.

Gravitational interactions between planets in transiting exoplanetary systems lead to variations in the times of transit that are diagnostic of the planetary masses and the dynamical state of the system. Here we show that synodic "chopping" contributions to these transit timing variations (TTVs) can be used to uniquely measure the masses of planets without full dynamical analyses involving direct integration of the equations of motion. We present simple analytic formulae for the chopping signal, which are valid (generally <10% error) for modest eccentricities e <~ 0.1. Importantly, these formulae primarily depend on the mass of the perturbing planet, and therefore the chopping signal can be used to break the mass/free-eccentricity degeneracy which can appear for systems near first order mean motion resonances. Using a harmonic analysis, we apply these TTV formulae to a number of Kepler systems which had been previously analyzed with full dynamical analyses. We show that when chopping is measured, the masses of both planets can be determined uniquely, in agreement with previous results, but without the need for numerical orbit integrations. This demonstrates how mass measurements from TTVs may primarily arise from an observable chopping signal. The formula for chopping can also be used to predict the number of transits and timing precision required for future observations, such as those made by TESS or PLATO, in order to infer planetary masses through analysis of TTVs.

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

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.

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

We present the results of a search through the photometric database of eclipsing Kepler binaries (Prsa et al. 2011; Slawson et al. 2011) looking for evidence of hierarchical triple star systems. The presence of a third star orbiting the binary can be inferred from eclipse timing variations. We apply a simple algorithm in an automated determination of the eclipse times for all 2157 binaries. The "calculated" eclipse times, based on a constant period model, are subtracted from those observed. The resulting O-C (observed minus calculated times) curves are then visually inspected for periodicities in order to find triple-star candidates. After eliminating false positives due to the beat frequency between the ~1/2-hour Kepler cadence and the binary period, 39 candidate triple systems were identified. The periodic O-C curves for these candidates were then fit for contributions from both the classical Roemer delay and so-called "physical" delay, in an attempt to extract a number of the system parameters of the triple. We discuss the limitations of the information that can be inferred from these O-C curves without further supplemental input, e.g., ground-based spectroscopy. Based on the limited range of orbital periods for the triple star systems to which this search is sensitive, we can extrapolate to estimate that at least 20% of all close binaries have tertiary companions.

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

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

The Kepler-11 planetary system contains six transiting planets ranging in size from 1.8 to 4.2 times the radius of Earth. Five of these planets orbit in a tightly-packed configuration with periods between 10 and 47 days. We perform a dynamical analysis of the system based upon transit timing variations observed in more than three years of \ik photometric data. Stellar parameters are derived using a combination of spectral classification and constraints on the star's density derived from transit profiles together with planetary eccentricity vectors provided by our dynamical study. Combining masses of the planets relative to the star from our dynamical study and radii of the planets relative to the star from transit depths together with deduced stellar properties yields measurements of the radii of all six planets, masses of the five inner planets, and an upper bound to the mass of the outermost planet, whose orbital period is 118 days. We find mass-radius combinations for all six planets that imply that substantial fractions of their volumes are occupied by constituents that are less dense than rock. The Kepler-11 system contains the lowest mass exoplanets for which both mass and radius have been measured.

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

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