We combine analytical understanding of resonant dynamics in two-planet systems with machine learning techniques to train a model capable of robustly classifying stability in compact multi-planet systems over long timescales of $10^9$ orbits. Our Stability of Planetary Orbital Configurations Klassifier (SPOCK) predicts stability using physically motivated summary statistics measured in integrations of the first $10^4$ orbits, thus achieving speed-ups of up to $10^5$ over full simulations. This computationally opens up the stability constrained characterization of multi-planet systems. Our model, trained on $\approx 100,000$ three-planet systems sampled at discrete resonances, generalizes both to a sample spanning a continuous period-ratio range, as well as to a large five-planet sample with qualitatively different configurations to our training dataset. Our approach significantly outperforms previous methods based on systems' angular momentum deficit, chaos indicators, and parametrized fits to numerical integrations. We use SPOCK to constrain the free eccentricities between the inner and outer pairs of planets in the Kepler-431 system of three approximately Earth-sized planets to both be below 0.05. Our stability analysis provides significantly stronger eccentricity constraints than currently achievable through either radial velocity or transit duration measurements for small planets, and within a factor of a few of systems that exhibit transit timing variations (TTVs). Given that current exoplanet detection strategies now rarely allow for strong TTV constraints (Hadden et al., 2019), SPOCK enables a powerful complementary method for precisely characterizing compact multi-planet systems. We publicly release SPOCK for community use.

We report the discovery and validation of four extrasolar planets hosted by the nearby, bright, Sun-like (G3V) star HD~108236 using data from the Transiting Exoplanet Survey Satellite (TESS). We present transit photometry, reconnaissance and precise Doppler spectroscopy as well as high-resolution imaging, to validate the planetary nature of the objects transiting HD~108236, also known as the TESS Object of Interest (TOI) 1233. The innermost planet is a possibly-rocky super-Earth with a period of $3.79523_{-0.00044}^{+0.00047}$ days and has a radius of $1.586\pm0.098$ $R_\oplus$. The outer planets are sub-Neptunes, with potential gaseous envelopes, having radii of $2.068_{-0.091}^{+0.10}$ $R_\oplus$, $2.72\pm0.11$ $R_\oplus$, and $3.12_{-0.12}^{+0.13}$ $R_\oplus$ and periods of $6.20370_{-0.00052}^{+0.00064}$ days, $14.17555_{-0.0011}^{+0.00099}$ days, and $19.5917_{-0.0020}^{+0.0022}$ days, respectively. With V and K$_{\rm s}$ magnitudes of 9.2 and 7.6, respectively, the bright host star makes the transiting planets favorable targets for mass measurements and, potentially, for atmospheric characterization via transmission spectroscopy. HD~108236 is the brightest Sun-like star in the visual (V) band known to host four or more transiting exoplanets. The discovered planets span a broad range of planetary radii and equilibrium temperatures, and share a common history of insolation from a Sun-like star ($R_\star = 0.888 \pm 0.017$ R$_\odot$, $T_{\rm eff} = 5730 \pm 50$ K), making HD 108236 an exciting, opportune cosmic laboratory for testing models of planet formation and evolution.

We present ground-based infrared transit observations for four dynamically interacting \textit{Kepler} planets, including Kepler-29b, Kepler-36c, KOI-1783.01, and Kepler-177c, obtained using the Wide-field Infrared Camera on the Hale 200" telescope at Palomar Observatory. By utilizing an engineered diffuser and custom guiding software, we mitigate time-correlated telluric and instrumental noise sources in these observations. We achieve an infrared photometric precision comparable to or better than that of space-based observatories such as the \textit{Spitzer Space Telescope}, and detect transits with greater than 3$\sigma$ significance for all planets. For Kepler-177c ($J=13.9$) our measurement uncertainties are only $1.2\times$ the photon noise limit and 1.9 times better than the predicted photometric precision for \textit{Spitzer} IRAC photometry of this same target. We find that a single transit observation obtained $4-5$ years after the end of the original \textit{Kepler} mission can reduce dynamical mass uncertainties by as much as a factor of three for these systems. Additionally, we combine our new observations of KOI-1783.01 with information from the literature to confirm the planetary nature of this system. We discuss the implications of our new mass and radius constraints in the context of known exoplanets with low incident fluxes, and we note that Kepler-177c may be a more massive analog to the currently known super-puffs given its core mass (3.8$\pm0.9M_\Earth$) and large gas-to-core ratio (2.8$\pm0.7$). Our demonstrated infrared photometric performance opens up new avenues for ground-based observations of transiting exoplanets previously thought to be restricted to space-based investigation.

