The minor planets on orbits that are dynamically stable in Neptune's 1:1 resonance on Gyr timescales were likely emplaced by Neptune's outward migration. We explore the intrinsic libration amplitude, eccentricity, and inclination distribution of Neptune's stable Trojans, using the detections and survey efficiency of the Outer Solar System Origins Survey (OSSOS) and Pan-STARRS1. We find that the libration amplitude of the stable Neptunian Trojan population can be well modeled as a Rayleigh distribution with a libration amplitude width $\sigma_{A_phi}$ of 15$^\circ$. When taken as a whole, the Neptune Trojan population can be acceptably modeled with a Rayleigh eccentricity distribution of width $\sigma_e$ of 0.045 and a typical sin(i) x gaussian inclination distribution with a width $\sigma_i$ of 14 +/- 2 degrees. However, these distributions are only marginally acceptable. This is likely because, even after accounting for survey detection biases, the known large Hr < 8 and small Hr >= 8 Neptune Trojans appear to have markedly different eccentricities and inclinations. We propose that like the classical Kuiper belt, the stable intrinsic Neptunian Trojan population have dynamically `hot' and dynamically `cold' components to its eccentricity/inclination distribution, with $\sigma_{e-cold}$ ~ 0.02 / $\sigma_{i-cold}$ ~ 6$^\circ$ and $\sigma_{e-hot}$~ 0.05 / $\sigma_{i-hot}$ ~ 18$^\circ$. In this scenario, the `cold' L4 Neptunian Trojan population lacks the Hr >= 8 member and has 13 +11/-6 `cold' Trojans with Hr < 8. On the other hand, the `hot' L4 Neptunian Trojan population has 136 +57/-48 Trojans with Hr < 10 -- a population 2.4 times greater than that of the L4 Jovian Trojans in the same luminosity range.

Observations show that 100-km-class Kuiper belt objects (KBOs) can be divided in (at least) two color groups, hereafter red (R, g-i<1.2) and very red (VR, g-i>1.2), reflecting a difference in their surface composition. This is thought to imply that KBOs formed over a relatively wide range of radial distance, r. The cold classicals at 42<r<47 au are predominantly VR and known Neptune Trojans at r=30 au are mostly R. Intriguingly, however, the dynamically hot KBOs show a mix of R and VR colors and no correlation of color with r. Here we perform migration/instability simulations where the Kuiper belt is populated from an extended planetesimal disk. We find that the color observations can be best understood if R objects formed at r<r* and VR objects at r>r*, with 30<r*<40 au. The proposed transition at 30<r*<40 au would explain why the VR objects in the dynamically hot population have smaller orbital inclinations than the R objects, because the orbital excitation from Neptune weakens for orbits starting beyond 30 au. Possible causes of the R-VR color bimodality are discussed.

The surface characterization of Trans-Neptunian Binaries (TNBs) is key to understanding the properties of the disk of planetesimals from which these objects formed. In the optical wavelengths, it has been demonstrated that most equal-sized component systems share similar colors, suggesting they have a similar composition. The color homogeneity of binary pairs contrasts with the overall diversity of colors in the Kuiper belt, which was interpreted as evidence that Trans-Neptunian Objects (TNOs) formed from a locally homogeneous and globally heterogeneous protoplanetary disk. In this paradigm, binary pairs must have formed early, before the dynamically hot TNOs were scattered out from their formation location. The latter inferences, however, relied on the assumption that the matching colors of the binary components imply matching composition. Here, we test this assumption by examining the component-resolved photometry of three TNBs found in the Outer Solar System Origins Survey: 505447 (2013 SQ99), 511551 (2014 UD225) and 506121 (2016 BP81), across the visible and J-band near-infrared wavelength range. We report similar colors within 2 sigma for the binary pairs suggestive of similar reflectance spectra and hence surface composition. This advocates for gravitational collapse of pebble clouds as a possible TNO formation route. We however stress that several similarly small TNOs, including at least one binary, have been shown to exhibit substantial spectral variability in the near-infrared, implying color equality of binary pairs is likely to be violated in some cases.

The Cold Classical Kuiper Belt, a class of small bodies in undisturbed orbits beyond Neptune, are primitive objects preserving information about Solar System formation. The New Horizons spacecraft flew past one of these objects, the 36 km long contact binary (486958) Arrokoth (2014 MU69), in January 2019. Images from the flyby show that Arrokoth has no detectable rings, and no satellites (larger than 180 meters diameter) within a radius of 8000 km, and has a lightly-cratered smooth surface with complex geological features, unlike those on previously visited Solar System bodies. The density of impact craters indicates the surface dates from the formation of the Solar System. The two lobes of the contact binary have closely aligned poles and equators, constraining their accretion mechanism.

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.

We use published models of the early Solar System evolution with a slow, long-range and grainy migration of Neptune to predict the orbital element distributions and the number of modern-day Centaurs. The model distributions are biased by the OSSOS survey simulator and compared with the OSSOS Centaur detections. We find an excellent match to the observed orbital distribution, including the wide range of orbital inclinations which was the most troublesome characteristic to fit in previous models. A dynamical model, in which the original population of outer disk planetesimals was calibrated from Jupiter Trojans, is used to predict that OSSOS should detect 11+/-4 Centaurs with semimajor axis a<30 au, perihelion distance q>7.5 au and diameter D>10 km (absolute magnitude H_r<13.7 for a 6% albedo). This is consistent with 15 actual OSSOS Centaur detections with H_r<13.7. The population of Centaurs is estimated to be 21,000+/-8,000 for D>10 km. The inner scattered disk at 50<a<200 au should contain (2.0+/-0.8)x10^7 D>10 km bodies and the Oort cloud should contain (5.0+/-1.9)x10^8 D>10 km comets. Population estimates for different diameter cutoffs can be obtained from the size distribution of Jupiter Trojans (N(>D) proportional to D^(-2.1) for 5<D<100 km). We discuss model predictions for the Large Synoptic Survey Telescope observations of Centaurs.

Most known trans-Neptunian objects (TNOs) gravitationally scattering off the giant planets have orbital inclinations consistent with an origin from the classical Kuiper belt, but a small fraction of these "scattering TNOs" have inclinations that are far too large (i > 45 deg) for this origin. These scattering outliers have previously been proposed to be interlopers from the Oort cloud or evidence of an undiscovered planet. Here we test these hypotheses using N-body simulations and the 69 centaurs and scattering TNOs detected in the Outer Solar Systems Origins Survey and its predecessors. We confirm that observed scattering objects cannot solely originate from the classical Kuiper belt, and we show that both the Oort cloud and a distant planet generate observable highly inclined scatterers. Although the number of highly inclined scatterers from the Oort Cloud is ~3 times less than observed, Oort cloud enrichment from the Sun's galactic migration or birth cluster could resolve this. Meanwhile, a distant, low-eccentricity 5 Earth-mass planet replicates the observed fraction of highly inclined scatterers, but the overall inclination distribution is more excited than observed. Furthermore, the distant planet generates a longitudinal asymmetry among detached TNOs that is less extreme than often presumed, and its direction reverses across the perihelion range spanned by known TNOs. More complete models that explore the dynamical origins of the planet are necessary to further study these features. With observational biases well-characterized, our work shows that the orbital distribution of detected scattering bodies is a powerful constraint on the unobserved distant solar system.

The migration of Neptune's resonances through the proto-Kuiper belt has been imprinted in the distribution of small bodies in the outer Solar System. Here we analyze five published Neptune migration models in detail, focusing on the high pericenter distance (high-q) trans-Neptunian Objects (TNOs) near Neptune's 5:2 and 3:1 mean-motion resonances, because they have large resonant populations, are outside the main classical belt, and are relatively isolated from other strong resonances. We compare the observationally biased output from these dynamical models with the detected TNOs from the Outer Solar System Origins Survey, via its Survey Simulator. All of the four new OSSOS detections of high-q non-resonant TNOs are on the Sunward side of the 5:2 and 3:1 resonances. We show that even after accounting for observation biases, this asymmetric distribution cannot be drawn from a uniform distribution of TNOs at 2sigma confidence. As shown by previous work, our analysis here tentatively confirms that the dynamical model that uses grainy slow Neptune migration provides the best match to the real high-q TNO orbital data. However, due to extreme observational biases, we have very few high-q TNO discoveries with which to statistically constrain the models. Thus, this analysis provides a framework for future comparison between the output from detailed, dynamically classified Neptune migration simulations and the TNO discoveries from future well-characterized surveys. We show that a deeper survey (to a limiting r-magnitude of 26.0) with a similar survey area to OSSOS could statistically distinguish between these five Neptune migration models.

The Colours of the Outer Solar System Origins Survey (Col-OSSOS) is acquiring near-simultaneous $g$, $r$, and $J$ photometry of unprecedented precision with the Gemini North Telescope, targeting nearly a hundred trans-Neptunian objects (TNOs) brighter than $m_r=23.6$ mag discovered in the Outer Solar System Origins Survey. Combining the optical and near-infrared photometry with the well-characterized detection efficiency of the Col-OSSOS target sample will provide the first flux-limited compositional dynamical map of the outer Solar System. In this paper, we describe our observing strategy and detail the data reduction processes we employ, including techniques to mitigate the impact of rotational variability. We present optical and near-infrared colors for 35 TNOs. We find two taxonomic groups for the dynamically excited TNOs, the neutral and red classes, which divide at $g-r \simeq 0.75$. Based on simple albedo and orbital distribution assumptions, we find that the neutral class outnumbers the red class, with a ratio of 4:1 and potentially as high as 11:1. Including in our analysis constraints from the cold classical objects, which are known to exhibit unique albedos and $r-z$ colors, we find that within our measurement uncertainty, our observations are consistent with the primordial Solar System protoplanetesimal disk being neutral-class-dominated, with two major compositional divisions in $grJ$ color space.

We present variability measurements and partial light curves of Trans-Neptunian Objects (TNOs) from a two-night pilot study using Hyper Suprime-Cam (HSC) on the Subaru Telescope (Maunakea, Hawai'i, USA). Subaru's large aperture (8-m) and HSC's large field of view (1.77 square degrees) allow us to obtain measurements of multiple objects with a range of magnitudes in each telescope pointing. We observed 65 objects with m_r = 22.6--25.5 mag in just six pointings, allowing 20--24 visits of each pointing over the two nights. Our sample, all discovered in the recent Outer Solar System Origins Survey (OSSOS), span absolute magnitudes H_r = 6.2--10.8 mag and thus investigates smaller objects than previous light curve projects have typically studied. Our data supports the existence of a correlation between light curve amplitude and absolute magnitude seen in other works, but does not support a correlation between amplitude and orbital inclination. Our sample includes a number of objects from different dynamical populations within the trans-Neptunian region, but we do not find any relationship between variability and dynamical class. We were only able to estimate periods for 12 objects in the sample and found that a longer baseline of observations is required for reliable period analysis. We find that 31 objects (just under half of our sample) have variability greater than 0.4 magnitudes during all of the observations; in smaller 1.25 hr, 1.85 hr and 2.45 hr windows, the median variability is 0.13, 0.16 and 0.19 mags, respectively. The fact that variability on this scale is common for small TNOs has important implications for discovery surveys (such as OSSOS or the Large Synoptic Survey Telescope) and color measurements.

Both physical and dynamical properties must be considered to constrain the origins of the dynamically excited distant Solar System populations. We present high-precision (g-r) colors for 25 small (Hr>5) dynamically excited Trans-Neptunian Objects (TNOs) and centaurs acquired as part of the Colours of the Outer Solar System Origins Survey (Col-OSSOS). We combine our dataset with previously published measurements and consider a set of 229 colors of outer Solar System objects on dynamically excited orbits. The overall color distribution is bimodal and can be decomposed into two distinct classes, termed `gray' and `red', that each has a normal color distribution. The two color classes have different inclination distributions: red objects have lower inclinations than the gray ones. This trend holds for all dynamically excited TNO populations. Even in the worst-case scenario, biases in the discovery surveys cannot account for this trend: it is intrinsic to the TNO population. Considering that TNOs are the precursors of centaurs, and that their inclinations are roughly preserved as they become centaurs, our finding solves the conundrum of centaurs being the only outer Solar System population identified so far to exhibit this property (Tegler et al. 2016). The different orbital distributions of the gray and red dynamically excited TNOs provide strong evidence that their colors are due to different formation locations in a disk of planetesimals with a compositional gradient.

