Torques from a mutually inclined perturber can change a transiting planet's impact parameter, resulting in variations in the transit shape and duration. Detection of and upper limits on changes in impact parameter yield valuable constraints on a planetary system's three dimensional architecture. Constraints for warm Jupiters are particularly interesting because they allow us to test origins theories that invoke a mutually inclined perturber. Because of warm Jupiters' high signal-to-noise transits, changes in impact parameter are feasible to detect. However, here we show that allowing the impact parameter to vary uniformly and independently from transit to transit leads to incorrect inferences about the change, propagating to incorrect inferences about the perturber. We demonstrate that an appropriate prior on the change in impact parameter mitigates this problem. We apply our approach to eight systems from the literature and find evidence for changes in impact parameter for warm Jupiter Kepler-46b. We conclude with our recommendations for light curve fitting, including when to fit impact parameters vs. transit durations.

Super-Earths and mini-Neptunes exhibit great diversity in their compositional and orbital properties. Their bulk densities span a large range, from those dense enough to be purely rocky to those needing a substantial contribution from volatiles to their volumes. Their orbital configurations range from compact, circular multi-transiting systems like Kepler-11 to systems like our Solar System's terrestrial planets with wider spacings and modest but significant eccentricities and mutual inclinations. Here we investigate whether a continuum of formation conditions resulting from variation in the amount of solids available in the inner disk can account for the diversity of orbital and compositional properties observed for super Earths, including the apparent dichotomy between single transiting and multiple transiting system. We simulate in situ formation of super-Earths via giant impacts and compare to the observed Kepler sample. We find that intrinsic variations among disks in the amount of solids available for in situ formation can account for the orbital and compositional diversity observed among Kepler's transiting planets. Our simulations can account for the planets' distributions of orbital period ratios, transit duration ratios, and transit multiplicity; higher eccentricities for single than multi transiting planets; smaller eccentricities for larger planets; scatter in the mass-radius relation, including lower densities for planets with masses measured with TTVs than RVs; and similarity in planets' sizes and spacings within each system. Our findings support the theory that variation among super-Earth and mini-Neptune properties is primarily locked in by different in situ formation conditions, rather than arising stochastically through subsequent evolution.

Resolved debris disk features (e.g., warps, offsets, edges and gaps, azimuthal asymmetries, radially thickened rings, scale heights) contain valuable information about the underlying planetary systems, such as the posited planet's mass, semi-major axis, and other orbital parameters. Most existing models assume a single planet is sculpting the disk feature, but recent observations of mature planetary systems (e.g., by radial velocity surveys or \textit{Kepler}) have revealed that many planets reside in multi-planet systems. Here we investigate if/how planet properties inferred from single-planet models are compromised when multiple planets reside in the system. For each disk feature, we build a two-planet model that includes a planet b with fixed parameters and a planet c with a full range of possible parameters. We investigate these two-planet systems and summarize the configurations for which assuming a single planet (i.e., planet b) leads to significantly flawed inferences of that planet's properties. We find that although disk features are usually primarily dominated by a single planet, when using single-planet models we are at risk of misinterpreting planet properties by orders of magnitude in extreme cases. Specifically, we are at high risk of misinterpreting planet properties from disk warps; at moderate risk from disk edges and gaps, radially thickened rings, and scale height features; and at low risk from host star-disk center offsets and azimuthal asymmetries. We summarize situations where we can infer the need to use a multi-planet model instead of a single-planet one from disk morphology dissimilarities.

Resonant dynamics plays a significant role in the past evolution and current state of our outer Solar System. The population ratios and spatial distribution of Neptune's resonant populations are direct clues to understanding the history of our planetary system. The orbital structure of the objects in Neptune's 2:1 mean-motion resonance (\emph{twotinos}) has the potential to be a tracer of planetary migration processes. Different migration processes produce distinct architectures, recognizable by well-characterized surveys. However, previous characterized surveys only discovered a few twotinos, making it impossible to model the intrinsic twotino population. With a well-designed cadence and nearly 100\% tracking success, the Outer Solar System Origins Survey (OSSOS) discovered 838 trans-Neptunian objects, of which 34 are securely twotinos with well-constrained libration angles and amplitudes. We use the OSSOS twotinos and the survey characterization parameters via the OSSOS Survey Simulator to inspect the intrinsic population and orbital distributions of twotino. The estimated twotino population, 4400$^{+1500}_{-1100}$ with $H_r<8.66$ (diameter$\sim$100km) at 95\% confidence, is consistent with the previous low-precision estimate. We also constrain the width of the inclination distribution to a relatively narrow value of $\sigma_i$=6$^\circ$$^{+1}_{-1}$, and find the eccentricity distribution is consistent with a Gaussian centered on $e_\mathrm{c}=0.275$ with a width $e_\mathrm{w}=0.06$. We find a single-slope exponential luminosity function with $\alpha=0.6$ for the twotinos. Finally, we for the first time meaningfully constrain the fraction of symmetric twotinos, and the ratio of the leading asymmetric islands; both fractions are in a range of 0.2--0.6. These measurements rule out certain theoretical models of Neptune's migration history.

Warm, large exoplanets with 10-100 day orbital periods pose a major challenge to our understanding of how planetary systems form and evolve. Although high eccentricity tidal migration has been invoked to explain their proximity to their host stars, a handful reside in or near orbital resonance with nearby planets, suggesting a gentler history of in situ formation or disk migration. Here we confirm and characterize a pair of warm, large exoplanets discovered by the TESS Mission orbiting K-dwarf TOI-216. Our analysis includes additional transits and transit exclusion windows observed via ground-based follow-up. We find two families of solutions, one corresponding to a sub-Saturn-mass planet accompanied by a Neptune-mass planet and the other to a Jupiter in resonance with a sub-Saturn-mass planet. We prefer the second solution based on the orbital period ratio, the planet radii, the lower free eccentricities, and libration of the 2:1 resonant argument, but cannot rule out the first. The free eccentricities and mutual inclination are compatible with stirring by other, undetected planets in the system, particularly for the second solution. We discuss prospects for better constraints on the planets' properties and orbits through follow-up, including transits observed from the ground.