Low-mass planets have an extraordinarily diverse range of bulk compositions, from primarily rocky worlds to those with deep gaseous atmospheres. As techniques for measuring the masses of exoplanets are advancing the field towards the terrestrial regime, from ultra-short orbital periods to Venus-like distances, we identify the bounds on planet compositions, where sizes and incident fluxes inform bulk planet properties. In some cases, measurement precisions of planet masses and sizes are approaching the theoretical uncertainties in planet models. An emerging picture explains aspects of the diversity of low-mass planets although some problems remain; do extreme low density low-mass planets challenge models of atmospheric mass loss? Are planet sizes strictly separated by bulk composition? Why do some stellar characterizations differ between observational techniques? As the TESS mission begins, low-mass exoplanets around the nearest stars will soon be discovered and characterized with unprecedented precision, permitting more detailed planetary modeling and atmospheric characterization of low mass exoplanets than ever before.

K2-19b and c were among the first planets discovered by NASA's K2 mission and together stand in stark contrast with the physical and orbital properties of the solar system planets. The planets are between the size of Uranus and Saturn at 7.0$\pm$0.2 R_E and 4.1$\pm$0.2 R_E, respectively, and reside a mere 0.1% outside the nominal 3:2 mean-motion resonance. They represent a different outcome of the planet formation process than the solar system, as well as the vast majority of known exoplanets. We measured the physical and orbital properties of these planets using photometry from K2, Spitzer, and ground-based telescopes, along with radial velocities from Keck/HIRES. Through a joint photodynamical model, we found that the planets have moderate eccentricities of $e \approx0.20$ and well-aligned apsides $\Delta \varpi \approx 0$ deg. The planets occupy a strictly non-resonant configuration: the resonant angles circulate rather than librate. This defies the predictions of standard formation pathways that invoke convergent or divergent migration, both of which predict $\Delta \varpi \approx 180$ deg and eccentricities of a few percent or less. We measured masses of $M_{p,b}$ = 32.4$\pm$1.7 M_E and $M_{p,c}$ = 10.8$\pm$0.6 M_E. Our measurements, with 5% fractional uncertainties, are among the most precise of any sub-Jovian exoplanet. Mass and size reflect a planet's core/envelope structure. Despite having a relatively massive core of $M_{core} \approx15$ $M_E$, K2-19b is envelope-rich, with an envelope mass fraction of roughly 50%. This planet poses a challenge to standard models core-nucleated accretion, which predict that cores $\gtrsim 10$ $M_E$ will quickly accrete gas and trigger runaway accretion when the envelope mass exceeds that of the core.

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

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

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

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

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

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

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

We present 50 nights of polarimetric observations of HD 189733 in $B$ band using the POLISH2 aperture-integrated polarimeter at the Lick Observatory Shane 3-m telescope. This instrument, commissioned in 2011, is designed to search for Rayleigh scattering from short-period exoplanets due to the polarized nature of scattered light. Since these planets are spatially unresolvable from their host stars, the relative contribution of the planet-to-total system polarization is expected to vary with an amplitude of order 10 parts per million (ppm) over the course of the orbit. Non-zero and also variable at the 10 ppm level, the inherent polarization of the Lick 3-m telescope limits the accuracy of our measurements and currently inhibits conclusive detection of scattered light from this exoplanet. However, the amplitude of observed variability conservatively sets a $3 \sigma$ upper limit to the planet-induced polarization of the system of 58 ppm in $B$ band, which is consistent with a previous upper limit from the POLISH instrument at the Palomar Observatory 5-m telescope (Wiktorowicz 2009). A physically-motivated Rayleigh scattering model, which includes the depolarizing effects of multiple scattering, is used to conservatively set a $3 \sigma$ upper limit to the geometric albedo of HD 189733b of $A_g < 0.37$. This value is consistent with the value $A_g = 0.226 \pm 0.091$ derived from occultation observations with HST STIS (Evans et al. 2013), but it is inconsistent with the large $A_g = 0.61 \pm 0.12$ albedo reported by (Berdyugina et al. 2011).

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

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

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

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