The ESA Euclid mission is a space telescope that will survey ~15,000 square degrees of the sky, primarily to study the distant universe (constraining cosmological parameters through the lensing of galaxies). It is also expected to observe ~150,000 Solar System Objects (SSOs), primarily in poorly understood high inclination populations, as it will mostly avoid +/-15 degrees from the ecliptic plane. With a launch date of 2022 and a 6 year survey, Euclid and LSST will operate at the same time, and have complementary capabilities. We propose a LSST mini-survey to coordinate quasi-simultaneous observations between these two powerful observatories, when possible, with the primary aim of greatly improving the orbits of SSOs discovered by these facilities. As Euclid will operate from a halo orbit around the Sun-Earth L2 Lagrangian point, there will be significant parallax between observations from Earth and Euclid (0.01 AU). This means that simultaneous observations will give an independent distance measurement to SSOs, giving additional constraints on orbits compared to single Euclid visits.

Context. Centaurs are icy objects in transition between the transneptunian region and the inner solar system, orbiting the Sun in the giant planet region. Some Centaurs display cometary activity, which cannot be sustained by the sublimation of water ice in this part of the solar system, and has been hypothesized to be due to the crystallization of amorphous water ice. Aims. In this work, we look at Centaurs discovered by the Outer Solar System Origins Survey (OSSOS) and search for cometary activity. Tentative detections would improve understanding of the origins of activity among these objects. Methods. We search for comae and structures by fitting and subtracting both Point Spread Functions (PSF) and Trailed point-Spread Functions (TSF) from the OSSOS images of each Centaur. When available, Col-OSSOS images were used to search for comae too. Results. No cometary activity is detected in the OSSOS sample. We track the recent orbital evolution of each new Centaur to confirm that none would actually be predicted to be active, and we provide size estimates for the objects. Conclusions. The addition of 20 OSSOS objects to the population of 250 known Centaurs is consistent with the currently understood scenario, in which drastic drops in perihelion distance induce changes in the thermal balance prone to trigger cometary activity in the giant planet region.