We present a statistical analysis of the first 300 stars observed by the Gemini Planet Imager Exoplanet Survey (GPIES). This subsample includes six detected planets and three brown dwarfs; from these detections and our contrast curves we infer the underlying distributions of substellar companions with respect to their mass, semi-major axis, and host stellar mass. We uncover a strong correlation between planet occurrence rate and host star mass, with stars M $>$ 1.5 $M_\odot$ more likely to host planets with masses between 2-13 M$_{\rm Jup}$ and semi-major axes of 3-100 au at 99.92% confidence. We fit a double power-law model in planet mass (m) and semi-major axis (a) for planet populations around high-mass stars (M $>$ 1.5M$_\odot$) of the form $\frac{d^2 N}{dm da} \propto m^\alpha a^\beta$, finding $\alpha$ = -2.4 $\pm$ 0.8 and $\beta$ = -2.0 $\pm$ 0.5, and an integrated occurrence rate of $9^{+5}_{-4}$% between 5-13 M$_{\rm Jup}$ and 10-100 au. A significantly lower occurrence rate is obtained for brown dwarfs around all stars, with 0.8$^{+0.8}_{-0.5}$% of stars hosting a brown dwarf companion between 13-80 M$_{\rm Jup}$ and 10-100 au. Brown dwarfs also appear to be distributed differently in mass and semi-major axis compared to giant planets; whereas giant planets follow a bottom-heavy mass distribution and favor smaller semi-major axes, brown dwarfs exhibit just the opposite behaviors. Comparing to studies of short-period giant planets from the RV method, our results are consistent with a peak in occurrence of giant planets between ~1-10 au. We discuss how these trends, including the preference of giant planets for high-mass host stars, point to formation of giant planets by core/pebble accretion, and formation of brown dwarfs by gravitational instability.

The recent discovery of a staggering diversity of planets beyond the Solar System has brought with it a greatly expanded search space for habitable worlds. The Kepler exoplanet survey has revealed that most planets in our interstellar neighborhood are larger than Earth and smaller than Neptune. Collectively termed super-Earths and mini-Neptunes, some of these planets may have the conditions to support liquid water oceans, and thus Earth-like biology, despite differing in many ways from our own planet. In addition to their quantitative abundance, super-Earths are relatively large and are thus more easily detected than true Earth twins. As a result, super-Earths represent a uniquely powerful opportunity to discover and explore a panoply of fascinating and potentially habitable planets in 2020 - 2030 and beyond.

We test the high-eccentricity tidal migration scenario for Kepler-419b, a member of the eccentric warm Jupiter class of planets whose origin is debated. Kepler-419 hosts two known planets (b,c). However, in its current configuration, planet c cannot excite the eccentricity of planet b enough to undergo high-eccentricity tidal migration. We investigate whether the presence of an undiscovered fourth body could explain the orbit of Kepler-419b. We explore the parameter space of this potential third giant planet using a suite of N-body simulations with a range of initial conditions. From the results of these simulations, coupled with observational constraints, we can rule out this mechanism for much of the parameter space of initial object d conditions. However, for a small range of parameters (masses between 0.5 and 7 $m_{\rm{Jup}}$, semi-major axes between 4 and 7.5 AU, eccentricities between 0.18 and 0.35, and mutual inclinations near 0$^{\circ}$) an undiscovered object d could periodically excite the eccentricity of Kepler-419b without destabilizing the system over 1 Gyr while producing currently undetectable radial velocity and transit timing variation signals.

We determine the orbital eccentricities of individual small Kepler planets, through a combination of asteroseismology and transit light-curve analysis. We are able to constrain the eccentricities of 51 systems with a single transiting planet, which supplement our previous measurements of 66 planets in multi-planet systems. Through a Bayesian hierarchical analysis, we find evidence that systems with only one detected transiting planet have a different eccentricity distribution than systems with multiple detected transiting planets. The eccentricity distribution of the single-transiting systems is well described by the positive half of a zero-mean Gaussian distribution with a dispersion $\sigma_e = 0.32 \pm 0.06$, while the multiple-transit systems are consistent with $\sigma_e = 0.083^{+0.015}_{-0.020}$. A mixture model suggests a fraction of $0.76^{+0.21}_{-0.12}$ of single-transiting systems have a moderate eccentricity, represented by a Rayleigh distribution that peaks at $0.26^{+0.04}_{-0.06}$. This finding may reflect differences in the formation pathways of systems with different numbers of transiting planets. We investigate the possibility that eccentricities are "self-excited" in closely packed planetary systems, as well as the influence of long-period giant companion planets. We find that both mechanisms can qualitatively explain the observations. We do not find any evidence for a correlation between eccentricity and stellar metallicity, as has been seen for giant planets. Neither do we find any evidence that orbital eccentricity is linked to the detection of a companion star. Along with this paper we make available all of the parameters and uncertainties in the eccentricity distributions, as well as the properties of individual systems, for use in future studies.

A question driving many studies is whether the thousands of exoplanets known today typically formed where we observe them or formed further out in the disk and migrated in. Early discoveries of giant exoplanets orbiting near their host stars and exoplanets in or near mean motion resonances were interpreted as evidence for migration and its crucial role in the beginnings of planetary systems. long-scale migration has been invoked to explain systems of planets in mean motion resonant chains consisting of three or more planets linked by integer period ratios. However, recent studies have reproduced specific resonant chains in systems via short-scale migration, and eccentricity damping has been shown to capture planets into resonant chains. We investigate whether the observed resonant chains in Kepler-80, Kepler-223, Kepler-60, and TRAPPIST-1 can be established through long-scale migration, short-scale migration, and/or only eccentricity damping by running suites of N-body simulations. We find that, for each system, all three mechanisms are able to reproduce the observed resonant chains. long-scale migration is not the only plausible explanation for resonant chains in these systems, and resonant chains are potentially compatible with in situ formation.

