A planet's orbital orientation relative to an observer's line of sight determines the chord length for a transiting planet, i.e., the projected distance a transiting planet travels across the stellar disc. For a given circular orbit, the chord length determines the transit duration. Changes in the orbital inclination, the direction of the ascending node, or both, can alter this chord length and thus result in transit duration variations (TDVs). Variation of the full orbital inclination vector can even lead to de-transiting or newly transiting planets for a system. We use Laplace-Lagrange secular theory to estimate the fastest nodal eigenfrequencies for over 100 short-period planetary systems. The highest eigenfrequency is an indicator of which systems should show the strongest TDVs. We further explore five cases (TRAPPIST-1, Kepler-11, K2-138, Kepler-445, and Kepler-334) using direct N-body simulations to characterize possible TDVs and to explore whether de-transiting planets could be possible for these systems. A range of initial conditions are explored, with each realization being consistent with the observed transits. We find that tens of percent of multiplanet systems have fast enough eigenfrequencies to expect large TDVs on decade timescales. Among the directly integrated cases, we find that de-transiting planets could occur on decade timescales and TDVs of 10 minutes per decade should be common.

Many of the multi-planet systems discovered to date have been notable for their compactness, with neighbouring planets closer together than any in the Solar System. Interestingly, planet-hosting stars have a wide range of ages, suggesting that such compact systems can survive for extended periods of time. We have used numerical simulations to investigate how quickly systems go unstable in relation to the spacing between planets, focusing on hypothetical systems of Earth-mass planets on evenly-spaced orbits (in mutual Hill radii). In general, the further apart the planets are initially, the longer it takes for a pair of planets to undergo a close encounter. We recover the results of previous studies, showing a linear trend in the initial planet spacing between 3 and 8 mutual Hill radii and the logarithm of the stability time. Investigating thousands of simulations with spacings up to 13 mutual Hill radii reveals distinct modulations superimposed on this relationship in the vicinity of first and second-order mean motion resonances of adjacent and next-adjacent planets. We discuss the impact of this structure and the implications on the stability of compact multi-planet systems. Applying the outcomes of our simulations, we show that isolated systems of up to five Earth-mass planets can fit in the habitable zone of a Sun-like star without close encounters for at least $10^9$ orbits.

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.

The requirement that planetary systems be dynamically stable is often used to vet new discoveries or set limits on unconstrained masses or orbital elements. This is typically carried out via computationally expensive N-body simulations. We show that characterizing the complicated and multi-dimensional stability boundary of tightly packed systems is amenable to machine learning methods. We find that training an XGBoost machine learning algorithm on physically motivated features yields an accurate classifier of stability in packed systems. On the stability timescale investigated ($10^7$ orbits), it is 3 orders of magnitude faster than direct N-body simulations. Optimized machine learning classifiers for dynamical stability may thus prove useful across the discipline, e.g., to characterize the exoplanet sample discovered by the upcoming Transiting Exoplanet Survey Satellite (TESS). This proof of concept motivates investing computational resources to train algorithms capable of predicting stability over longer timescales and over broader regions of phase space.

We report the orbital distribution of the Trans-Neptunian objects (TNOs) discovered during the High Ecliptic Latitude (HiLat) extension of the Canada-France Ecliptic Plane Survey (CFEPS), conducted from June 2006 to July 2009. The HiLat component was designed to address one of the shortcomings of ecliptic surveys (like CFEPS), their lack of sensitivity to high-inclination objects. We searched 701~deg$^2$ of sky ranging from 12$^\circ$ to 85$^\circ$ ecliptic latitude and discovered \lKBO TNOs, with inclinations between 15$^\circ$ to 104$^\circ$. This survey places a very strong constraint on the inclination distribution of the hot component of the classical Kuiper Belt, ruling out any possibility of a large intrinsic fraction of highly inclined orbits. Using the parameterization of \citet{2001AJ....121.2804B}, the HiLat sample combined with CFEPS imposes a width $14^\circ \le \sigma \le 15.5^\circ$, with a best match for $\sigma = 14.5^\circ$. HiLat discovered the first retrograde TNO, 2008~KV$_{42}$, with an almost polar orbit with inclination 104$^\circ$, and (418993), a scattering object with perihelion in the region of Saturn's influence, with $a \sim 400$~AU and $i = 68^\circ$.

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 present the sixth catalog of Kepler candidate planets based on nearly 4 years of high precision photometry. This catalog builds on the legacy of previous catalogs released by the Kepler project and includes 1493 new Kepler Objects of Interest (KOIs) of which 554 are planet candidates, and 131 of these candidates have best fit radii <1.5 R_earth. This brings the total number of KOIs and planet candidates to 7305 and 4173 respectively. We suspect that many of these new candidates at the low signal-to-noise limit may be false alarms created by instrumental noise, and discuss our efforts to identify such objects. We re-evaluate all previously published KOIs with orbital periods of >50 days to provide a consistently vetted sample that can be used to improve planet occurrence rate calculations. We discuss the performance of our planet detection algorithms, and the consistency of our vetting products. The full catalog is publicly available at the NASA Exoplanet Archive.