We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

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.

All surveys include observational biases, which makes it impossible to directly compare properties of discovered trans-Neptunian Objects (TNOs) with dynamical models. However, by carefully keeping track of survey pointings on the sky, detection limits, tracking fractions, and rate cuts, the biases from a survey can be modelled in Survey Simulator software. A Survey Simulator takes an intrinsic orbital model (from, for example, the output of a dynamical Kuiper belt emplacement simulation) and applies the survey biases, so that the biased simulated objects can be directly compared with real discoveries. This methodology has been used with great success in the Outer Solar System Origins Survey (OSSOS) and its predecessor surveys. In this chapter, we give four examples of ways to use the OSSOS Survey Simulator to gain knowledge about the true structure of the Kuiper Belt. We demonstrate how to statistically compare different dynamical model outputs with real TNO discoveries, how to quantify detection biases within a TNO population, how to measure intrinsic population sizes, and how to use upper limits from non-detections. We hope this will provide a framework for dynamical modellers to statistically test the validity of their models.

The scattering trans-Neptunian Objects (TNOs) can be measured to smaller sizes than any other distant small-body population. We use the largest sample yet obtained, 68 discoveries, primarily by the Outer Solar System Origins Survey (OSSOS), to constrain the slope of its luminosity distribution, with sensitivity to much fainter absolute $H$ magnitudes than previous work. Using the analysis technique in Shankman et al. (2016), we confirm that a single slope for the $H$-distribution is not an accurate representation of the scattering TNOs and Centaurs, and that a break in the distribution is required, in support of previous conclusions. A bright-end slope of $\alpha_b=0.9$ transitioning to a faint-end slope $\alpha_f$ of 0.4-0.5 with a differential number contrast $c$ from 1 (a knee) to 10 (a divot) provides an acceptable match to our data. We find that break magnitudes $H_b$ of 7.7 and 8.3, values both previously suggested for dynamically hot Kuiper belt populations, are equally non-rejectable for a range of $\alpha_f$ and $c$ in our statistical analysis. Our preferred divot $H$-distribution transitions to $\alpha_f=0.5$ with a divot of contrast $c=3$ at $H_b=8.3$, while our preferred knee $H$-distribution transitions to $\alpha_f=0.4$ at $H_b=7.7$. The intrinsic population of scattering TNOs required to match the OSSOS detections is $3\times10^6$ for $H_r<12$, and $9\times10^4$ for $H_r<8.66$ ($D\gtrsim100$~km), with Centaurs having an intrinsic population two orders of magnitude smaller.

Several different classes of trans-Neptunian objects (TNOs) have been identified based on their optical and near-infrared colors. As part of the Colours of the Outer Solar System Origins Survey, we have obtained $g$, $r$, and $z$ band photometry of 26 TNOs using Subaru and Gemini Observatories. Previous color surveys have not utilized $z$ band reflectance, and the inclusion of this band reveals significant surface reflectance variations between sub-populations. The colors of TNOs in $g-r$ and $r-z$ show obvious structure, and appear consistent with the previously measured bi-modality in $g-r$. The distribution of colors of the two dynamically excited surface types can be modeled using the two-component mixing models from Fraser \& Brown (2012). With the combination of $g-r$ and $r-z$, the dynamically excited classes can be separated cleanly into red and neutral surface classes. In $g - r$ and $r - z$, the two dynamically excited surface groups are also clearly distinct from the cold classical TNO surfaces, which are red, with $g-r\gtrsim$0.85 and $r-z\lesssim$0.6, while all dynamically excited objects with similar $g-r$ colors exhibit redder $r-z$ colors. The $z$ band photometry makes it possible for the first time to differentiate the red excited TNO surfaces from the red cold classical TNO surfaces. The discovery of different $r-z$ colors for these cold classical TNOs makes it possible to search for cold classical surfaces in other regions of the Kuiper belt and to completely separate cold classical TNOs from the dynamically excited population, which overlaps in orbital parameter space.

The accumulating, but small, set of large semi-major axis trans-Neptunian objects (TNOs) shows an apparent clustering in the orientations of their orbits. This clustering must either be representative of the intrinsic distribution of these TNOs, or else arise as a result of observation biases and/or statistically expected variations for such a small set of detected objects. The clustered TNOs were detected across different and independent surveys, which has led to claims that the detections are therefore free of observational bias. This apparent clustering has led to the so-called "Planet 9" hypothesis that a super-Earth currently resides in the distant solar system and causes this clustering. The Outer Solar System Origins Survey (OSSOS) is a large program that ran on the Canada-France-Hawaii Telescope from 2013--2017, discovering more than 800 new TNOs. One of the primary design goals of OSSOS was the careful determination of observational biases that would manifest within the detected sample. We demonstrate the striking and non-intuitive biases that exist for the detection of TNOs with large semi-major axes. The eight large semi-major axis OSSOS detections are an independent dataset, of comparable size to the conglomerate samples used in previous studies. We conclude that the orbital distribution of the OSSOS sample is consistent with being detected from a uniform underlying angular distribution.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

The cold classical Kuiper belt objects have low inclinations and eccentricities and are the only Kuiper belt population suspected to have formed in situ. Compared with the dynamically excited populations, which exhibit a broad range of colours and a low binary fraction of ~10% cold classical Kuiper belt objects typically have red optical colours with ~30% of the population found in binary pairs; the origin of these differences remains unclear. We report the detection of a population of blue-coloured, tenuously bound binaries residing among the cold classical Kuiper belt objects. Here we show that widely separated binaries could have survived push-out into the cold classical region during the early phases of Neptune's migration. The blue binaries may be contaminants, originating at ~38 au, and could provide a unique probe of the formative conditions in a region now nearly devoid of objects. The idea that the blue objects, which are predominantly binary, are the products of push-out requires that the planetesimals formed entirely as multiples. Plausible formation routes include planetesimal formation via pebble accretion and subsequent binary production through dynamic friction and binary formation during the collapse of a cloud of solids.