The HR 8799 system uniquely harbors four young super-Jupiters whose orbits can provide insights into the system's dynamical history and constrain the masses of the planets themselves. Using the Gemini Planet Imager (GPI), we obtained down to one milliarcsecond precision on the astrometry of these planets. We assessed four-planet orbit models with different levels of constraints and found that assuming the planets are near 1:2:4:8 period commensurabilities, or are coplanar, does not worsen the fit. We added the prior that the planets must have been stable for the age of the system (40 Myr) by running orbit configurations from our posteriors through $N$-body simulations and varying the masses of the planets. We found that only assuming the planets are both coplanar and near 1:2:4:8 period commensurabilities produces dynamically stable orbits in large quantities. Our posterior of stable coplanar orbits tightly constrains the planets' orbits, and we discuss implications for the outermost planet b shaping the debris disk. A four-planet resonance lock is not necessary for stability up to now. However, planet pairs d and e, and c and d, are each likely locked in two-body resonances for stability if their component masses are above $6~M_{\rm{Jup}}$ and $7~M_{\rm{Jup}}$, respectively. Combining the dynamical and luminosity constraints on the masses using hot-start evolutionary models and a system age of $42 \pm 5$~Myr, we found the mass of planet b to be $5.8 \pm 0.5~M_{\rm{Jup}}$, and the masses of planets c, d, and e to be $7.2_{-0.7}^{+0.6}~M_{\rm{Jup}}$ each.

Today, Pluto and Eris are the largest and most massive Trans-Neptunian Objects respectively. They are believed to be the last remnants of a population of planetesimals that has been reduced by >99% since the time of its formation. This reduction implies a primordial population of hundreds or thousands of Pluto-mass objects, and a mass-number distribution that could have extended to hundreds of Lunas, dozens of Mars, and several Earths. Such lost protoplanets would have left signatures in the dynamics of the present-day Trans-Neptunian Populations, and we statistically limit their primordial number by considering the survival of ultra-wide binary TNOs, the Cold Classical Kuiper belt, and the resonant population. We find that if the primordial mass-number distribution extended to masses greater than Pluto (~1e-3 Earth masses), it must have turned downwards to be no more top-heavy than roughly equal mass per log size, a significant deviation from the distribution observed between 1e-5 and 1e-3 Earth masses. We compare these limits to the predicted mass-number distribution of various planetesimal and proto-planet growth models. The limits derived here provide a test for future models of planetesimal formation.

The Outer Solar System Origins Survey (OSSOS), a wide-field imaging program in 2013-2017 with the Canada-France-Hawaii Telescope, surveyed 155 deg$^{2}$ of sky to depths of $m_r = 24.1$-25.2. We present 838 outer Solar System discoveries that are entirely free of ephemeris bias. This increases the inventory of trans-Neptunian objects (TNOs) with accurately known orbits by nearly 50%. Each minor planet has 20-60 Gaia/Pan-STARRS-calibrated astrometric measurements made over 2-5 oppositions, which allows accurate classification of their orbits within the trans-Neptunian dynamical populations. The populations orbiting in mean-motion resonance with Neptune are key to understanding Neptune's early migration. Our 313 resonant TNOs, including 132 plutinos, triple the available characterized sample and include new occupancy of distant resonances out to semi-major axis $a \sim 130$ au. OSSOS doubles the known population of the non-resonant Kuiper belt, providing 436 TNOs in this region, all with exceptionally high-quality orbits of $a$ uncertainty $\sigma_{a} \leq 0.1\%$; they show the belt exists from $a \gtrsim 37$ au, with a lower perihelion bound of $35$ au. We confirm the presence of a concentrated low-inclination $a\simeq 44$ au "kernel" population and a dynamically cold population extending beyond the 2:1 resonance. We finely quantify the survey's observational biases. Our survey simulator provides a straightforward way to impose these biases on models of the trans-Neptunian orbit distributions, allowing statistical comparison to the discoveries. The OSSOS TNOs, unprecedented in their orbital precision for the size of the sample, are ideal for testing concepts of the history of giant planet migration in the Solar System.

We discuss the detection in the Outer Solar System Origins Survey (OSSOS) of two objects in Neptune's distant 9:1 mean motion resonance at semimajor axis $a\approx~130$~au. Both objects are securely resonant on 10~Myr timescales, with one securely in the 9:1 resonance's leading asymmetric libration island and the other in either the symmetric or trailing asymmetric island. These objects are the largest semimajor axis objects with secure resonant classifications, and their detection in a carefully characterized survey allows for the first robust resonance population estimate beyond 100~au. The detection of these objects implies a 9:1 resonance population of $1.1\times10^4$ objects with $H_r<8.66$ ($D~\gtrsim~100$~km) on similar orbits (95\% confidence range of $\sim0.4-3\times10^4$). Integrations over 4~Gyr of an ensemble of clones spanning these objects' orbit fit uncertainties reveal that they both have median resonance occupation timescales of $\sim1$~Gyr. These timescales are consistent with the hypothesis that these objects originate in the scattering population but became transiently stuck to Neptune's 9:1 resonance within the last $\sim1$~Gyr of solar system evolution. Based on simulations of a model of the current scattering population, we estimate the expected resonance sticking population in the 9:1 resonance to be 1000-4500 objects with $H_r<8.66$; this is marginally consistent with the OSSOS 9:1 population estimate. We conclude that resonance sticking is a plausible explanation for the observed 9:1 population, but we also discuss the possibility of a primordial 9:1 population, which would have interesting implications for the Kuiper belt's dynamical history.

Hot Jupiters were the first exoplanets to be discovered around main sequence stars and astonished us with their close-in orbits. They are a prime example of how exoplanets have challenged our textbook, solar-system inspired story of how planetary systems form and evolve. More than twenty years after the discovery of the first hot Jupiter, there is no consensus on their predominant origin channel. Three classes of hot Jupiter creation hypotheses have been proposed: in situ formation, disk migration, and high eccentricity tidal migration. Although no origin channel alone satisfactorily explains all the evidence, two major origins channels together plausibly account for properties of hot Jupiters themselves and their connections to other exoplanet populations.