The first discoveries of exoplanets around Sun-like stars have fueled efforts to find ever smaller worlds evocative of Earth and other terrestrial planets in the Solar System. While gas-giant planets appear to form preferentially around metal-rich stars, small planets (with radii less than four Earth radii) can form under a wide range of metallicities. This implies that small, including Earth-size, planets may have readily formed at earlier epochs in the Universe's history when metals were far less abundant. We report Kepler spacecraft observations of KOI-3158, a metal-poor Sun-like star from the old population of the Galactic thick disk, which hosts five planets with sizes between Mercury and Venus. We used asteroseismology to directly measure a precise age of 11.2+/-1.0 Gyr for the host star, indicating that KOI-3158 formed when the Universe was less than 20% of its current age and making it the oldest known system of terrestrial-size planets. We thus show that Earth-size planets have formed throughout most of the Universe's 13.8-billion-year history, providing scope for the existence of ancient life in the Galaxy.

The Kepler mission discovered 2842 exoplanet candidates with 2 years of data. We provide updates to the Kepler planet candidate sample based upon 3 years (Q1-Q12) of data. Through a series of tests to exclude false-positives, primarily caused by eclipsing binary stars and instrumental systematics, 855 additional planetary candidates have been discovered, bringing the total number known to 3697. We provide revised transit parameters and accompanying posterior distributions based on a Markov Chain Monte Carlo algorithm for the cumulative catalogue of Kepler Objects of Interest. There are now 130 candidates in the cumulative catalogue that receive less than twice the flux the Earth receives and more than 1100 have a radius less than 1.5 Rearth. There are now a dozen candidates meeting both criteria, roughly doubling the number of candidate Earth analogs. A majority of planetary candidates have a high probability of being bonafide planets, however, there are populations of likely false-positives. We discuss and suggest additional cuts that can be easily applied to the catalogue to produce a set of planetary candidates with good fidelity. The full catalogue is publicly available at the NASA Exoplanet Archive.

The chemical composition of stars hosting small exoplanets (with radii less than four Earth radii) appears to be more diverse than that of gas-giant hosts, which tend to be metal-rich. This implies that small, including Earth-size, planets may have readily formed at earlier epochs in the Universe's history when metals were more scarce. We report Kepler spacecraft observations of Kepler-444, a metal-poor Sun-like star from the old population of the Galactic thick disk and the host to a compact system of five transiting planets with sizes between those of Mercury and Venus. We validate this system as a true five-planet system orbiting the target star and provide a detailed characterization of its planetary and orbital parameters based on an analysis of the transit photometry. Kepler-444 is the densest star with detected solar-like oscillations. We use asteroseismology to directly measure a precise age of 11.2+/-1.0 Gyr for the host star, indicating that Kepler-444 formed when the Universe was less than 20% of its current age and making it the oldest known system of terrestrial-size planets. We thus show that Earth-size planets have formed throughout most of the Universe's 13.8-billion-year history, leaving open the possibility for the existence of ancient life in the Galaxy. The age of Kepler-444 not only suggests that thick-disk stars were among the hosts to the first Galactic planets, but may also help to pinpoint the beginning of the era of planet formation.