We report the discovery of the minor planet 2013 SY$_{99}$, on an exceptionally distant, highly eccentric orbit. With a perihelion of 50.0 au, 2013 SY$_{99}$'s orbit has a semi-major axis of $730 \pm 40$ au, the largest known for a high-perihelion trans-Neptunian object (TNO), well beyond those of (90377) Sedna and 2012 VP$_{113}$. Yet, with an aphelion of $1420 \pm 90$ au, 2013 SY$_{99}$'s orbit is interior to the region influenced by Galactic tides. Such TNOs are not thought to be produced in the current known planetary architecture of the Solar System, and they have informed the recent debate on the existence of a distant giant planet. Photometry from the Canada-France-Hawaii Telescope, Gemini North and Subaru indicate 2013 SY$_{99}$ is $\sim 250$ km in diameter and moderately red in colour, similar to other dynamically excited TNOs. Our dynamical simulations show that Neptune's weak influence during 2013 SY$_{99}$'s perihelia encounters drives diffusion in its semi-major axis of hundreds of astronomical units over 4 Gyr. The overall symmetry of random walks in semi-major axis allow diffusion to populate 2013 SY$_{99}$'s orbital parameter space from the 1000-2000 au inner fringe of the Oort cloud. Diffusion affects other known TNOs on orbits with perihelia of 45 to 49 au and semi-major axes beyond 250 au, providing a formation mechanism that implies an extended population, gently cycling into and returning from the inner fringe of the Oort cloud.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

We report the discovery and orbit of a new dwarf planet candidate, 2015 RR$_{245}$, by the Outer Solar System Origins Survey (OSSOS). 2015 RR$_{245}$'s orbit is eccentric ($e=0.586$), with a semi-major axis near 82 au, yielding a perihelion distance of 34 au. 2015 RR$_{245}$ has $g-r = 0.59 \pm 0.11$ and absolute magnitude $H_{r} = 3.6 \pm 0.1$; for an assumed albedo of $p_V = 12$% the object has a diameter of $\sim670$ km. Based on astrometric measurements from OSSOS and Pan-STARRS1, we find that 2015 RR$_{245}$ is securely trapped on ten-Myr timescales in the 9:2 mean-motion resonance with Neptune. It is the first TNO identified in this resonance. On hundred-Myr timescales, particles in 2015 RR$_{245}$-like orbits depart and sometimes return to the resonance, indicating that 2015 RR$_{245}$ likely forms part of the long-lived metastable population of distant TNOs that drift between resonance sticking and actively scattering via gravitational encounters with Neptune. The discovery of a 9:2 TNO stresses the role of resonances in the long-term evolution of objects in the scattering disk, and reinforces the view that distant resonances are heavily populated in the current Solar System. This object further motivates detailed modelling of the transient sticking population.

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 report the discovery, tracking and detection circumstances for 85 trans-Neptunian objects (TNOs) from the first 42 deg$^{2}$ of the Outer Solar System Origins Survey (OSSOS). This ongoing $r$-band Solar System survey uses the 0.9 deg$^{2}$ field-of-view MegaPrime camera on the 3.6 m Canada-France-Hawaii Telescope. Our orbital elements for these TNOs are precise to a fractional semi-major axis uncertainty $<0.1\%$. We achieve this precision in just two oppositions, as compared to the normal 3-5 oppositions, via a dense observing cadence and innovative astrometric technique. These discoveries are free of ephemeris bias, a first for large trans-Neptunian surveys. We also provide the necessary information to enable models of TNO orbital distributions to be tested against our TNO sample. We confirm the existence of a cold "kernel" of objects within the main cold classical Kuiper belt, and infer the existence of an extension of the "stirred" cold classical Kuiper belt to at least several AU beyond the 2:1 mean motion resonance with Neptune. We find that the population model of Petit et al. (2011) remains a plausible representation of the Kuiper belt. The full survey, to be completed in 2017, will provide an exquisitely characterized sample of important resonant TNO populations, ideal for testing models of giant planet migration during the early history of the Solar System.

The first two observational sky "blocks" of the Outer Solar System Origins Survey (OSSOS) have significantly increased the number of well-characterized observed trans-Neptunian objects (TNOs) in Neptune's mean motion resonances. We describe the 31 securely resonant TNOs detected by OSSOS so far, and we use them to independently verify the resonant population models from the Canada-France Ecliptic Plane Survey (CFEPS; Gladman et al. 2012), with which we find broad agreement. We confirm that the 5:2 resonance is more populated than models of the outer Solar System's dynamical history predict; our minimum population estimate shows that the high eccentricity (e>0.35) portion of the resonance is at least as populous as the 2:1 and possibly as populated as the 3:2 resonance. One OSSOS block was well-suited to detecting objects trapped at low libration amplitudes in Neptune's 3:2 resonance, a population of interest in testing the origins of resonant TNOs. We detected three 3:2 objects with libration amplitudes below the cutoff modeled by CFEPS; OSSOS thus offers new constraints on this distribution. The OSSOS detections confirm that the 2:1 resonance has a dynamically colder inclination distribution than either the 3:2 or 5:2 resonances. Using the combined OSSOS and CFEPS 2:1 detections, we constrain the fraction of 2:1 objects in the symmetric mode of libration to be 0.2-0.85; we also constrain the fraction of leading vs. trailing asymmetric librators, which has been theoretically predicted to vary depending on Neptune's migration history, to be 0.05-0.8. Future OSSOS blocks will improve these constraints.

The Gaia mission started its regular observing program in the summer of 2014, and since then it is regularly obtaining observations of asteroids. This paper draws the outline of the data processing for Solar System objects, and in particular on the daily "short-term" processing, from the on-board data acquisition to the ground-based processing. We illustrate the tools developed to compute predictions of asteroid observations, we discuss the procedures implemented by the daily processing, and we illustrate some tests and validations of the processing of the asteroid observations. Our findings are overall consistent with the expectations concerning the performances of Gaia and the effectiveness of the developed software for data reduction.