We present the discovery of a brown dwarf companion to the debris disk host star HR 2562. This object, discovered with the Gemini Planet Imager (GPI), has a projected separation of 20.3$\pm$0.3 au (0.618$\pm$0.004") from the star. With the high astrometric precision afforded by GPI, we have confirmed common proper motion of HR 2562B with the star with only a month time baseline between observations to more than $5\sigma$. Spectral data in $J$, $H$, and $K$ bands show morphological similarity to L/T transition objects. We assign a spectral type of L7$\pm$3 to HR 2562B, and derive a luminosity of $\log$(L$_{\rm bol}$/L$_{\odot}$)=-4.62$\pm$0.12, corresponding to a mass of 30$\pm$15 M$_{\rm Jup}$ from evolutionary models at an estimated age of the system of 300-900 Myr. Although the uncertainty in the age of the host star is significant, the spectra and photometry exhibit several indications of youth for HR 2562B. The source has a position angle consistent with an orbit in the same plane as the debris disk recently resolved with Herschel. Additionally, it appears to be interior to the debris disk. Though the extent of the inner hole is currently too uncertain to place limits on the mass of HR 2562B, future observations of the disk with higher spatial resolution may be able to provide mass constraints. This is the first brown dwarf-mass object found to reside in the inner hole of a debris disk, offering the opportunity to search for evidence of formation above the deuterium burning limit in a circumstellar disk.

Transitional disks, protoplanetary disks with deep and wide central gaps, may be the result of planetary sculpting. By comparing numerical planet-opening-gap models with observed gaps, we find systems of 3-6 giant planets are needed in order to open gaps with the observed depths and widths. We explore the dynamical stability of such multi-planet systems using N-body simulations that incorporate prescriptions for gas effects. We find they can be stable over a typical disk lifetime, with the help of eccentricity damping from the residual gap gas that facilitates planets locking into mean motion resonances. However, in order to account for the occurrence rate of transitional disks, the planet sculpting scenario demands gap-opening-friendly disk conditions, in particular, a disk viscosity $\alpha\lesssim0.001$. In addition, the demography of giant planets at $\sim 3-30$ AU separations, poorly constrained by current data, has to largely follow occurrence rates extrapolated outward from radial velocity surveys, not the lower occurrence rates extrapolated inward from direct imaging surveys. Even with the most optimistic occurrence rates, transitional disks cannot be a common phase that most gas disks experience at the end of their life, as popularly assumed, simply because there are not enough planets to open these gaps. Finally, as consequences of demanding almost all giant planets at large separations participate in transitional disk sculpting, the majority of such planets must form early and end up in a chain of mean motion resonances at the end of disk lifetime.

The giant impact phase of terrestrial planet formation establishes connections between super-Earths' orbital properties (semimajor axis spacings, eccentricities, mutual inclinations) and interior compositions (the presence or absence of gaseous envelopes). Using N-body simulations and analytic arguments, we show that spacings derive not only from eccentricities, but also from inclinations. Flatter systems attain tighter spacings, a consequence of an eccentricity equilibrium between gravitational scatterings, which increase eccentricities, and mergers, which damp them. Dynamical friction by residual disk gas plays a critical role in regulating mergers and in damping inclinations and eccentricities. Systems with moderate gas damping and high solid surface density spawn gas-enveloped super-Earths with tight spacings, small eccentricities, and small inclinations. Systems in which super-Earths coagulate without as much ambient gas, in disks with low solid surface density, produce rocky planets with wider spacings, larger eccentricities, and larger mutual inclinations. A combination of both populations can reproduce the observed distributions of spacings, period ratios, transiting planet multiplicities, and transit duration ratios exhibited by Kepler super-Earths. The two populations, both formed in situ, also help to explain observed trends of eccentricity vs. planet size, and bulk density vs. method of mass measurement (radial velocities vs. transit timing variations). Simplifications made in this study --- including the limited timespan of the simulations, and the approximate treatments of gas dynamical friction and gas depletion history --- should be improved upon in future work to enable a detailed quantitative comparison to the observations.

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

Directly detecting thermal emission from young extrasolar planets allows measurement of their atmospheric composition and luminosity, which is influenced by their formation mechanism. Using the Gemini Planet Imager, we discovered a planet orbiting the \$sim$20 Myr-old star 51 Eridani at a projected separation of 13 astronomical units. Near-infrared observations show a spectrum with strong methane and water vapor absorption. Modeling of the spectra and photometry yields a luminosity of L/LS=1.6-4.0 x 10-6 and an effective temperature of 600-750 K. For this age and luminosity, "hot-start" formation models indicate a mass twice that of Jupiter. This planet also has a sufficiently low luminosity to be consistent with the "cold- start" core accretion process that may have formed Jupiter.

We present $H$-band observations of $\beta$ Pic with the Gemini Planet Imager's (GPI's) polarimetry mode that reveal the debris disk between ~0.3" (~6 AU) and ~1.7" (~33 AU), while simultaneously detecting $\beta$ Pic $b$. The polarized disk image was fit with a dust density model combined with a Henyey-Greenstein scattering phase function. The best fit model indicates a disk inclined to the line of sight ($\phi=85.27{\deg}^{+0.26}_{-0.19}$) with a position angle $\theta_{PA}=30.35{\deg}^{+0.29}_{-0.28}$ (slightly offset from the main outer disk, $\theta_{PA}\approx29{\deg}$), that extends from an inner disk radius of $23.6^{+0.9}_{-0.6}$ AU to well outside GPI's field of view. In addition, we present an updated orbit for $\beta$ Pic $b$ based on new astrometric measurements taken in GPI's spectroscopic mode spanning 14 months. The planet has a semi-major axis of $a=9.2^{+1.5}_{-0.4}$AU, with an eccentricity $e\leq 0.26$. The position angle of the ascending node is $\Omega=31.75{\deg}\pm0.15$, offset from both the outer main disk and the inner disk seen in the GPI image. The orbital fit constrains the stellar mass of $\beta$ Pic to $1.60\pm0.05 M_{\odot}$. Dynamical sculpting by $\beta$ Pic $b$ cannot easily account for the following three aspects of the inferred disk properties: 1) the modeled inner radius of the disk is farther out than expected if caused by $\beta$ Pic b; 2) the mutual inclination of the inner disk and $\beta$ Pic $b$ is $4{\deg}$, when it is expected to be closer to zero; and 3) the aspect ratio of the disk ($h_0 = 0.137^{+0.005}_{-0.006}$) is larger than expected from interactions with $\beta$ Pic $b$ or self-stirring by the disk's parent bodies.

