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.

The exploration of the binary Pluto-Charon and its small satellites during the New Horizons flyby in 2015 revealed not only widespread geologic and compositional diversity across Pluto, but surprising complexity, a wide range of surface unit ages, evidence for widespread activity stretching across billion of years to the near-present, as well as numerous atmospheric puzzles, and strong atmospheric coupling with its surface. New Horizons also found an unexpected diversity of landforms on its binary companion, Charon. Pluto's four small satellites yielded surprises as well, including their unexpected rapid and high obliquity rotation states, high albedos, and diverse densities. Here we briefly review the findings made by New Horizons and the case for a follow up mission to investigate the Pluto system in more detail. As the next step in the exploration of this spectacular planet-satellite system, we recommend an orbiter to study it in considerably more detail, with new types of instrumentation, and to observe its changes with time. We further call for the in-depth study of Pluto orbiter missions as a precursor to the 2023 Planetary Science Decadal Survey.

It is often claimed that asteroids' sharing of orbits is the reason they were re-classified from planets to non-planets. A critical review of the literature from the 19th Century to the present shows this is factually incorrect. The literature shows the term asteroid was broadly recognized as a subset of planet for 150 years. On-going discovery of asteroids resulted in a de facto stretching of the concept of planet to include the ever-smaller bodies. Scientists found utility in this taxonomic identification as it provided categories needed to argue for the leading hypothesis of planet formation, Laplace's nebular hypothesis. In the 1950s, developments in planet formation theory found it no longer useful to maintain taxonomic identification between asteroids and planets, Ceres being the primary exception. At approximately the same time, there was a flood of publications on the geophysical nature of asteroids showing them to be geophysically different than the large planets. This is when the terminology in asteroid publications calling them planets abruptly plunged from a high level of usage where it had hovered during the period 1801 - 1957 to a low level that held constant thereafter. This marks the point where the community effectively formed consensus that asteroids should be taxonomically distinct from planets. The evidence demonstrates this consensus formed on the basis of geophysical differences between asteroids and planets, not the sharing of orbits. We suggest attempts to build consensus around planetary taxonomy not rely on the non-scientific process of voting, but rather through precedent set in scientific literature and discourse, by which perspectives evolve with additional observations and information, just as they did in the case of asteroids.

NASA's New Horizons spacecraft will conduct a close flyby of the cold classical Kuiper Belt Object (KBO) designated (486958) 2014 MU69 on January 1, 2019. At a heliocentric distance of 44 AU, "MU69" will be the most distant object ever visited by a spacecraft. To enable this flyby, we have developed an extremely high precision orbit fitting and uncertainty processing pipeline, making maximal use of the Hubble Space Telescope's Wide Field Camera 3 (WFC3) and pre-release versions of the ESA Gaia Data Release 2 (DR2) catalog. This pipeline also enabled successful predictions of a stellar occultation by MU69 in July 2017. We describe how we process the WFC3 images to match the Gaia DR2 catalog, extract positional uncertainties for this extremely faint target (typically 140 photons per WFC3 exposure), and translate those uncertainties into probability distribution functions for MU69 at any given time. We also describe how we use these uncertainties to guide New Horizons, plan stellar occultions of MU69, and derive MU69's orbital evolution and long-term stability.

Pluto's landscape is shaped by the cycles of the volatile ices covering its surface. In particular, the Sputnik Planitia (SP) ice sheet displays a large diversity of terrains, with bright and dark plains, pits, topographic depressions and evidences of recent and past glacial flows. Outside SP, New Horizons also revealed numerous N2 ice deposits, in Tombaugh Regio and at mid-northern latitudes. These observations suggest a complex history involving volatile and glacial processes on different timescales. We present numerical simulations of volatile transport on Pluto performed with a model able to simulate the N2 cycle over millions of years (Myrs), taking into account the changes of obliquity and orbital parameters as experienced by Pluto. Results show that over one obliquity cycle, the latitudes of SP between 25{\deg}S-30{\deg}N are dominated by N2 condensation, while the latitudes between 30-50{\deg}N are dominated by N2 sublimation. We find that a net amount of 1 km of ice has sublimed at the northern edge of SP during the last 2 Myrs. By comparing these results with the observed geology of SP, we can relate the formation of the pits and the brightness of the ice to the ice flux occurring at the annual timescale, while the glacial flows at its eastern edge and the erosion of the water ice mountains all around the ice sheet are related to the astronomical timescale. We also perform simulations with a glacial flow scheme which shows that SP is currently at its minimum extent. We also explore the stability of N2 ice outside SP. Results show that it is not stable at the poles but rather in the equatorial regions, in particular in depressions, where thick deposits may persist over tens of Myrs, before being trapped in SP. Finally, another key result is that the minimum and maximum surface pressures obtained over the simulated Myrs remain in the range of mm-Pa and Pa, respectively.

On July 14th 2015, NASA's New Horizons mission gave us an unprecedented detailed view of the Pluto system. The complex compositional diversity of Pluto's encounter hemisphere was revealed by the Ralph/LEISA infrared spectrometer on board of New Horizons. We present compositional maps of Pluto defining the spatial distribution of the abundance and textural properties of the volatiles methane and nitrogen ices and non-volatiles water ice and tholin. These results are obtained by applying a pixel-by-pixel Hapke radiative transfer model to the LEISA scans. Our analysis focuses mainly on the large scale latitudinal variations of methane and nitrogen ices and aims at setting observational constraints to volatile transport models. Specifically, we find three latitudinal bands: the first, enriched in methane, extends from the pole to 55deg N, the second dominated by nitrogen, continues south to 35deg N, and the third, composed again mainly of methane, reaches 20deg N. We demonstrate that the distribution of volatiles across these surface units can be explained by differences in insolation over the past few decades. The latitudinal pattern is broken by Sputnik Planitia, a large reservoir of volatiles, with nitrogen playing the most important role. The physical properties of methane and nitrogen in this region are suggestive of the presence of a cold trap or possible volatile stratification. Furthermore our modeling results point to a possible sublimation transport of nitrogen from the northwest edge of Sputnik Planitia toward the south.