We measure the absolute magnitude, $H$, distribution, $dN(H) \propto 10^{\alpha H}$ of the scattering Trans-Neptunian Objects (TNOs) as a proxy for their size-frequency distribution. We show that the H-distribution of the scattering TNOs is not consistent with a single-slope distribution, but must transition around $H_g \sim 9$ to either a knee with a shallow slope or to a divot, which is a differential drop followed by second exponential distribution. Our analysis is based on a sample of 22 scattering TNOs drawn from three different TNO surveys, the Canada-France Ecliptic Plane Survey (CFEPS, Petit et al. 2011), Alexandersen et al. (2014), and the Outer Solar System Origins Survey (OSSOS, Bannister et al. 2016), all of which provide well characterized detection thresholds, combined with a cosmogonic model for the formation of the scattering TNO population. Our measured absolute magnitude distribution result is independent of the choice of cosmogonic model. Based on our analysis, we estimate that number of scattering TNOs is (2.4-8.3)$\times 10^5$ for $H_r < 12$. A divot $H$-distribution is seen in a variety of formation scenarios and may explain several puzzles in Kuiper Belt science. We find that a divot $H$-distribution simultaneously explains the observed scattering TNO, Neptune Trojan, Plutino, and Centaur $H$-distributions while simultaneously predicting a large enough scattering TNO population to act as the sole supply of the Jupiter-Family Comets.

Comets and chondrites show non-monotonic behaviour of their Deuterium to Hydrogen (D/H) ratio as a function of their formation location from the Sun. This is difficult to explain with a classical protoplanetary disk model that has a decreasing temperature structure with radius from the Sun. We want to understand if a protoplanetary disc with a dead zone, a region of zero or low turbulence, can explain the measured D/H values in comets and chondrites. We use time snapshots of a vertically layered disk model with turbulent surface layers and a dead zone at the midplane. The disc has a non-monotonic temperature structure due to increased heating from self-gravity in the outer parts of the dead zone. We couple this to a D/H ratio evolution model in order to quantify the effect of such thermal profiles on D/H enrichment in the nebula. We find that the local temperature peak in the disk can explain the diversity in the D/H ratios of different chondritic families. This disk temperature profile leads to a non-monotonic D/H enrichment evolution, allowing these families to acquire their different D/H values while forming in close proximity. The formation order we infer for these families is compatible with that inferred from their water abundances. However, we find that even for very young disks, the thermal profile reversal is too close to the Sun to be relevant for comets.

The Canada-France Ecliptic Plane Survey discovered four trans-Neptunian objects with semi-major axes near the 5:1 resonance, revealing a large and previously undetected intrinsic population. Three of these objects are currently resonant with Neptune, and the fourth is consistent with being an object that escaped the resonance at some point in the past. The non-resonant object may be representative of a detached population that is stable at slightly lower semi-major axes than the 5:1 resonance. We generated clones of these objects by resampling the astrometric uncertainty and examined their behavior over a 4.5 Gyr numerical simulation. The majority of the clones of the three resonant objects (>90%) spend a total of 10^7 years in resonance during their 4.5 Gyr integrations; most clones experience multiple periods of resonance capture. Our dynamical integrations reveal an exchange between the 5:1 resonance, the scattering objects, and other large semi-major axis resonances, especially the 4:1, 6:1, and 7:1. The multiple capture events and relatively short resonance lifetimes after capture suggest that these objects are captured scattering objects that stick in the 5:1 resonance. These 5:1 resonators may be representative of a temporary population, requiring regular contributions from a source population. We examined the dynamical characteristics (inclination, eccentricity, resonant island, libration amplitude) of the detected objects and their clones in order to provide an empirical model of the orbit structure of the 5:1 resonance. This resonance is dynamically hot and includes primarily symmetric librators. Given our orbit model, the intrinsic population necessary for the detection of these three objects in the 5:1 resonance is 1900(+3300 -1400, 95% confidence) with H_g< 8 and e > 0.5.

The Trans-Neptunian Objects (TNOs) may preserve evidence of planet building in their orbital and size-distributions. While all populations show steep size-distributions for large objects, recently relative deficit of Neptunian Trojans and scattering objects with diameters $D<100\,\mathrm{km}$ were detected. We have investigated this deficit with a 32 square degree survey, detecting 77 TNOs to a limiting $r$-band magnitude of 24.6. Our Plutinos sample (18 objects in 3:2 mean motion resonance with Neptune) also shows a deficit of $D<100\,\mathrm{km}$ objects. We reject a single power-law size-distribution and find that the Plutinos favour a divot. The Plutinos are thus added the list of populations with a deficit of $D<100\,\mathrm{km}$ objects. The fact that three independent samples of three different populations show this trend suggests that it is a real feature, possibly shared by all hot TNO populations as a remnant of "born big" planetesimal formation processes. We surmise the existence of $9000\pm3000$ Plutinos with $H_r\leq8.66$ and $37000^{+12000}_{-10000}$ Plutinos with $H_r\leq10.0$. Our survey also discovered one temporary Uranian Trojan, one temporary Neptunian Trojans and one stable Neptunian Trojan, from which we derive populations of $110^{+500}_{-100}$, $210^{+900}_{-200}$ and $150^{+600}_{-140}$ , respectively, with $H_r\leq10.0$. The Neptunian Trojans are thus less numerous than the main belt asteroids, which has over 700 asteroids with $H_r\leq10.0$. With such numbers, the temporary Neptunian Trojans cannot be previously stable Trojans that happen to be escaping the resonance now; they must be captured from another reservoir. With three 3:1 and one 4:1 resonators, we add to the growing evidence that the high-order resonances are more populated than expected.

Several chapters in this book present various aspects of Michel Henon's scientific acheivements that spread over a large range of subjects, and yet managed to make deep contributions to most of them. The authors of these chapters make a much better job at demonstrating the big advancements that Michel Henon allowed in these fields than I could ever do. Here I rather present some facets of his personnality that most appealed to me. Michel Henon was a reserved person, almost shy, so it was not obvious for a young student to grasp the profoundness of his insight and what a marvelous advisor he could be. The two most prominent aspects of his mind, in my view, were his ability to simplify any scientific question to its core complexity, and to find the fun and amusing part in his everyday work, even in the tiniest details of his scientific investigations.

The formation mechanisms of the ice giants Uranus and Neptune, and the origin of their elemental and isotopic compositions, have long been debated. The density of solids in the outer protosolar nebula is too low to explain their formation, and spectroscopic observations show that both planets are highly enriched in carbon, very poor in nitrogen, and the ices from which they originally formed might had deuterium-to-hydrogen ratios lower than the predicted cometary value, unexplained properties observed in no other planets. Here we show that all these properties can be explained naturally if Uranus and Neptune both formed at the carbon monoxide iceline. Due to the diffusive redistribution of vapors, this outer region of the protosolar nebula intrinsically has enough surface density to form both planets from carbon-rich solids but nitrogen-depleted gas, in abundances consistent with their observed values. Water rich interiors originating mostly from transformed CO ices reconcile the D/H value of Uranus and Neptune's building blocks with the cometary value. Finally, Our scenario generalizes a well known hypothesis that Jupiter formed on an iceline (water snowline) for the two ice giants, and might be a first step towards generalizing this mechanism for other giant planets.