Planets with sizes between those of Earth and Neptune divide into two populations: purely rocky bodies whose atmospheres contribute negligibly to their sizes, and larger gas-enveloped planets possessing voluminous and optically thick atmospheres. We show that whether a planet forms rocky or gas-enveloped depends on the solid surface density of its parent disk. Assembly times for rocky cores are sensitive to disk solid surface density. Lower surface densities spawn smaller planetary embryos; to assemble a core of given mass, smaller embryos require more mergers between bodies farther apart and therefore exponentially longer formation times. Gas accretion simulations yield a rule of thumb that a rocky core must be at least 2$M_\oplus$ before it can acquire a volumetrically significant atmosphere from its parent nebula. In disks of low solid surface density, cores of such mass appear only after the gas disk has dissipated, and so remain purely rocky. Higher surface density disks breed massive cores more quickly, within the gas disk lifetime, and so produce gas-enveloped planets. We test model predictions against observations, using planet radius as an observational proxy for gas-to-rock content and host star metallicity as a proxy for disk solid surface density. Theory can explain the observation that metal-rich stars host predominantly gas-enveloped planets.

It is well-known that the light curve of a transiting planet contains information about the planet's orbital period and size relative to the host star. More recently, it has been demonstrated that a tight constraint on an individual planet's eccentricity can sometimes be derived from the light curve via the "photoeccentric effect," the effect of a planet's eccentricity on the shape and duration of its light curve. This has only been studied for large planets and high signal-to-noise scenarios, raising the question of how well it can be measured for smaller planets or low signal-to-noise cases. We explore the limits of the photoeccentric effect over a wide range of planet parameters. The method hinges upon measuring $g$ directly from the light curve, where $g$ is the ratio of the planet's speed (projected on the plane of the sky) during transit to the speed expected for a circular orbit. We find that when the signal-to-noise in the measurement of $g$ is $< 10$, the ability to measure eccentricity with the photoeccentric effect decreases. We develop a "rule of thumb" that for per-point relative photometric uncertainties $\sigma = \{ 10^{-3}, 10^{-4}, 10^{-5} \}$, the critical values of planet-star radius ratio are $R_p / R_\star \approx \{ 0.1, 0.05, 0.03 \}$ for Kepler-like 30-minute integration times. We demonstrate how to predict the best-case uncertainty in eccentricity that can be found with the photoeccentric effect for any light curve. This clears the path to study eccentricities of individual planets of various sizes in the Kepler sample and future transit surveys.

Gas giant planets orbiting within 0.1 AU of their host stars, unlikely to have formed in situ, are evidence for planetary migration. It is debated whether the typical hot Jupiter smoothly migrated inward from its formation location through the proto-planetary disk or was perturbed by another body onto a highly eccentric orbit, which tidal dissipation subsequently shrank and circularized during close stellar passages. Socrates and collaborators predicted that the latter class of model should produce a population of super-eccentric proto-hot Jupiters readily observable by Kepler. We find a paucity of such planets in the Kepler sample, inconsistent with the theoretical prediction with 96.9% confidence. Observational effects are unlikely to explain this discrepancy. We find that the fraction of hot Jupiters with orbital period P > 3 days produced by the star-planet Kozai mechanism does not exceed (at two-sigma) 44%. Our results may indicate that disk migration is the dominant channel for producing hot Jupiters with P > 3 days. Alternatively, the typical hot Jupiter may have been perturbed to a high eccentricity by interactions with a planetary rather than stellar companion and began tidal circularization much interior to 1 AU after multiple scatterings. A final alternative is that tidal circularization occurs much more rapidly early in the tidal circularization process at high eccentricities than later in the process at low eccentricities, contrary to current tidal theories.

The orbits of giant extrasolar planets often have surprisingly small semi-major axes, large eccentricities, or severe misalignments between their normals and their host stars' spin axes. In some formation scenarios invoking Kozai-Lidov oscillations, an external planetary companion drives a planet onto an orbit having these properties. The mutual inclinations for Kozai-Lidov oscillations can be large and have not been confirmed observationally. Here we deduce that observed eccentric warm Jupiters with eccentric giant companions have mutual inclinations that oscillate between 35-65 deg. Our inference is based on the pairs' observed apsidal separations, which cluster near 90 deg. The near-orthogonality of periapse directions is effected by the outer companion's quadrupolar and octupolar potentials. These systems may be undergoing a stalled version of tidal migration that produces warm Jupiters over hot Jupiters, and provide evidence for a population of multi-planet systems that are not flat and have been sculpted by Kozai-Lidov oscillations.

Numerous telescopes and techniques have been used to find and study extrasolar planets, but none has been more successful than NASA's Kepler Space Telescope. Kepler has discovered the majority of known exoplanets, the smallest planets to orbit normal stars, and the worlds most likely to be similar to our home planet. Most importantly, Kepler has provided our first look at typical characteristics of planets and planetary systems for planets with sizes as small as and orbits as large as those of the Earth.