KOI-152 is among the first known systems of multiple transiting planetary candidates (Steffen et al. 2010) ranging in size from 3.5 to 7 times the size of the Earth, in a compact configuration with orbital periods near a 1:2:4:6 chain of commensurability, from 13.5 days to 81.1 days. All four planets exhibit transit timing variations with periods that are consistent with the distance of each planet to resonance with its neighbors. We perform a dynamical analysis of the system based on transit timing measurements over 1282 days of \textit{Kepler} photometry. Stellar parameters are obtained with a combination of spectral classification and the stellar density constraints provided by light curve analysis and orbital eccentricity solutions from our dynamical study. Our models provide tight constraints on the masses of all four transiting bodies, demonstrating that they are planets and that they orbit the same star. All four of KOI-152's transiting planets have low densities given their sizes, consistent with other studies of compact multiplanet transiting systems. The largest of the four, KOI-152.01, has the lowest bulk density yet determined amongst sub-Saturn mass planets.

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.

We study the radial and vertical stability of dust grains launched with all charge-to-mass ratios at arbitrary distances from rotating planets with complex magnetic fields. We show that the aligned dipole magnetic field model analyzed by Jontof-Hutter and Hamilton (2012) is an excellent approximation in most cases, but that fundamentally new physics arises with the inclusion of non-axisymmetric magnetic field terms. In particular, large numbers of distant negatively-charged dust grains, stable in a magnetic dipole, can be driven to escape by a more complex field. We trace the origin of the instability to overlapping Lorentz resonances which are extremely powerful when the gravitational and electromagnetic forces on a dust grain are comparable. These resonances enable a dust grain to tap the spin energy of the planet to power its escape. We also explore the relatively minor influence of different launch speeds and the far more important effects of variable grain charge. Only the latter are capable of significantly affecting the micron-sized grains that dominate visible and infrared images of faint dust rings. Finally, we present full stability maps for Earth, Jupiter, Saturn, Uranus, and Neptune with magnetic fields modeled out to octupole order. Not surprisingly, dust in the tortured magnetic fields of Uranus and Neptune show the greatest instability.

We study the stability of charged dust grains orbiting a planet and subject to gravity and the electromagnetic force. Our numerical models cover a broad range of launch distances from the planetary surface to beyond synchronous orbit, and the full range of charge-to-mass ratios from ions to rocks. Treating the spinning planetary magnetic field as an aligned dipole, we map regions of radial and vertical instability where dust grains are driven to escape or crash into the planet. We derive the boundaries between stable and unstable trajectories analytically, and apply our models to Jupiter, Saturn and the Earth, whose magnetic fields are reasonably well represented by aligned dipoles.

Using data mostly assembled by previous authors, we consider the linear correlation between the apparent radiative efficiency $\epsilon_{\gamma}$ (defined as the ratio of isotropic equivalent radiative output to inferred isotropic equivalent kinetic energy of the blast) and $E_{peak}^{\alpha}$ where $1.4<\alpha<2$, for 17 of 22 GRBs (Lloyd-Ronning and Zhang, 2004). We note in a quantitative manner that this is consistent with the hypothesis that $\epsilon_{\gamma}$ and $E_{peak}$ are influenced by viewing angle. We suggest a more general theoretically derived expression for this correlation that could be tested with a richer data set. If the reduction in both $\epsilon_{\gamma}$ and $E_{peak}$ is due to viewing angle effects, then the actual radiative efficiency is $\sim 7$. We also find preliminary evidence (with a small sample) for a separate class of weak GRB afterglows.