Remote sensing observations meet some limitations when used to study the bulk atmospheric composition of the giant planets of our solar system. A remarkable example of the superiority of in situ probe measurements is illustrated by the exploration of Jupiter, where key measurements such as the determination of the noble gases abundances and the precise measurement of the helium mixing ratio have only been made available through in situ measurements by the Galileo probe. This paper describes the main scientific goals to be addressed by the future in situ exploration of Saturn placing the Galileo probe exploration of Jupiter in a broader context and before the future probe exploration of the more remote ice giants. In situ exploration of Saturn's atmosphere addresses two broad themes that are discussed throughout this paper: first, the formation history of our solar system and second, the processes at play in planetary atmospheres. In this context, we detail the reasons why measurements of Saturn's bulk elemental and isotopic composition would place important constraints on the volatile reservoirs in the protosolar nebula. We also show that the in situ measurement of CO (or any other disequilibrium species that is depleted by reaction with water) in Saturn's upper troposphere would constrain its bulk O/H ratio. We highlight the key measurements required to distinguish competing theories to shed light on giant planet formation as a common process in planetary systems with potential applications to most extrasolar systems. In situ measurements of Saturn's stratospheric and tropospheric dynamics, chemistry and cloud-forming processes will provide access to phenomena unreachable to remote sensing studies. Different mission architectures are envisaged, which would benefit from strong international collaborations.

Amateur contributions to professional publications have increased exponentially over the last decades in the field of Planetary Astronomy. Here we review the different domains of the field in which collaborations between professional and amateur astronomers are effective and regularly lead to scientific publications. We discuss the instruments, detectors, softwares and methodologies typically used by amateur astronomers to collect the scientific data in the different domains of interest. Amateur contributions to the monitoring of planets and interplanetary matter, characterization of asteroids and comets, as well as the determination of the physical properties of Kuiper Belt Objects and exoplanets are discussed.

The C--to--O ratio is a crucial determinant of the chemical properties of planets. The recent observation of WASP 12b, a giant planet with a C/O value larger than that estimated for its host star, poses a conundrum for understanding the origin of this elemental ratio in any given planetary system. In this paper, we propose a mechanism for enhancing the value of C/O in the disk through the transport and distribution of volatiles. We construct a model that computes the abundances of major C and O bearing volatiles under the influence of gas drag, sublimation, vapor diffusion, condensation and coagulation in a multi--iceline 1+1D protoplanetary disk. We find a gradual depletion in water and carbon monoxide vapors inside the water's iceline with carbon monoxide depleting slower than water. This effect increases the gaseous C/O and decreases the C/H ratio in this region to values similar to those found in WASP 12b's day side atmosphere. Giant planets whose envelopes were accreted inside the water's iceline should then display C/O values larger than those of their parent stars, making them members of the class of so-called ``carbon-rich planets''.

Context. The observation of carbon-rich disks have motivated several studies questioning the influence of the C/O ratio on their gas phase composition in order to establish the connection between the metallicity of hot-Jupiters and that of their parent stars. Aims. We to propose a method that allows the characterization of the adopted C/O ratio in protoplanetary disks independently from the determination of the host star composition. Titanium and vanadium chemistries are investigated because they are strong optical absorbers and also because their oxides are known to be sensitive to the C/O ratio in some exoplanet atmospheres. Methods. We use a commercial package based on the Gibbs energy minimization technique to compute the titanium and vanadium equilibrium chemistries in protoplanetary disks for C/O ratios ranging from 0.05 to 10. Our calculations are performed for pressures ranging from 1e-6 to 1e-2 bar, and for temperatures ranging from 50 to 2000 K. Results. We find that the vanadium nitride/vanadium oxide and titanium hydride/titanium oxide gas phase ratios strongly depend on the C/O ratio in the hot parts of disks (T > 1000 K). Our calculations suggest that, in these regions, these ratios can be used as tracers of the C/O value in protoplanetary disks.

Trojan objects share a planet's orbit, never straying far from the triangular Lagrangian points, 60 degrees ahead of (L4) or behind (L5) the planet. We report the detection of a Uranian Trojan; in our numerical integrations, 2011 QF99 oscillates around the Uranian L4 Lagrange point for >70 kyr and remains co-orbital for ~1 Myr before becoming a Centaur. We constructed a Centaur model, supplied from the transneptunian region, to estimate temporary co-orbital capture frequency and duration (to factor of two accuracy), finding that at any time 0.4% and 2.8% of the population will be Uranian and Neptunian co-orbitals, respectively. The co-orbital fraction (~2.4%) among Centaurs in the IAU Minor Planet Centre database is thus as expected under transneptunian supply.

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

Via joint analysis of a calibrated telescopic survey, which found scattering Kuiper Belt objects, and models of their expected orbital distribution, we measure the form of the scattering object's size distribution. Ruling out a single power-law at greater than 99% confidence, we constrain the form of the size distribution and find that, surprisingly, our analysis favours a very sudden decrease (a divot) in the number distribution as diameters decrease below 100 km, with the number of smaller objects then rising again as expected via collisional equilibrium. Extrapolating at this collisional equilibrium slope produced enough kilometer-scale scattering objects to supply the nearby Jupiter-Family comets. Our interpretation is that this divot feature is a preserved relic of the size distribution made by planetesimal formation, now "frozen in" to portions of the Kuiper Belt sharing a "hot" orbital inclination distribution, explaining several puzzles in Kuiper Belt science. Additionally, we show that to match today's scattering-object inclination distribution, the supply source that was scattered outward must have already been vertically heated to of order 10 degrees.

We propose a scenario that explains the apparent nitrogen deficiency in comets in a way consistent with the fact that the surfaces of Pluto and Triton are dominated by nitrogen-rich ice. We use a statistical thermodynamic model to investigate the composition of the successive multiple guest clathrates that may have formed during the cooling of the primordial nebula from the most abundant volatiles present in the gas phase. These clathrates agglomerated with the other ices (pure condensates or stoichiometric hydrates) and formed the building blocks of comets. We report that molecular nitrogen is a poor clathrate former, when we consider a plausible gas phase composition of the primordial nebula. This implies that its trapping into cometesimals requires a low disk temperature ($\sim$20 K) in order to allow the formation of its pure condensate. We find that it is possible to explain the lack of molecular nitrogen in comets as a consequence of their postformation internal heating engendered by the decay of short-lived radiogenic nuclides. This scenario is found consistent with the presence of nitrogen-rich ice covers on Pluto and Triton. Our model predicts that comets should present xenon-to-water and krypton-to-water ratios close to solar xenon-to-oxygen and krypton-to-oxygen ratios, respectively. In contrast, the argon-to-water ratio is predicted to be depleted by a factor of $\sim$300 in comets compared to solar argon-to-oxygen, as a consequence of poor trapping efficiency and radiogenic heating.

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.