Observations of the Pluto-Charon system, acquired with the ALMA interferometer on June 12-13, 2015, have yielded a detection of the CO(3-2) and HCN(4-3) rotational transitions from Pluto, providing a strong confirmation of the presence of CO, and the first observation of HCN, in Pluto's atmosphere. The CO and HCN lines probe Pluto's atmosphere up to ~450 km and ~900 km altitude, respectively. The CO detection yields (i) a much improved determination of the CO mole fraction, as 515+/-40 ppm for a 12 ubar surface pressure (ii) clear evidence for a well-marked temperature decrease (i.e., mesosphere) above the 30-50 km stratopause and a best-determined temperature of 70+/-2 K at 300 km, in agreement with recent inferences from New Horizons / Alice solar occultation data. The HCN line shape implies a high abundance of this species in the upper atmosphere, with a mole fraction >1.5x10-5 above 450 km and a value of 4x10-5 near 800 km. The large HCN abundance and the cold upper atmosphere imply supersaturation of HCN to a degree (7-8 orders of magnitude) hitherto unseen in planetary atmospheres, probably due to the slow kinetics of condensation at the low pressure and temperature conditions of Pluto's upper atmosphere. HCN is also present in the bottom ~100 km of the atmosphere, with a 10-8 - 10-7 mole fraction; this implies either HCN saturation or undersaturation there, depending on the precise stratopause temperature. The HCN column is (1.6+/-0.4)x10^14 cm-2, suggesting a surface-referred net production rate of ~2x10^7 cm-2s-1. Although HCN rotational line cooling affects Pluto's atmosphere heat budget, the amounts determined in this study are insufficient to explain the well-marked mesosphere and upper atmosphere's ~70 K temperature. We finally report an upper limit on the HC3N column density (< 2x10^13 cm-2) and on the HC15N / HC14N ratio (< 1/125).

The observation of flat rotation curves in asymptotic region of galaxies implies that the dark matter density profile decreases exponentially with the gravitational potential in this region. It is curious that this behavior is identical to what one would have from a Boltzmann gas in thermal equilibrium. This suggests that an {\it effective} Boltzmann gas description of dark matter is possible for the halo, which follows even though no requirement of thermal equilibrium is assumed for the dark matter particles. We examine some qualitative and quantitative consequences of such an effective description. This is done by studying a simple model that takes into account gravitational self-interactions of dark matter and its gravitational interactions with the baryonic core. Solutions to the dynamics are determined by two parameters, one of which is the ratio of the dark matter mass to the effective temperature of the gas. Crude estimates can be made for these parameters based on observations.

After deducing the density profiles and gravitational potential functions of eight galaxies from the rotation velocity data from THINGS, we find that the density decreases exponentially with the potential in substantial regions of the haloes. Such behavior is in agreement with that of a single-component isothermal Boltzmann gas, and suggests that an effective description in terms of a Boltzmann gas is possible for dark matter in these regions. This could be an indication that dark matter self-interactions are sufficient in strength and number to lead to thermal equilibrium in these regions. We write down the dynamics and boundary conditions for a Boltzmann gas description and examine some of its qualitative and quantitative consequences. Solutions to the dynamical system are determined by three dimensionfull parameters, and provide reasonable fits to the rotational velocity data in the regions where the Boltzmann-like behavior was found. Unlike in the usual approach to curve fitting, we do not assume a specific form for the dark matter density profile and we do not require a detailed knowledge of the baryonic content of the galaxy.

We present new optical emission-line images of the young SNR 1E 0102-7219 (E0102) in the SMC obtained with the HST Advanced Camera for Surveys (ACS). E0102 is a member of the oxygen-rich class of SNRs showing strong oxygen, neon , and other metal-line emissions in its optical and X-ray spectra, and an absence of H and He. The progenitor of E0102 may have been a Wolf-Rayet star that underwent considerable mass loss prior to exploding as a Type Ib/c or IIL/b SN. The ejecta in this SNR are fast-moving (V > 1000 km/s) and emit as they are compressed and heated in the reverse shock. In 2003, we obtained optical [O III], H-alpha, and continuum images with the ACS Wide Field Camera. The [O III] image captures the full velocity range of the ejecta, and shows considerable high-velocity emission projected in the middle of the SNR that was Doppler-shifted out of the narrow F502N bandpass of a previous Wide Field and Planetary Camera 2 image from 1995. Using these two epochs separated by ~8.5 years, we measure the transverse expansion of the ejecta around the outer rim in this SNR for the first time at visible wavelengths. From proper-motion measurements of 12 ejecta filaments, we estimate a mean expansion velocity for the bright ejecta of ~2000 km/s and an inferred kinematic age for the SNR of \~2050 +/- 600 years. The age we derive from HST data is about twice that inferred by Hughes et al.(2000) from X-ray data, though our 1-sigma error bars overlap. Our proper-motion age is consistent with an independent optical kinematic age derived by Eriksen et al.(2003) using spatially resolved [O III] radial-velocity data. We derive an expansion center that lies very close to X-ray and radio hotspots, which could indicate the presence of a compact remnant (neutron star or black hole).