It is debated whether the two hot Jupiter populations --- those on orbits misaligned from their host star's spin axis and those well-aligned --- result from two migration channels or from two tidal realignment regimes. Here I demonstrate that equilibrium tides raised by a planet on its star can account for three observed spin-orbit alignment trends: the aligned orbits of hot Jupiters orbiting cool stars, the planetary mass cut-off for retrograde planets, and the stratification by planet mass of cool host stars' rotation frequencies. The first trend can be caused by strong versus weak magnetic braking (the Kraft break), rather than realignment of the star's convective envelope versus the entire star. The second trend can result from a small effective stellar moment of inertia participating in the tidal realignment in hot stars, enabling massive retrograde planets to partially realign to become prograde. The third trend is attributable to higher mass planets more effectively counteracting braking to spin up their stars. Both hot and cool star require a small effective stellar moment of inertia participating in the tidal realignment, e.g., an outer layer weakly coupled to the interior. I demonstrate via Monte Carlo that this model can match the observed trends and distributions of sky-projected misalignments and stellar rotation frequencies. I discuss implications for inferring hot Jupiter migration mechanisms from obliquities, emphasizing that even hot stars do not constitute a pristine sample.

We establish the three-dimensional architecture of the Kepler-419 (previously KOI-1474) system to be eccentric yet with a low mutual inclination. Kepler-419b is a warm Jupiter at semi-major axis a = 0.370 +0.007/-0.006 AU with a large eccentricity e=0.85 +0.08/-0.07 measured via the "photoeccentric effect." It exhibits transit timing variations induced by the non-transiting Kepler-419c, which we uniquely constrain to be a moderately eccentric (e=0.184 +/- 0.002), hierarchically-separated (a=1.68 +/- 0.03 AU) giant planet (7.3 +/- 0.4 MJup). We combine sixteen quarters of Kepler photometry, radial-velocity (RV) measurements from the HIgh Resolution Echelle Spectrometer (HIRES) on Keck, and improved stellar parameters that we derive from spectroscopy and asteroseismology. From the RVs, we measure the mass of inner planet to be 2.5+/-0.3MJup and confirm its photometrically-measured eccentricity, refining the value to e=0.83+/-0.01. The RV acceleration is consistent with the properties of the outer planet derived from TTVs. We find that, despite their sizable eccentricities, the planets are coplanar to within 9 +8/-6 degrees, and therefore the inner planet's large eccentricity and close-in orbit are unlikely to be the result of Kozai migration. Moreover, even over many secular cycles, the inner planet's periapse is most likely never small enough for tidal circularization. Finally, we present and measure a transit time and impact parameter from four simultaneous ground-based light curves from 1m-class telescopes, demonstrating the feasibility of ground-based follow-up of Kepler giant planets exhibiting large TTVs.

Kepler's immense photometric precision to date was maintained through satellite stability and precise pointing. In this white paper, we argue that image modeling--fitting the Kepler-downlinked raw pixel data--can vastly improve the precision of Kepler in pointing-degraded two-wheel mode. We argue that a non-trivial modeling effort may permit continuance of photometry at 10-ppm-level precision. We demonstrate some baby steps towards precise models in both data-driven (flexible) and physics-driven (interpretably parameterized) modes. We demonstrate that the expected drift or jitter in positions in the two-weel era will help with constraining calibration parameters. In particular, we show that we can infer the device flat-field at higher than pixel resolution; that is, we can infer pixel-to-pixel variations in intra-pixel sensitivity. These results are relevant to almost any scientific goal for the repurposed mission; image modeling ought to be a part of any two-wheel repurpose for the satellite. We make other recommendations for Kepler operations, but fundamentally advocate that the project stick with its core mission of finding and characterizing Earth analogs. [abridged]

In today's mailing, Hogg et al. propose image modeling techniques to maintain 10-ppm-level precision photometry in Kepler data with only two working reaction wheels. While these results are relevant to many scientific goals for the repurposed mission, all modeling efforts so far have used a toy model of the Kepler telescope. Because the two-wheel performance of Kepler remains to be determined, we advocate for the consideration of an alternate strategy for a >1 year program that maximizes the science return from the "low-torque" fields across the ecliptic plane. Assuming we can reach the precision of the original Kepler mission, we expect to detect 800 new planet candidates in the first year of such a mission. Our proposed strategy has benefits for transit timing variation and transit duration variation studies, especially when considered in concert with the future TESS mission. We also expect to help address the first key science goal of Kepler: the frequency of planets in the habitable zone as a function of spectral type.

We present the discovery and characterization of a giant planet orbiting the young Sun-like star Kepler-63 (KOI-63, $m_{\rm Kp} = 11.6$, $T_{\rm eff} = 5576$ K, $M_\star = 0.98\, M_\odot$). The planet transits every 9.43 days, with apparent depth variations and brightening anomalies caused by large starspots. The planet's radius is $6.1 \pm 0.2 R_{\earth}$, based on the transit light curve and the estimated stellar parameters. The planet's mass could not be measured with the existing radial-velocity data, due to the high level of stellar activity, but if we assume a circular orbit we can place a rough upper bound of $120 M_{\earth}$ (3$\sigma$). The host star has a high obliquity ($\psi$ = $104^{\circ}$), based on the Rossiter-McLaughlin effect and an analysis of starspot-crossing events. This result is valuable because almost all previous obliquity measurements are for stars with more massive planets and shorter-period orbits. In addition, the polar orbit of the planet combined with an analysis of spot-crossing events reveals a large and persistent polar starspot. Such spots have previously been inferred using Doppler tomography, and predicted in simulations of magnetic activity of young Sun-like stars.