The low-inclination component of the Classical Kuiper Belt is host to a population of extremely widely-separated binaries. These systems are similar to other Trans-Neptunian binaries (TNBs) in that the primary and secondary components of each system are of roughly equal size. We have performed an astrometric monitoring campaign of a sample of seven wide-separation, long-period TNBs and present the first-ever well-characterized mutual orbits for each system. The sample contains the most eccentric (2006 CH69, e=0.9) and the most widely-separated, weakly bound (2001 QW322, a/Rh~0.22) binary minor planets known, and also contains the system with lowest-measured mass of any TNB (2000 CF105, M~1.85E17 kg). Four systems orbit in a prograde sense, and three in a retrograde sense. They have a different mutual inclination distribution compared to all other TNBs, preferring low mutual-inclination orbits. These systems have geometric r-band albedos in the range of 0.09-0.3, consistent with radiometric albedo estimates for larger solitary low-inclination Classical Kuiper Belt objects, and we limit the plausible distribution of albedos in this region of the Kuiper Belt. We find that gravitational collapse binary formation models produce a similar orbital distribution to that currently observed, which along with a confluence of other factors supports formation of the cold Classical Kuiper Belt in situ through relatively rapid gravitational collapse rather than slow hierarchical accretion. We show that these binary systems are sensitive to disruption via collisions, and their existence suggests that the size distribution of TNOs at small sizes remains relatively shallow.

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.

Hot temperature minerals have been detected in a large number of comets and were also identified in the samples of Comet Wild 2 that were returned by the Stardust mission. Meanwhile, observations of the distribution of hot minerals in young stellar systems suggest that these materials were produced in the inner part of the primordial nebula and have been transported outward in the formation zone of comets. We investigate the possibility that photophoresis provides a viable mechanism to transport high-temperature materials from the inner solar system to the regions in which the comets were forming. We use a grid of time-dependent disk models of the solar nebula to quantify the distance range at which hot minerals can be transported from the inner part of the disk toward its outer regions as a function of their size and density. The particles considered here are in the form of aggregates that presumably were assembled from hot mineral individual grains ranging down to submicron sizes and formed by condensation within the hottest portion of the solar nebula. Our particle-transport model includes the photophoresis, radiation pressure, and gas drag. Depending on the postulated disk parameters and the density of particles, 10-2 to 10-1 m aggregates can reach heliocentric distances up to 35 AU in the primordial nebula over very short timescales (no more than a few hundred thousand years). 10-3 m particles follow the same trajectory as the larger ones but their maximum migration distance does not exceed 26 AU and is reached at later epochs in the disks. On the other hand, 10-5 to 10-4 m aggregates are continuously pushed outward during the evolution of the solar nebula. Our simulations suggest that irrespective of the employed solar nebula model, photophoresis is a mechanism that can explain the presence of hot temperature minerals in the formation region of comets.

The D/H enrichment observed in Saturn's satellite Enceladus is remarkably similar to the values observed in the nearly-isotropic comets. Given the predicted strong variation of D/H with heliocentric distance in the solar nebula, this observation links the primordial source region of the nearly-isotropic comets with the formation location of Enceladus. That is, comets from the nearly-isotropic class were most likely fed into their current reservoir, the Oort cloud, from a source region near the formation location of Enceladus. Dynamical simulations of the formation of the Oort cloud indicate that Uranus and Neptune are, primarily, responsible for the delivery of material into the Oort cloud. In addition, Enceladus formed from material that condensed from the solar nebula near the location at which Saturn captured its gas envelope, most likely at or near Saturn's current location in the solar system. The coupling of these lines of evidence appears to require that Uranus and Neptune were, during the epoch of the formation of the Oort cloud, much closer to the current location of Saturn than they are currently. Such a configuration is consistent with the Nice model of the evolution of the outer solar system. Further measurements of the D/H enrichment in comets, particularly in ecliptic comets, will provide an excellent discriminator among various models of the formation of the outer solar system.

Favored theories of giant planet formation center around two main paradigms, namely the core accretion model and the gravitational instability model. These two formation scenarios support the hypothesis that the giant planet metallicities should be higher or equal to that of the parent star. Meanwhile, spectra of the transiting hot Jupiter HD189733b suggest that carbon and oxygen abundances range from depleted to enriched with respect to the star. Here, using a model describing the formation sequence and composition of planetesimals in the protoplanetary disk, we determine the range of volatile abundances in the envelope of HD189733b that is consistent with the 20--80 Earth-masses of heavy elements estimated to be present in the planet's envelope. We then compare the inferred carbon and oxygen abundances to those retrieved from spectroscopy and we find a range of supersolar values that directly fit both spectra and internal structure models. In some cases, we find that the apparent contradiction between the subsolar elemental abundances and the {mass of heavy elements predicted in HD189733b by internal structure models} can be explained by the presence of large amounts of carbon molecules in the form of polycyclic aromatic hydrocarbons and soots in the upper layers of the envelope, as suggested by recent photochemical models. A diagnostic test that would confirm the presence of these compounds in the envelope is the detection of acetylene.

We report here the discovery of an new near-equal mass Trans-Neptunian Binaries (TNBs) L5c02 and the and the putative detection of a second TNB (L4k12) among the year two and three detections of the Canada-France-Eclipic Plane Survey (CFEPS). These new binaries (internal designation L4k12 and L5c02) have moderate separations of 0.4" and 0.6" respectively. The follow-up observation confirmed the binarity of L5c02, but L4k12 are still lack of more followup observations. L4k12 has a heliocentric orbital inclination of ~ 35?degree, marking this system as having the highest heliocentric orbital inclination among known near-equal mass binaries. Both systems are members of the classical main Kuiper belt population. Based on the sample of objects searched we determine that the fraction of near-equal mass wide binaries with separations > 0.4" is 1.5% to 20% in the cold main classical Kuiper belt and, if our detection the binarity L4k12 holds, 3% to 43% in the hot main classical objects are binary. In this manuscript we describe our detection process, the sample of objects surveyed, our confirmation observations.

The difference between the measured atmospheric abundances of neon, argon, krypton and xenon for Venus, the Earth and Mars is striking. Because these abundances drop by at least two orders of magnitude as one moves outward from Venus to Mars, the study of the origin of this discrepancy is a key issue that must be explained if we are to fully understand the different delivery mechanisms of the volatiles accreted by the terrestrial planets. In this work, we aim to investigate whether it is possible to quantitatively explain the variation of the heavy noble gas abundances measured on Venus, the Earth and Mars, assuming that cometary bombardment was the main delivery mechanism of these noble gases to the terrestrial planets. To do so, we use recent dynamical simulations that allow the study of the impact fluxes of comets upon the terrestrial planets during the course of their formation and evolution. Assuming that the mass of noble gases delivered by comets is proportional to rate at which they collide with the terrestrial planets, we show that the krypton and xenon abundances in Venus and the Earth can be explained in a manner consistent with the hypothesis of cometary bombardment. In order to explain the krypton and xenon abundance differences between the Earth and Mars, we need to invoke the presence of large amounts of CO2-dominated clathrates in the Martian soil that would have efficiently sequestered these noble gases.