Earth-scale planets in the classical habitable zone (HZ) are more likely to be habitable if they possess active geophysics. Without a constant internal energy source, planets cool as they age, eventually terminating tectonic activity and rendering the planet sterile to life. However, for planets orbiting low-mass stars, the presence of an outer companion could generate enough tidal heat in the HZ planet to prevent such cooling. The range of mass and orbital parameters for the companion that give adequate long-term heating of the inner HZ planet, while avoiding very early total desiccation, is probably substantial. We locate the ideal location for the outer of a pair of planets, under the assumption that the inner planet has the same incident flux as Earth, orbiting example stars: a generic late M dwarf ($T_{eff}=2670 K$) and the M9V/L0 dwarf DEN1048. Thus discoveries of Earth-scale planets in the HZ zone of old small stars should be followed by searches for outer companion planets that might be essential for current habitability.

Classical secular theory can be a powerful tool to describe the qualitative character of multi-planet systems and offer insight into their histories. The eigenmodes of the secular behavior, rather than current orbital elements, can help identify tidal effects, early planet-planet scattering, and dynamical coupling among the planets, for systems in which mean-motion resonances do not play a role. Although tidal damping can result in aligned major axes after all but one eigenmode have damped away, such alignment may simply be fortuitous. An example of this is 55 Cancri (orbital solution of Fischer et al., 2008) where multiple eigenmodes remain undamped. Various solutions for 55 Cancri are compared, showing differing dynamical groupings, with implications for the coupling of eccentricities and for the partitioning of damping among the planets. Solutions for orbits that include expectations of past tidal evolution with observational data, must take into account which eigenmodes should be damped, rather than expecting particular eccentricities to be near zero. Classical secular theory is only accurate for low eccentricity values, but comparison with other results suggests that it can yield useful qualitative descriptions of behavior even for moderately large eccentricity values, and may have advantages for revealing underlying physical processes and, as large numbers of new systems are discovered, for triage to identify where more comprehensive dynamical studies should have priority.

Although warm jupiters are generally too far from their stars for tides to be important, the presence of an inner planetary companion to a warm jupiter can result in tidal evolution of the system. Insight into the process and its effects comes form classical secular theory of planetary perturbations. The lifetime of the inner planet may be shorter than the age of the system, because the warm jupiter maintains its eccentricity and hence promotes tidal migration into the star. Thus a warm jupiter observed to be alone in its system might have previously cleared away any interior planets. Before its demise, even if an inner planet is of terrestrial scale, it may promote damping of the warm jupiter's eccentricity. Thus any inferences of the initial orbit of an observed warm jupiter must include the possibility of a greater initial eccentricity than would be estimated by assuming it had always been alone. Tidal evolution involving multiple planets also enhances the internal heating of the planets, which readily exceeds that of stellar radiation for the inner planet, and may be great enough to affect the internal structure of warm jupiters. Secular theory gives insight into the tidal processes, providing, among other things, a way to constrain eccentricities of transiting planets based on estimates of the tidal parameter Q.

In a multiplanet system, tides acting on the inner planet can significantly affect the orbital evolution of the entire system. While tides usually damp eccentricities, a novel mechanism identified by Correia et al. (2012) tends to raise eccentricities as a result of the tides' effect on the inner planet's rotation. Our analytical description of this spin-driven tidal (SDT) effect shows that, while the inner planet's eccentricity undergoes pumping, the process is more completely described by an exchange of strength between the two eigenmodes of the dynamical system. Our analysis allows derivation of criteria for two-planet coplanar systems where the SDT effect can reverse tidal damping, and may preclude the effect's being significant for realistic systems. For the specific case quantified by Correia et al., the effect is strong because of the large adopted tidal time lag, which may not be appropriate for the assumed Saturn-like inner planet. On the other hand, the effective Q for any given planet in exotic circumstances is very uncertain, so the SDT effect could play a role in planetary evolution.