We analyze new/archival VLT/NaCo and Gemini/NICI high-contrast imaging of the young, self-luminous planet $\beta$ Pictoris b in seven near-to-mid IR photometric filters, using advanced image processing methods to achieve high signal-to-noise, high precision measurements. While $\beta$ Pic b's near-IR colors mimick that of a standard, cloudy early-to-mid L dwarf, it is overluminous in the mid-infrared compared to the field L/T dwarf sequence. Few substellar/planet-mass objects -- i.e. $\kappa$ And b and 1RXJ 1609B -- match $\beta$ Pic b's $JHK_{s}L^\prime$ photometry, and its 3.1 $\mu m$ and 5 $\mu m$ photometry are particularly difficult to reproduce. Atmosphere models adopting cloud prescriptions and large ($\sim$ 60 $\mu m$) dust grains fail to reproduce the $\beta$ Pic b spectrum. However, models incorporating thick clouds similar to those found for HR 8799 bcde but also with small (a few microns) modal particle sizes yield fits consistent with the data within uncertainties. Assuming solar abundance models, thick clouds, and small dust particles ($<a>$ = 4 $\mu m$) we derive atmosphere parameters of log(g) = 3.8 $\pm$ 0.2 and $T_{eff}$ = 1575--1650 $K$, an inferred mass of 7$^{+4}_{-3}$ $M_{J}$, and a luminosity of log(L/L$_{\odot}$) $\sim$ -3.80 $\pm$ 0.02. The best-estimated planet radius, $\approx$ 1.65 $\pm$ 0.06 $R_{J}$, is near the upper end of allowable planet radii for hot-start models given the host star's age and likely reflects challenges with constructing accurate atmospheric models. Alternatively, these radii are comfortably consistent with hot-start model predictions if $\beta$ Pic b is younger than $\approx$ 7 Myr, consistent with a late formation, well after its host star's birth $\sim$ 12$^{+8}_{-4}$ Myr ago.

Gas giants orbiting interior to the ice line are thought to have been displaced from their formation locations by processes that remain debated. Here we uncover several new metallicity trends, which together may indicate that two competing mechanisms deliver close-in giant planets: gentle disk migration, operating in environments with a range of metallicities, and violent planet-planet gravitational interactions, primarily triggered in metal-rich systems in which multiple giant planets can form. First, we show with 99.1% confidence that giant planets with semi-major axes between 0.1 and 1 AU orbiting metal-poor stars ([Fe/H]<0) are confined to lower eccentricities than those orbiting metal-rich stars. Second, we show with 93.3% confidence that eccentric proto-hot Jupiters undergoing tidal circularization primarily orbit metal-rich stars. Finally, we show that only metal-rich stars host a pile-up of hot Jupiters, helping account for the lack of such a pile-up in the overall Kepler sample. Migration caused by stellar perturbers (e.g. stellar Kozai) is unlikely to account for the trends. These trends further motivate follow-up theoretical work addressing which hot Jupiter migration theories can also produce the observed population of eccentric giant planets between 0.1 and 1 AU.

The exoplanets known as hot Jupiters---Jupiter-sized planets with periods less than 10 days---likely are relics of dynamical processes that shape all planetary system architectures. Socrates et al. (2012) argued that high eccentricity migration (HEM) mechanisms proposed for situating these close-in planets should produce an observable population of highly eccentric proto-hot Jupiters that have not yet tidally circularized. HEM should also create failed-hot Jupiters, with periapses just beyond the influence of fast circularization. Using the technique we previously presented for measuring eccentricities from photometry (the "photoeccentric effect"), we are distilling a collection of eccentric proto- and failed-hot Jupiters from the Kepler Objects of Interest (KOI). Here we present the first, KOI-1474.01, which has a long orbital period (69.7340 days) and a large eccentricity e = 0.81+0.10/-0.07, skirting the proto-hot Jupiter boundary. Combining Kepler photometry, ground-based spectroscopy, and stellar evolution models, we characterize host KOI-1474 as a rapidly-rotating F-star. Statistical arguments reveal that the transiting candidate has a low false-positive probability of 3.1%. KOI-1474.01 also exhibits transit timing variations of order an hour. We explore characteristics of the third-body perturber, which is possibly the "smoking-gun" cause of KOI-1474.01's large eccentricity. Using the host-star's rotation period, radius, and projected rotational velocity, we find KOI-1474.01's orbit is marginally consistent with aligned with the stellar spin axis, although a reanalysis is warranted with future additional data. Finally, we discuss how the number and existence of proto-hot Jupiters will not only demonstrate that hot Jupiters migrate via HEM, but also shed light on the typical timescale for the mechanism.

Exoplanet orbital eccentricities offer valuable clues about the history of planetary systems. Eccentric, Jupiter-sized planets are particularly interesting: they may link the "cold" Jupiters beyond the ice line to close-in hot Jupiters, which are unlikely to have formed in situ. To date, eccentricities of individual transiting planets primarily come from radial velocity measurements. Kepler has discovered hundreds of transiting Jupiters spanning a range of periods, but the faintness of the host stars precludes radial velocity follow-up of most. Here we demonstrate a Bayesian method of measuring an individual planet's eccentricity solely from its transit light curve using prior knowledge of its host star's density. We show that eccentric Jupiters are readily identified by their short ingress/egress/total transit durations -- part of the "photoeccentric" light curve signature of a planet's eccentricity --- even with long-cadence Kepler photometry and loosely-constrained stellar parameters. A Markov Chain Monte Carlo exploration of parameter posteriors naturally marginalizes over the periapse angle and automatically accounts for the transit probability. To demonstrate, we use three published transit light curves of HD 17156 b to measure an eccentricity of e = 0.71 +0.16/-0.09, in good agreement with the discovery value e = 0.67+/-0.08 based on 33 radial-velocity measurements. We present two additional tests using actual Kepler data. In each case the technique proves to be a viable method of measuring exoplanet eccentricities and their confidence intervals. Finally, we argue that this method is the most efficient, effective means of identifying the extremely eccentric, proto hot Jupiters predicted by Socrates et al. (2012).