We have performed an ecliptic imaging survey of the Kuiper belt with our deepest and widest field achieving a limiting flux of m(g') = 26.4, with a sky coverage of 3.0 square-degrees. This is the largest coverage of any other Kuiper belt survey to this depth. We detect 72 objects, two of which have been previously observed. We have improved the Bayesian maximum likelihood fitting technique presented in Gladman et al. (1998) to account for calibration and sky density variations and have used this to determine the luminosity function of the Kuiper belt. Combining our detections with previous surveys, we find the luminosity function is well represented by a single power-law with slope alpha = 0.65 +/- 0.05 and an on ecliptic sky density of 1 object per square-degree brighter than m(R)=23.42 +/- 0.13. Assuming constant albedos, this slope suggests a differential size-distribution slope of 4.25 +/- 0.25, which is steeper than the Dohnanyi slope of 3.5 expected if the belt is in a state of collisional equilibrium. We find no evidence for a roll-over or knee in the luminosity function and reject such models brightward of m(R) ~ 24.6.

Protoplanetary disks start their lives with a dust free inner region where the temperatures are higher than the sublimation temperature of solids. As the star illuminates the innermost particles, which are immersed in gas at the sublimation edge, these particles are subject to a photophoretic force. We examine the motion of dust particles at the inner edge of protoplanetary disks due to photophoretic drag. We find that, in a laminar disk, photophoretic drag increases the size of the inner hole after accretion onto the central body has become subdued. This region within the hole becomes an optically transparent zone containing gas and large dusty particles (>>10 cm), but devoid of, or strongly depleted in, smaller dust aggregates. Photophoresis can clear the inner disk of dust out to 10 AU in less than 1 Myr. At late times, the edge reaches a stable equilibrium between inward drift and photophoretic outward drift, at a distance of some tens of AU. Eventually, the edge will move inwards again as the disk disperses, shifting the equilibrium position back from about 40 AU to below 30 AU in 1-2 Myr in the disk model. In a turbulent disk, diffusion can delay the clearing of a disk by photophoresis. Smaller and/or age-independent holes of radii of a few AU are also possible outcomes of turbulent diffusion counteracting photophoresis. This outward and then inward moving edge marks a region of high dust concentration. This density enhancement, and the efficient transport of particles from close to the star to large distances away, can explain features of comets such as high measured ratios of crystalline to amorphous silicates, and has a large number of other applications.

We present the data acquisition strategy and characterization procedures for the Canada-France Ecliptic Plane Survey (CFEPS), a sub-component of the Canada-France-Hawaii Telescope Legacy Survey. The survey began in early 2003 and as of summer 2005 has covered 430 square degrees of sky within a few degrees of the ecliptic. Moving objects beyond the orbit of Uranus are detected to a magnitude limit of $m_R$=23 -- 24 (depending on the image quality). To track as large a sample as possible and avoid introducing followup bias, we have developed a multi-epoch observing strategy that is spread over several years. We present the evolution of the uncertainties in ephemeris position and orbital elements as the objects progress through the epochs. We then present a small 10-object sample that was tracked in this manner as part of a preliminary survey starting a year before the main CFEPS project. We describe the CFEPS survey simulator, to be released in 2006, which allows theoretical models of the Kuiper Belt to be compared with the survey discoveries since CFEPS has a well-documented pointing history with characterized detection efficiencies as a function of magnitude and rate of motion on the sky. Using the pre-survey objects we illustrate the usage of the simulator in modeling the classical Kuiper Belt.

Using the cometary nucleus model developed by Espinasse et al. (1991), we calculate the thermodynamical evolution of Comet 9P/Tempel 1 over a period of 360 years. Starting from an initially amorphous cometary nucleus which incorporates an icy mixture of H2O and CO, we show that, at the time of Deep Impact collision, the crater is expected to form at depths where ice is in its crystalline form. Hence, the subsurface exposed to space should not be primordial. We also attempt an order-of-magnitude estimate of the heating and material ablation effects on the crater activity caused by the 370 Kg projectile released by the DI spacecraft. We thus show that heating effects play no role in the evolution of crater activity. We calculate that the CO production rate from the impacted region should be about 300-400 times higher from the crater resulting from the impact with a 35 m ablation than over the unperturbed nucleus in the immediate post-impact period. We also show that the H2O production rate is decreased by several orders of magnitude at the crater base just after ablation.

Given the large orbital separation and high satellite-to-primary mass ratio of all known Kuiper Belt Object (KBO) binaries, it is important to reassess their stability as bound pairs with respect to several disruptive mechanisms. Besides the classical shattering and dispersing of the secondary due to a high-velocity impact, we considered the possibility that the secondary is kicked off its orbit by a direct collision of a small impactor, or that it is gravitationally perturbed due to the close approach of a somewhat larger TNO. Depending on the values for the size/mass/separation of the binaries that we used, 2 or 3 of the 8 pairs can be dispersed in a timescale shorter than the age of the solar system in the current rarefied environment. A contemporary formation scenario could explain why we still observe these binaries, but no convincing mechanism has been proposed to date. The primordial formation scenarios, which seem to be the only viable ones, must be revised to increase the formation efficiency in order to account for this high dispersal rate. Objects like the large-separatioKBO binary n2001 QW$_{322}$ must have been initially an order of magnitude more numerous. If the KBO binaries are indeed primordial, then we show that the mass depletion of the Kuiper belt cannot result from collisional grinding, but must rather be due to dynamical ejection.

By telescopic tracking, we have established that the orbit of the trans-neptunian object (2000 CR$_{105}$) has a perihelion of $\simeq$44 AU, and is thus outside the domain controlled by strong gravitational close encounters with Neptune. Because this object is on a very large, eccentric orbit (with semimajor axis $a\simeq$216 AU and eccentricity $e\simeq$0.8) this object must have been placed on this orbit by a gravitational perturbation which is {\it not} direct gravitational scattering off of any of the giant planets (on their current orbits). The existence of this object may thus have profound cosmogonic implications for our understanding of the formation of the outer Solar System. We discuss some viable scenarios which could have produced it, including long-term diffusive chaos and scattering off of other massive bodies in the outer Solar System. This discovery implies that there must be a large population of trans-neptunian objects in an `extended scattered disk' with perihelia above the previously-discussed 38 AU boundary.