The transneptunian objects (TNOs) trapped in mean-motion resonances with Neptune were likely emplaced there during planet migration late in the giant-planet formation process. We perform detailed modelling of the resonant objects detected in the Canada-France Ecliptic Plane Survey (CFEPS) in order to provide population estimates and, for some resonances, constrain the complex internal orbital element distribution. Detection biases play a critical role because phase relationships with Neptune make object discovery more likely at certain longitudes. This paper discusses the 3:2, 5:2, 2:1, 3:1, 5:1, 4:3, 5:3, 7:3, 5:4, and 7:4 mean-motion resonances, all of which had CFEPS detections, along with our upper limit on 1:1 Neptune Trojans (which is consistent with their small population estimated elsewhere). For the plutinos (TNOs in the 3:2 resonance) we refine the orbital element distribution given in Kavelaars et al. (2009) and show that steep H-magnitude distributions (N(H) proportional to 10aH, with a=0.8-0.9) are favoured in the range Hg=8-9, and confirm that this resonance does not share the inclination distribution of the classical Kuiper Belt. We give the first population estimate for the 5:2 resonance and find that, to within the uncertainties, the population is equal to that of the 3:2 (13,000 TNOs with Hg < 9.16), whereas the 2:1 population is smaller by a factor of 3-4 compared to the other two resonances. We also measure significant populations inhabiting the 4:3, 5:3, 7:3, 5:4, 7:4, 3:1, and 5:1 resonances, with Hg < 9.16 (D >100 km) populations in the thousands. We compare our intrinsic population and orbital-element distributions with several published models of resonant-TNO production; the most striking discrepancy is that resonances beyond the 2:1 are in reality more heavily populated than in published models.

Tidal damping of one of the orbits in a planetary system can lead to aligned major-axes (the so-called "fixed-point" condition), but currently aligned major axes do not necessarily imply such a history. An example is the nominal orbital solution for the 61 Virginis system where two orbits librate about alignment, but evaluation of the eigenmodes of the secular theory shows it could not be the result of tidal damping but rather of initial conditions. Nevertheless, the amplitudes of the eigenmodes suggest that this system may have undergone some degree of tidal damping.

We report the orbital distribution of the trans-neptunian objects (TNOs) discovered during the Canada-France Ecliptic Plane Survey, whose discovery phase ran from early 2003 until early 2007. The follow-up observations started just after the first discoveries and extended until late 2009. We obtained characterized observations of 321 sq.deg. of sky to depths in the range g ~ 23.5--24.4 AB mag. We provide a database of 169 TNOs with high-precision dynamical classification and known discovery efficiency. Using this database, we find that the classical belt is a complex region with sub-structures that go beyond the usual splitting of inner (interior to 3:2 mean-motion resonance [MMR]), outer (exterior to 2:1 MMR), and main (in between). The main classical belt (a=40--47 AU) needs to be modeled with at least three components: the `hot' component with a wide inclination distribution and two `cold' components (stirred and kernel) with much narrower inclination distributions. The hot component must have a significantly shallower absolute magnitude (Hg) distribution than the other two components. With 95% confidence, there are 8000+1800-1600 objects in the main belt with Hg <= 8.0, of which 50% are from the hot component, 40% from the stirred component and 10% from the kernel; the hot component's fraction drops rapidly with increasing Hg. Because of this, the apparent population fractions depend on the depth and ecliptic latitude of a trans-neptunian survey. The stirred and kernel components are limited to only a portion of the main belt, while we find that the hot component is consistent with a smooth extension throughout the inner, main and outer regions of the classical belt; the inner and outer belts are consistent with containing only hot-component objects. The Hg <= 8.0 TNO population estimates are 400 for the inner belt and 10,000 for the outer belt within a factor of two.

In a multi-planet system, a gradual change in one planet's semi-major axis will affect the eccentricities of all the planets, as angular momentum is distributed via secular interactions. If tidal dissipation in the planet is the cause of the change in semi-major axis, it also damps that planet's eccentricity, which in turn also contributes to the evolution of all the eccentricities. Formulae quantifying the combined effects on the whole system due to semi-major axis changes, as well as eccentricity damping, are derived here for a two-planet system. The CoRoT 7 system is considered as an example.

We explore the origin and orbital evolution of the Kuiper belt in the framework of a recent model of the dynamical evolution of the giant planets, sometimes known as the Nice model. This model is characterized by a short, but violent, instability phase, during which the planets were on large eccentricity orbits. One characteristic of this model is that the proto-planetary disk must have been truncated at roughly 30 to 35 AU so that Neptune would stop migrating at its currently observed location. As a result, the Kuiper belt would have initially been empty. In this paper we present a new dynamical mechanism which can deliver objects from the region interior to ~35 AU to the Kuiper belt without excessive inclination excitation. Assuming that the last encounter with Uranus delivered Neptune onto a low-inclination orbit with a semi-major axis of ~27 AU and an eccentricity of ~0.3, and that subsequently Neptune's eccentricity damped in ~1 My, our simulations reproduce the main observed properties of the Kuiper belt at an unprecedented level.