Neptune's dynamical history shaped the current orbits of Kuiper belt objects (KBOs), leaving clues to the planet's orbital evolution. In the "classical" region, a population of dynamically "hot" high-inclination KBOs overlies a flat "cold" population with distinct physical properties. Simulations of qualitatively different histories for Neptune -including smooth migration on a circular orbit or scattering by other planets to a high eccentricity - have not simultaneously produced both populations. We explore a general Kuiper belt assembly model that forms hot classical KBOs interior to Neptune and delivers them to the classical region, where the cold population forms in situ. First, we present evidence that the cold population is confined to eccentricities well below the limit dictated by long-term survival. Therefore Neptune must deliver hot KBOs into the long-term survival region without excessively exciting the eccentricities of the cold population. Imposing this constraint, we explore the parameter space of Neptune's eccentricity and eccentricity damping, migration, and apsidal precession. We rule out much of parameter space, except where Neptune is scattered to a moderately eccentric orbit (e > 0.15) and subsequently migrates a distance Delta aN=1-6 AU. Neptune's moderate eccentricity must either damp quickly or be accompanied by fast apsidal precession. We find that Neptune's high eccentricity alone does not generate a chaotic sea in the classical region. Chaos can result from Neptune's interactions with Uranus, exciting the cold KBOs and placing additional constraints. Finally, we discuss how to interpret our constraints in the context of the full, complex dynamical history of the solar system.

The current dynamical structure of the Kuiper belt was shaped by the orbital evolution of the giant planets, especially Neptune, during the era following planet formation, when the giant planets may have undergone planet-planet scattering and/or planetesimal-driven migration. Numerical simulations of this process, while reproducing many properties of the belt, fail to generate the high inclinations and eccentricities observed for some objects while maintaining the observed dynamically "cold" population. We present the first of a three-part parameter study of how different dynamical histories of Neptune sculpt the planetesimal disk. Here we identify which dynamical histories allow an in situ planetesimal disk to remain dynamically cold, becoming today's cold Kuiper belt population. We find that if Neptune undergoes a period of elevated eccentricity and/or inclination, it secularly excites the eccentricities and inclinations of the planetesimal disk. We demonstrate that there are several well-defined regimes for this secular excitation, depending on the relative timescales of Neptune's migration, the damping of Neptune's orbital inclination and/or eccentricity, and the secular evolution of the planetesimals. We model this secular excitation analytically in each regime, allowing for a thorough exploration of parameter space. Neptune's eccentricity and inclination can remain high for a limited amount of time without disrupting the cold classical belt. In the regime of slow damping and slow migration, if Neptune is located (for example) at 20 AU, then its eccentricity must stay below 0.18 and its inclination below 6{\deg}.

The vertical warp in the debris disk Beta Pictoris -- an inclined inner disk extending into a flat outer disk -- has long been interpreted as the signpost of a planet on an inclined orbit. Direct images spanning 2004-2010 have revealed Beta Pictoris b, a planet with a mass and orbital distance consistent with this picture. However, it was recently reported that the orbit of planet b is aligned with the flat outer disk, not the inclined inner disk, and thus lacks the inclination to warp the disk. We explore three scenarios for reconciling the apparent misalignment of the directly imaged planet Beta Pictoris b with the warped inner disk of Beta Pictoris: observational uncertainty, an additional planet, and damping of planet b's inclination. We find that, at the extremes of the uncertainties, the orbit of Beta Pictoris b has the inclination necessary to produce the observed warp. We also find that if planet b were aligned with the flat outer disk, it would prevent another planet from creating a warp with the observed properties; therefore planet b itself must be responsible for the warp. Finally, planet b's inclination could have been damped by dynamical friction and still produce the observed disk morphology, but the feasibility of damping depends on disk properties and the presence of other planets. More precise observations of the orbit of planet b and the position angle of the outer disk will allow us to distinguish between the first and third scenario.

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

The radial velocity (RV) method for detecting extrasolar planets has been the most successful to date. The RV signal imprinted by a few Earth-mass planet around a cool star is at the limit of the typical single measurement uncertainty obtained using state-of-the-art spectrographs. This requires relying on statistics in order to unearth signals buried below noise. Artifacts introduced by observing cadences can produce spurious signals or mask genuine planets that should be easily detected otherwise. Here we discuss a particularly confusing statistical degeneracy resulting from the yearly aliasing of the first eccentric harmonic of an already-detected planet. This problem came sharply into focus after the recent announcement of the detection of a 3.1 Earth mass planet candidate in the habitable zone of the nearby low mass star GJ 581. The orbital period of the new candidate planet (GJ 581g) corresponds to an alias of the first eccentric harmonic of a previously reported planet, GJ 581d. Although the star is stable, the combination of the observing cadence and the presence of multiple planets can cause period misinterpretations. In this work, we determine whether the detection of GJ 581g is justified given this degeneracy. We also discuss the implications of our analysis for the recent Bayesian studies of the same data set, which failed to confirm the existence of the new planet. Performing a number of statistical tests, we show that, despite some caveats, the existence of GJ 581g remains the most likely orbital solution to the currently available RV data.

Radial velocity measurements of stellar reflex motion have revealed many extrasolar planets, but gaps in the observations produce aliases, spurious frequencies that are frequently confused with the planets' orbital frequencies. In the case of Gl 581 d, the distinction between an alias and the true frequency was the distinction between a frozen, dead planet and a planet possibly hospitable to life (Udry et al. 2007; Mayor et al. 2009). To improve the characterization of planetary systems, we describe how aliases originate and present a new approach for distinguishing between orbital frequencies and their aliases. Our approach harnesses features in the spectral window function to compare the amplitude and phase of predicted aliases with peaks present in the data. We apply it to confirm prior alias distinctions for the planets GJ 876 d and HD 75898 b. We find that the true periods of Gl 581 d and HD 73526 b/c remain ambiguous. We revise the periods of HD 156668 b and 55 Cnc e, which were afflicted by daily aliases. For HD 156668 b, the correct period is 1.2699 days and minimum mass is (3.1 +/- 0.4) Earth masses. For 55 Cnc e, the correct period is 0.7365 days -- the shortest of any known planet -- and minimum mass is (8.3 +/- 0.3) Earth masses. This revision produces a significantly improved 5-planet Keplerian fit for 55 Cnc, and a self-consistent dynamical fit describes the data just as well. As radial velocity techniques push to ever-smaller planets, often found in systems of multiple planets, distinguishing true periods from aliases will become increasingly important.