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.

Centaurs are minor planets thought to have originated in the outer Solar System region known as the Kuiper Belt. Active Centaurs enigmatically display comet-like features (e.g., tails, comae) even though they orbit in the gas giant region where it is too cold for water to readily sublimate. Only 18 active Centaurs have been identified since 1927 and, consequently, the underlying activity mechanism(s) have remained largely unknown up to this point. Here we report the discovery of activity emanating from Centaur 2014 OG392, based on archival images we uncovered plus our own new observational evidence acquired with the Dark Energy Camera (Cerro Tololo Inter-American Observatory Blanco 4 m telescope), the Inamori-Magellan Areal Camera & Spectrograph (Las Campanas Observatory 6.5 m Walter Baade Telescope) and the Large Monolithic Imager (Lowell Observatory 4.3 m Discovery Channel Telescope). We detect a coma as far as 400,000 km from 2014 OG392, and our novel analysis of sublimation processes and dynamical lifetime suggest carbon dioxide and/or ammonia are the most likely candidates for causing activity on this and other active Centaurs. We find 2014 OG392 is optically red, but CO2 and NH3 are spectrally neutral in this wavelength regime so the reddening agent is as yet unidentified.

Can multicellular life be distinguished from single cellular life on an exoplanet? We hypothesize that abundant upright photosynthetic multicellular life (trees) will cast shadows at high sun angles that will distinguish them from single cellular life and test this using Earth as an exoplanet. We first test the concept using Unmanned Arial Vehicles (UAVs) at a replica moon landing site near Flagstaff, Arizona and show trees have both a distinctive reflectance signature (red edge) and geometric signature (shadows at high sun angles) that can distinguish them from replica moon craters. Next, we calculate reflectance signatures for Earth at several phase angles with POLDER (Polarization and Directionality of Earth's reflectance) satellite directional reflectance measurements and then reduce Earth to a single pixel. We compare Earth to other planetary bodies (Mars, the Moon, Venus, and Uranus) and hypothesize that Earths directional reflectance will be between strongly backscattering rocky bodies with no weathering (like Mars and the Moon) and cloudy bodies with more isotropic scattering (like Venus and Uranus). Our modelling results put Earth in line with strongly backscattering Mars, while our empirical results put Earth in line with more isotropic scattering Venus. We identify potential weaknesses in both the modeled and empirical results and suggest additional steps to determine whether this technique could distinguish upright multicellular life on exoplanets.

We present here the c-o colors for identified Flora, Vesta, Nysa-Polana, Themis, and Koronis family members within the historic data set (2015-2018) of the Asteroid Terrestrial-impact Last Alert System (ATLAS). The Themis and Koronis families are known to be relatively pure C- and S-type Bus-DeMeo taxonomic families, respectively, and the extracted color data from the ATLAS broadband c- and o-filters of these two families is used to demonstrate that the ATLAS c-o color is a sufficient parameter to distinguish between the C- and S-type taxonomies. The Vesta and Nysa-Polana families are known to display a mixture of taxonomies possibly due to Vesta's differentiated parent body origin and Nysa-Polana actually consisting of two nested families with differing taxonomies. Our data show that the Flora family also displays a large degree of taxonomic mixing and the data reveal a substantial H-magnitude dependence on color. We propose and exclude several interpretations for the observed taxonomic mix. Additionally, we extract rotation periods of all of the targets reported here and find good agreement with targets that have previously reported periods.

The {\it Transiting Exoplanet Survey Satellite} (TESS) searches for planets transiting bright and nearby stars using high-cadence, large-scale photometric observations. Full Frame Images provided by the TESS mission include large number of serendipitously observed main-belt asteroids. Due to the cadence of the published Full Frame Images we are sensitive to periods as long as of order tens of days, a region of phase space that is generally not accessible through traditional observing. This work represents a much less biased measurement of the period distribution in this period range. We have derived rotation periods for 300~main-belt asteroids and have partial lightcurves for a further 7277 asteroids, including 43 with periods $P > 100$ h; this large number of slow rotators is predicted by theory. Of these slow rotators we find none requiring significant internal strength to resist rotational reshaping. We find our derived rotation periods to be in excellent agreement with results in the Lightcurve Database for the 55~targets that overlap. Over the nominal two-year lifetime of the mission, we expect the detection of around 85,000 unique asteroids with rotation period solutions for around 6000 asteroids. We project that the systematic analysis of the entire TESS data set will increase the number of known slow-rotating asteroids (period > 100~h) by a factor of 10. Comparing our new period determinations with previous measurements in the literature, we find that the rotation period of asteroid (2320) Blarney has decreased by at least 20\% over the past decade, potentially due to surface activity or subcatastrophic collisions.

The Mission Accessible Near-Earth Object Survey (MANOS) aims to observe and characterize small (mean absolute magnitude H ~ 25 mag) Near-Earth Objects (NEOs) that are accessible by spacecraft (mean $\Delta v$ ~ 5.7 km/s) and that make close approaches with the Earth (mean Minimum Orbital Intersection Distance MOID ~ 0.03 AU). We present here the first results of the MANOS visible spectroscopic survey. The spectra were obtained from August 2013 to March 2018 at Lowell Observatory's Discovery Channel 4.3 meter telescope, and both Gemini North and South facilities. In total, 210 NEOs have been observed and taxonomically classified. Our taxonomic distribution shows significant variations with respect to surveys of larger objects. We suspect these to be due to a dependence of Main Belt source regions on object size. Compared to previous surveys of larger objects (Binzel et al. 2019, 2004; Perna et al. 2018), we report a lower fraction of S+Q-complex asteroids of 43.8 $\pm$ 4.6%. We associate this decrease with a lack of Phocaea family members at very small size. We also report higher fractions of X-complex and A-type asteroids of 23.8 $\pm$ 3.3% and 3.8 $\pm$ 1.3% respectively due to an increase of Hungaria family objects at small size. We find a strong correlation between the Q/S ratio and perihelion distance. We suggest this correlation is due to planetary close encounters with Venus playing a major role in turning asteroids from S to Q-type. This hypothesis is supported by a similar correlation between the Q/S ratio and Venus MOID.

Through the lens of the LSST Science Collaborations' experience, this paper advocates for new and improved ways to fund large, complex collaborations at the interface of data science and astrophysics as they work in preparation for and on peta-scale, complex surveys, of which LSST is a prime example. We advocate for the establishment of programs to support both research and infrastructure development that enables innovative collaborative research on such scales.

We present Searches After Gravitational-waves Using ARizona Observatories (SAGUARO), a comprehensive effort dedicated to the discovery and characterization of optical counterparts to gravitational wave (GW) events. SAGUARO utilizes ground-based facilities ranging from 1.5m to 10m in diameter, located primarily in the Northern Hemisphere. We provide an overview of SAGUARO's telescopic resources, pipeline for transient detection, and database for candidate visualization. We describe SAGUARO's discovery component, which utilizes the $5$~deg$^2$ field-of-view optical imager on the Mt. Lemmon 1.5m telescope, reaching limits of $\approx 21.3$~AB mag while rapidly tiling large areas. We also describe the follow-up component of SAGUARO, used for rapid vetting and monitoring of optical candidates. With the onset of Advanced LIGO/Virgo's third observing run, we present results from the first three SAGUARO searches following the GW events S190408an, S190425z and S190426c, which serve as a valuable proof-of-concept of SAGUARO. We triggered and searched 15, 60 and 60 deg$^{2}$ respectively, 17.6, 1.4 and 41.8 hrs after the initial GW alerts. We covered 7.8, 3.0 and 5.1\% of the total probability within the GW event localizations, reaching 3$\sigma$ limits of 19.8, 21.3 and 20.8 AB mag, respectively. Although no viable counterparts associated with these events were found, we recovered 6 known transients and ruled out 5 potential candidates. We also present Large Binocular Telescope spectroscopy of PS19eq/SN2019ebq, a promising kilonova candidate that was later determined to be a supernova. With the ability to tile large areas and conduct detailed follow-up, SAGUARO represents a significant addition to GW counterpart searches.

Thermal infrared observations are the most effective way to measure asteroid diameter and albedo for a large number of near-Earth objects. Major surveys like NEOWISE, NEOSurvey, ExploreNEOs, and NEOLegacy find a small fraction of high albedo objects that do not have clear analogs in the current meteorite population. About 8% of Spitzer-observed near-Earth objects have nominal albedo solutions greater than 0.5. This may be a result of lightcurve variability leading to an incorrect estimate of diameter or inaccurate absolute visual magnitudes. For a sample of 23 high albedo NEOs we do not find that their shapes are significantly different from the McNeill et al. (2019) near-Earth object shape distribution. We performed a Monte Carlo analysis on 1505 near-Earth objects observed by Spitzer, sampling the visible and thermal fluxes of all targets to determine the likelihood of obtaining a high albedo erroneously. Implementing the McNeill shape distribution for near-Earth objects, we provide an upper-limit on the geometric albedo of 0.5+/-0.1 for the near-Earth population.

A combination band due to a mechanism whereby a photon excites two or more vibrational modes ({\it e.g.} a bend and a stretch) of an individual molecule is commonly seen in laboratory and astronomical spectroscopy. Here, we present evidence of a much less commonly seen combination band $-$ one where a photon simultaneously excites two adjacent molecules in an ice. In particular, we present near-infrared spectra of laboratory CO/N$_2$ ice samples where we identify a band at 4467.5 cm$^{-1}$ (2.239 $\mu$m) that results from single photons exciting adjacent pairs of CO and N$_2$ molecules. We also present a near-infrared spectrum of Neptune's largest satellite Triton taken with the Gemini-South 8.1 meter telescope and the Immersion Grating Infrared Spectrograph (IGRINS) that shows this 4467.5 cm$^{-1}$ (2.239 $\mu$m) CO-N$_2$ combination band. The existence of the band in a spectrum of Triton indicates that CO and N$_2$ molecules are intimately mixed in the ice rather than existing as separate regions of pure CO and pure N$_2$ deposits. Our finding is important because CO and N$_2$ are the most volatile species on Triton and so dominate seasonal volatile transport across its surface. Our result will place constraints on the interaction between the surface and atmosphere of Triton. 1

Though pairs of dynamically associated asteroids in the Main Belt have been identified and studied for over a decade, very few pair systems have been identified in the near-Earth asteroid population. We present data and analysis that supports the existence of two genetically related pairs in near-Earth space. The members of the individual systems, 2015 EE7 -- 2015 FP124 and 2017 SN16 -- 2018 RY7, are found to be of the same spectral taxonomic class, and both pairs are interpreted to have volatile-poor compositions. In conjunction with dynamical arguments, this suggests that these two systems formed via YORP spin-up and/or dissociation of a binary precursor. Backwards orbital integrations suggest a separation age of <10 kyr for the pair 2017 SN16 -- 2018 RY7, making these objects amongst the youngest multiple asteroid systems known to date. A unique separation age was not realized for 2015 EE7 -- 2015 FP124 due to large uncertainties associated with these objects' orbits. Determining the ages of such young pairs is of great value for testing models of space weathering and asteroid spin-state evolution. As the NEO catalog continues to grow with current and future discovery surveys, it is expected that more NEO pairs will be found, thus providing an ideal laboratory for studying time dependent evolutionary processes that are relevant to asteroids throughout the Solar System.

We present archival observations demonstrating that main belt asteroid (6478) Gault has an extensive history of comet-like activity. Outbursts have taken place during multiple epochs since 2013 and at distances extending as far as 2.68 au, nearly aphelion. (6478) Gault is a member of the predominately S-type (i.e., volatile-poor) Phocaea family; no other main belt object of this type has ever shown more than a single activity outburst. Furthermore, our data suggest that this is the longest duration of activity caused by a body spinning near the rotational breakup barrier. If activity is indeed unrelated to volatiles, as appears to be the case, (6478) Gault represents a new class of object, perpetually active due to rotational spin-up.

Unexpected dynamic phenomena have surprised solar system observers in the past and have led to important discoveries about solar system workings. Observations at the initial stages of these events provide crucial information on the physical processes at work. We advocate for long-term/permanent programs on ground-based and space-based telescopes of all sizes - including Extremely Large Telescopes (ELTs) - to conduct observations of high-priority dynamic phenomena, based on a predefined set of triggering conditions. These programs will ensure that the best initial dataset of the triggering event are taken; separate additional observing programs will be required to study the temporal evolution of these phenomena. While not a comprehensive list, the following are notional examples of phenomena that are rare, that cannot be anticipated, and that provide high-impact advances to our understandings of planetary processes. Examples include: new cryovolcanic eruptions or plumes on ocean worlds; impacts on Jupiter, Saturn, Uranus, or Neptune; extreme eruptions on Io; convective superstorms on Saturn, Uranus, or Neptune; collisions within the asteroid belt or other small-body populations; discovery of an interstellar object passing through our solar system (e.g. 'Oumuamua); and responses of planetary atmospheres to major solar flares or coronal mass ejections.

As part of our multi-observatory, multi-filter campaign, we present \rmi color observations of 82 Near-Earth Objects (NEOs) obtained with the RATIR instrument on the 1.5m robotic telescope at the San Pedro Martir's National Observatory in Mexico. Our project is particularly focused on rapid response observations of small ($\lesssim 850$ m) NEOs. The rapid response and the use of spectrophotometry allows us to constrain the taxonomic classification of NEOs with high efficiency. Here we present the methodology of our observations and our result, suggesting that the ratio of C-type to S-type asteroids in a size range of $\sim$30-850m is 1.1, which is in accordance with our previous results. We also find that 10$\%$ of all NEOs in our sample are neither C- nor S-type asteroids

We present here multi-band photometry for over 2000 Main-belt asteroids. For each target we report the probabilistic taxonomy using the measured V-R and V-I colors in combination with a machine-learning generated decision surface in color-color space. Through this method we classify >85% of our targets as one the four main Bus-DeMeo complexes: S-, C-, X-, or D-type. Roughly one third of our targets have a known associated dynamic family with 69 families represented in our data. Within uncertainty our results show no discernible difference in taxonomic distribution between family members and non-family members. Nine of the 69 families represented in our observed sample had 20 or more members present and therefore we investigate the taxonomy of these families in more detail and find excellent agreement with literature. Out of these 9 well-sampled families, our data show that the Themis, Koronis, Hygiea, Massalia, and Eunomia families display a high degree of taxonomic homogeneity and that the Vesta, Flora, Nysa-Polana, and Eos families show a significant level of mixture in taxonomies. Using the taxonomic purity and the degree of dispersion in observed colors for each of the 9 well-sampled collisional families we also speculate which of these families potentially originated from a differentiated parent body and/or is a family with a possible undetermined nested family. Additionally, we obtained sufficient photometric data for 433 of our targets to extract reliable rotation periods and observe no obvious correlation between rotation properties and family membership.

Imaging and resolved spectroscopy reveal varying environmental conditions in our dynamic solar system. Many key advances have focused on how these conditions change over time. Observatory-level commitments to conduct annual observations of solar system bodies would establish a long-term legacy chronicling the evolution of dynamic planetary atmospheres, surfaces, and rings. Science investigations will use these temporal datasets to address potential biosignatures, circulation and evolution of atmospheres from the edge of the habitable zone to the ice giants, orbital dynamics and planetary seismology with ring systems, exchange between components in the planetary system, and the migration and processing of volatiles on icy bodies, including Ocean Worlds. The common factor among these diverse investigations is the need for a very long campaign duration, and temporal sampling at an annual cadence.

In the absence of dense photometry for a large population of Near Earth Objects (NEOs), the best method of obtaining a shape distribution comes from sparse photometry and partial lightcurves. We have used 867 partial lightcurves obtained by Spitzer to determine a shape distribution for sub-kilometre NEOs. From this data we find a best fit average elongation $\frac{b}{a}=0.72 \pm 0.08$. We compare this result with a shape distribution obtained from 1869 NEOs in the same size range observed by Pan-STARRS 1 and find the Spitzer-obtained elongation to be in excellent agreement with this PS1 value of $\frac{b}{a}=0.70 \pm 0.10$. These values are also in agreement with literature values for $1<D<10$ km objects in the main asteroid belt, however, there is a size discrepancy between the two datasets. Using a smaller sample of NEOs in the size range $1<D<5$ km from PS1 data, we obtain an average axis ratio $b/a = 0.70 \pm 0.12$. This is more elongated than the shape distribution for main belt objects in the same size regime, although the current uncertainties are sizeable and this should be verified using a larger data set. As future large surveys come online it will be possible to observe smaller main belt asteroids to allow for better comparisons of different sub-kilometre populations.

A foundational goal of the Large Synoptic Survey Telescope (LSST) is to map the Solar System small body populations that provide key windows into understanding of its formation and evolution. This is especially true of the populations of the Outer Solar System -- objects at the orbit of Neptune $r > 30$AU and beyond. In this whitepaper, we propose a minimal change to the LSST cadence that can greatly enhance LSST's ability to discover faint distant Solar System objects across the entire wide-fast-deep (WFD) survey area. Specifically, we propose that the WFD cadence be constrained so as to deliver least one sequence of $\gtrsim 10$ visits per year taken in a $\sim 10$ day period in any combination of $g, r$, and $i$ bands. Combined with advanced shift-and-stack algorithms (Whidden et al. 2019) this modification would enable a nearly complete census of the outer Solar System to $\sim 25.5$ magnitude, yielding $4-8$x more KBO discoveries than with single-epoch baseline, and enabling rapid identification and follow-up of unusual distant Solar System objects in $\gtrsim 5$x greater volume of space. These increases would enhance the science cases discussed in Schwamb et al. (2018) whitepaper, including probing Neptune's past migration history as well as discovering hypothesized planet(s) beyond the orbit of Neptune (or at least placing significant constraints on their existence).

We propose an ecliptic Deep Drilling Field that will discover some 10,000~small and faint Kuiper Belt Objects (KBOs) --- primitive rocky/icy bodies that orbit at the outside of our Solar System and uniquely record the processes of planetary system formation and evolution. The primary goals are to measure the KBO size and shape distributions down to 25~km, a size that probes both the early and ongoing evolution of this population. These goals can be met with around 10~hours total of on-sky time (five separate fields that are observed for 2.1~hours each). Additional science will result from downstream observations that provide colors and orbit refinement, for a total time request of 40~hours over the ten year LSST main survey.

1I/`Oumuamua is the first confirmed interstellar body in our Solar System. Here we report on observations of `Oumuamua made with the Spitzer Space Telescope on 2017 November 21--22 (UT). We integrated for 30.2~hours at 4.5 micron (IRAC channel 2). We did not detect the object and place an upper limit on the flux of 0.3 uJy (3sigma). This implies an effective spherical diameter less than [98, 140, 440] meters and albedo greater than [0.2, 0.1, 0.01] under the assumption of low, middle, or high thermal beaming parameter eta, respectively. With an aspect ratio for `Oumuamua of 6:1, these results correspond to dimensions of [240:40, 341:57, 1080:180] meters, respectively. We place upper limits on the amount of dust, CO, and CO2 coming from this object that are lower than previous results; we are unable to constrain the production of other gas species. Both our size and outgassing limits are important because `Oumuamua's trajectory shows non-gravitational accelerations that are sensitive to size and mass and presumably caused by gas emission. We suggest that `Oumuamua may have experienced low-level post-perihelion volatile emission that produced a fresh, bright, icy mantle. This model is consistent with the expected eta value and implied high albedo value for this solution, but, given our strict limits on CO and CO2, requires another gas species --- probably H2O --- to explain the observed non-gravitational acceleration. Our results extend the mystery of `Oumuamua's origin and evolution.

Using the first 18 months of the Pan-STARRS 1 survey we have identified 33 candidate high-amplitude objects for follow-up observations and carried out observations of 22 asteroids. 4 of the observed objects were found to have observed amplitude $A_{obs}\geq 1.0$ mag. We find that these high amplitude objects are most simply explained by single rubble pile objects with some density-dependent internal strength, allowing them to resist mass shedding even at their highly elongated shapes. 3 further objects although below the cut-off for 'high-amplitude' had a combination of elongation and rotation period which also may require internal cohesive strength, depending on the density of the body. We find that none of the 'high-amplitude asteroids' identified here require any unusual cohesive strengths to resist rotational fission. 3 asteroids were sufficiently observed to allow for shape and spin pole models to be determined through light curve inversion. 45864 was determined to have retrograde rotation with spin pole axes $\lambda=218\pm 10^{\circ}, \beta=-82\pm 5^{\circ}$ and asteroid 206167 was found to have best fit spin pole axes $\lambda= 57 \pm 5^{\circ}$, $\beta=-67 \pm 5^{\circ}$. An additional object not initially measured with $A_{obs}>1.0$ mag, 49257, was determined to have a shape model which does suggest a high-amplitude object. Its spin pole axes were best fit for values $\lambda=112\pm 6^{\circ}, \beta=6\pm 5^{\circ}$. In the course of this project to date no large super-fast rotators ($P_{rot} < 2.2$ h) have been identified.

Over 4.5 years, the Mission Accessible Near-Earth Object Survey (MANOS) assembled 228 Near-Earth Object (NEO) lightcurves. We report rotational lightcurves for 82 NEOs, constraints on amplitudes and periods for 21 NEOs, lightcurves with no detected variability within the image signal to noise and length of our observing block for 30 NEOs, and 10 tumblers. We uncovered 2 ultra-rapid rotators with periods below 20s; 2016MA with a potential rotational periodicity of 18.4s, and 2017QG$_{18}$ rotating in 11.9s, and estimate the fraction of fast/ultra-rapid rotators undetected in our project plus the percentage of NEOs with a moderate/long periodicity undetectable during our typical observing blocks. We summarize the findings of a simple model of synthetic NEOs to infer the object morphologies distribution using the measured distribution of lightcurve amplitudes. This model suggests a uniform distribution of axis ratio can reproduce the observed sample. This suggests that the quantity of spherical NEOs (e.g., Bennu) is almost equivalent to the quantity of highly elongated objects (e.g., Itokawa), a result that can be directly tested thanks to shape models from Doppler delay radar imaging analysis. Finally, we fully characterized 2 NEOs as appropriate targets for a potential robotic/human mission: 2013YS$_{2}$ and 2014FA$_{7}$ due to their moderate spin periods and low $\Delta v$.

Gaia Data Release 2 includes observational data for 14,099 pre-selected asteroids. From the sparsely sampled G band photometry, we derive lower-limit lightcurve amplitudes for 11,665 main belt asteroids in order to provide constraints on the distribution of shapes in the asteroid main belt. Assuming a triaxial shape model for each asteroid, defined through the axial aspect ratios a > b and b=c, we find an average b/a=0.80+-0.04 for the ensemble, which is in agreement with previous results. By combining the Gaia data with asteroid properties from the literature, we investigate possible correlations of the aspect ratio with size, semi-major axis, geometric albedo, and intrinsic color. Based on our model simulations, we find that main belt asteroids greater than 50 km in diameter on average have higher b/a aspect ratios (are rounder) than smaller asteroids. We furthermore find significant differences in the shape distribution of main belt asteroids as a function of the other properties that do not affect the average aspect ratios. We conclude that a more detailed investigation of shape distribution correlations requires a larger data sample than is provided in Gaia Data Release 2.

We present lightcurves and derive periods and amplitudes for a subset of 38 near earth objects (NEOs) observed at 4.5 microns with the IRAC camera on the the Spitzer Space Telescope, many of them having no previously reported rotation periods. This subset was chosen from about 1800 IRAC NEO observations as having obvious periodicity and significant amplitude. For objects where the period observed did not sample the full rotational period, we derived lower limits to these parameters based on sinusoidal fits. Lightcurve durations ranged from 42 to 544 minutes, with derived periods from 16 to 400 minutes. We discuss the effects of lightcurve variations on the thermal modeling used to derive diameters and albedos from Spitzer photometry. We find that both diameters and albedos derived from the lightcurve maxima and minima agree with our previously published results, even for extreme objects, showing the conservative nature of the thermal model uncertainties. We also evaluate the NEO rotation rates, sizes, and their cohesive strengths.

We present a community-led assessment of the solar system investigations achievable with NASA's next-generation space telescope, the Wide Field InfraRed Survey Telescope (WFIRST). WFIRST will provide imaging, spectroscopic, and coronagraphic capabilities from 0.43-2.0 $\mu$m and will be a potential contemporary and eventual successor to JWST. Surveys of irregular satellites and minor bodies are where WFIRST will excel with its 0.28 deg$^2$ field of view Wide Field Instrument (WFI). Potential ground-breaking discoveries from WFIRST could include detection of the first minor bodies orbiting in the Inner Oort Cloud, identification of additional Earth Trojan asteroids, and the discovery and characterization of asteroid binary systems similar to Ida/Dactyl. Additional investigations into asteroids, giant planet satellites, Trojan asteroids, Centaurs, Kuiper Belt Objects, and comets are presented. Previous use of astrophysics assets for solar system science and synergies between WFIRST, LSST, JWST, and the proposed NEOCam mission are discussed. We also present the case for implementation of moving target tracking, a feature that will benefit from the heritage of JWST and enable a broader range of solar system observations.

Two planetary mass objects in the far outer Solar System --- collectively referred to here as Planet X --- have recently been hypothesized to explain the orbital distribution of distant Kuiper Belt Objects. Neither planet is thought to be exceptionally faint, but the sky locations of these putative planets are poorly constrained. Therefore, a wide area survey is needed to detect these possible planets. The Large Synoptic Survey Telescope (LSST) will carry out an unbiased, large area (around 18,000 deg$^2$), deep (limiting magnitude of individual frames of 24.5) survey (the "wide-fast-deep" survey) of the southern sky beginning in 2022, and is therefore an important tool to search for these hypothesized planets. Here we explore the effectiveness of LSST as a search platform for these possible planets. Assuming the current baseline cadence (which includes the wide-fast-deep survey plus additional coverage) we estimate that LSST will confidently detect or rule out the existence of Planet X in 61\% of the entire sky. At orbital distances up to $\sim$75 au, Planet X could simply be found in the normal nightly moving object processing; at larger distances, it will require custom data processing. We also discuss the implications of a non-detection of Planet X in LSST data.

The large and tidally-locked classical moons of Uranus display longitudinal and planetocentric trends in their surface compositions. Spectrally red material has been detected primarily on the leading hemispheres of the outer moons, Titania and Oberon. Furthermore, detected H2O ice bands are stronger on the leading hemispheres of the classical satellites, and the leading/trailing asymmetry in H2O ice band strengths decreases with distance from Uranus. We hypothesize that the observed distribution of red material and trends in H2O ice band strengths results from infalling dust from Uranian irregular satellites. These dust particles migrate inward on slowly decaying orbits, eventually reaching the classical satellite zone, where they collide primarily with the outer moons. The latitudinal distribution of dust swept up by these moons should be fairly even across their southern and northern hemispheres. However, red material has only been detected over the southern hemispheres of these moons (subsolar latitude 81 S). Consequently, to test whether irregular satellite dust impacts drive the observed enhancement in reddening, we have gathered new ground-based data of the now observable northern hemispheres of these moons (sub-observer latitudes, 17 to 35 N). Our results and analyses indicate that longitudinal and planetocentric trends in reddening and H2O ice band strengths are broadly consistent across both southern and northern latitudes of these moons, thereby supporting our hypothesis. Utilizing a suite of numerical best fit models, we investigate the composition of the reddening agent detected on these moons, finding that both complex organics and amorphous pyroxene match the spectral slopes of our data. We also present spectra that span 2.9 to 4.1 microns, a previously unexplored wavelength range in terms of spectroscopy for the Uranian moons.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

We present VRI spectrophotometry of 1003 Main-Belt Asteroids (MBAs) observed with the Sutherland, South Africa, node of the Korea Microlensing Telescope Network (KMTNet). All of the observed MBAs were serendipitously captured in KMTNet's large 2deg $\times$ 2deg field of view during a separate targeted near-Earth Asteroid study (Erasmus et al. 2017). Our broadband spectrophotometry is reliable enough to distinguish among four asteroid taxonomies and we confidently categorize 836 of the 1003 observed targets as either a S-, C-, X-, or D-type asteroid by means of a Machine Learning (ML) algorithm approach. Our data show that the ratio between S-type MBAs and (C+X+D)-type MBAs, with H magnitudes between 12 and 18 (12 km $\gtrsim$ diameter $\gtrsim$ 0.75 km), is almost exactly 1:1. Additionally, we report 0.5- to 3-hour (median: 1.3-hour) light-curve data for each MBA and we resolve the complete rotation periods and amplitudes for 59 targets. Two out of the 59 targets have rotation periods potentially below the theoretical zero cohesion boundary limit of 2.2 hours. We report lower limits for the rotation periods and amplitudes for the remaining targets. Using the resolved and unresolved light curves we determine the shape distribution for this population using a Monte Carlo simulation. Our model suggests a population with an average elongation $b/a = 0.74\pm0.07$ and also shows that this is independent of asteroid size and taxonomy.

1I/'Oumuamua was discovered by the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS 1) on 19 October 2017. Unlike all previously discovered minor planets this object was determined to have eccentricity $e > 1.0$, suggesting an interstellar origin. Since this discovery and within the limited window of opportunity, several photometric and spectroscopic studies of the object have been made. Using the measured light curve amplitudes and rotation periods we find that, under the assumption of a triaxial ellipsoid, a density range $1500 < \rho < 2800$ kg m$^{-3}$ matches the observations and no significant cohesive strength is required. We also determine that an aspect ratio of $6\pm 1:1$ is most likely after accounting for phase-angle effects and considering the potential effect of surface properties. This elongation is still remarkable but less than some other estimates.

The Large Synoptic Survey Telescope (LSST) is uniquely equipped to search for Solar System bodies due to its unprecedented combination of depth and wide field coverage. Over a ten-year period starting in 2022, LSST will generate the largest catalog of Solar System objects to date. The main goal of the LSST Solar System Science Collaboration (SSSC) is to facilitate the efforts of the planetary community to study the planets and small body populations residing within our Solar System using LSST data. To prepare for future survey cadence decisions and ensure that interesting and novel Solar System science is achievable with LSST, the SSSC has identified and prioritized key Solar System research areas for investigation with LSST in this roadmap. The ranked science priorities highlighted in this living document will inform LSST survey cadence decisions and aid in identifying software tools and pipelines needed to be developed by the planetary community as added value products and resources before the planned start of LSST science operations.

The majority of known asteroid diameters are derived from thermal-infrared observations. Diameters are derived using asteroid thermal models that approximate their surface temperature distributions and compare the measured thermal-infrared flux with model-dependent predictions. The most commonly used thermal model is the Near-Earth Asteroid Thermal Model (NEATM), which is usually perceived as superior to other models like the Fast-Rotating Model (FRM). We investigate the applicability of the NEATM and the FRM to thermal-infrared observations of Near-Earth Objects using synthetic asteroids with properties based on the real Near-Earth Asteroid (NEA) population. We find the NEATM to provide more accurate diameters and albedos than the FRM in most cases, with a few exceptions. The modeling results are barely affected by the physical properties of the objects, but we find a large impact of the solar phase angle on the modeling results. We conclude that the NEATM provides statistically more robust diameter estimates for NEAs observed at solar phase angles less than ~65{\deg}, while the FRM provides more robust diameter estimates for solar phase angles greater than ~65{\deg}. We estimate that <5% of all NEA diameters and albedos derived up to date are affected by systematic effects that are of the same order of magnitude as the typical thermal model uncertainties. We provide statistical correction functions for diameters and albedos derived using the NEATM and FRM as a function of solar phase angle.

The recently discovered minor body 1I/2017 U1 (`Oumuamua) is the first known object in our Solar System that is not bound by the Sun's gravity. Its hyperbolic orbit (eccentricity greater than unity) strongly suggests that it originated outside our Solar System; its red color is consistent with substantial space weathering experienced over a long interstellar journey. We carry out an simple calculation of the probability of detecting such an object. We find that the observed detection rate of 1I-like objects can be satisfied if the average mass of ejected material from nearby stars during the process of planetary formation is ~20 Earth masses, similar to the expected value for our Solar System. The current detection rate of such interstellar interlopers is estimated to be 0.2/year, and the expected number of detections over the past few years is almost exactly one. When the Large Synoptic Survey Telescope begins its wide, fast, deep all-sky survey the detection rate will increase to 1/year. Those expected detections will provide further constraints on nearby planetary system formation through a better estimate of the number and properties of interstellar objects.

Although the final observations of the Spitzer Warm Mission are currently scheduled for March 2019, it can continue operations through the end of the decade with no loss of photometric precision. As we will show, there is a strong science case for extending the current Warm Mission to December 2020. Spitzer has already made major impacts in the fields of exoplanets (including microlensing events), characterizing near Earth objects, enhancing our knowledge of nearby stars and brown dwarfs, understanding the properties and structure of our Milky Way galaxy, and deep wide-field extragalactic surveys to study galaxy birth and evolution. By extending Spitzer through 2020, it can continue to make ground-breaking discoveries in those fields, and provide crucial support to the NASA flagship missions JWST and WFIRST, as well as the upcoming TESS mission, and it will complement ground-based observations by LSST and the new large telescopes of the next decade. This scientific program addresses NASA's Science Mission Directive's objectives in astrophysics, which include discovering how the universe works, exploring how it began and evolved, and searching for life on planets around other stars.

Context. A lot of photometric data is produced by surveys such as Pan-STARRS, LONEOS, WISE or Catalina. These data are a rich source of information about the physical properties of asteroids. There are several possible approaches for utilizing these data. Lightcurve inversion is a typical method that works with individual asteroids. Our approach in this paper is statistical when we focused on large groups of asteroids like dynamical families and taxonomic classes, and the data were not sufficient for individual models. Aims. Our aim was to study the distributions of shape elongation $b/a$ and the spin axis latitude $\beta$ for various subpopulations of asteroids and to compare our results, based on Pan-STARRS1 survey, with statistics previously done using different photometric data (Lowell database, WISE data). Methods. We use the LEADER algorithm to compare the $b/a$ and $\beta$ distributions for different subpopulations of asteroids. The algorithm creates a cumulative distributive function (CDF) of observed brightness variations, and computes the $b/a$ and $\beta$ distributions using analytical basis functions that yield the observed CDF. A variant of LEADER is used to solve the joint distributions for synthetic populations to test the validity of the method. Results. When comparing distributions of shape elongation for groups of asteroids with different diameters $D$, we found that there are no differences for $D < 25$ km. We also constructed distributions for asteroids with different rotation periods and revealed that the fastest rotators with $P = 0 - 4$ h are more spheroidal than the population with $P = 4 - 8$ h.

We present here VRI spectrophotometry of 39 near-Earth asteroids (NEAs) observed with the Sutherland, South Africa, node of the Korea Microlensing Telescope Network (KMTNet). Of the 39 NEAs, 19 were targeted, but because of KMTNet's large 2 deg by 2 deg field of view, 20 serendipitous NEAs were also captured in the observing fields. Targeted observations were performed within 44 days (median: 16 days, min: 4 days) of each NEA's discovery date. Our broadband spectrophotometry is reliable enough to distinguish among four asteroid taxonomies and we were able to confidently categorize 31 of the 39 observed targets as either a S-, C-, X- or D-type asteroid by means of a Machine Learning (ML) algorithm approach. Our data suggest that the ratio between "stony" S-type NEAs and "not-stony" (C+X+D)-type NEAs, with H magnitudes between 15 and 25, is roughly 1:1. Additionally, we report ~1-hour light curve data for each NEA and of the 39 targets we were able to resolve the complete rotation period and amplitude for six targets and report lower limits for the remaining targets.

The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible.

We analyzed data from the first year of a survey for Near Earth Objects (NEOs) that we are carrying out with the Dark Energy Camera (DECam) on the 4-meter Blanco telescope at the Cerro Tololo Inter-American Observatory. We implanted synthetic NEOs into the data stream to derive our nightly detection efficiency as a function of magnitude and rate of motion. Using these measured efficiencies and the Solar System absolute magnitudes derived by the Minor Planet Center for the 1377 measurements of 235 unique NEOs detected, we directly derive, for the first time from a single observational data set, the NEO size distribution from 1 km down to 10 meters. We find that there are 10^6.6 NEOs larger than 10 meters. This result implies a factor of ten fewer small NEOs than some previous results, though our derived size distribution is in good agreement with several other estimates.

The JAXA Hayabusa-2 mission was approved in 2010 and launched on December 3, 2014. The spacecraft will arrive at the near-Earth asteroid 162173 Ryugu in 2018 where it will perform a survey, land and obtain surface material, then depart in Dec 2019 and return to Earth in Dec 2020. We observed Ryugu with the Herschel Space Observatory in Apr 2012 at far-IR thermal wavelengths, supported by several ground-based observations to obtain optical lightcurves. We reanalysed previously published Subaru-COMICS and AKARI-IRC observations and merged them with a Spitzer-IRS data set. In addition, we used a large set of Spitzer-IRAC observations obtained in the period Jan to May, 2013. The data set includes two complete rotational lightcurves and a series of ten "point-and-shoot" observations. The almost spherical shape of the target together with the insufficient lightcurve quality forced us to combine radiometric and lightcurve inversion techniques in different ways to find the object's key physical and thermal parameters. We find that the solution which best matches our data sets leads to this C class asteroid having a retrograde rotation with a spin-axis orientation of (lambda = 310-340 deg; beta = -40+/-15 deg) in ecliptic coordinates, an effective diameter (of an equal-volume sphere) of 850 to 880 m, a geometric albedo of 0.044 to 0.050 and a thermal inertia in the range 150 to 300 Jm-2s-0.5K-1. Based on estimated thermal conductivities of the top-layer surface in the range 0.1 to 0.6 WK-1m-1, we calculated that the grain sizes are approximately equal to between 1 and 10 mm. The finely constrained values for this asteroid serve as a `design reference model', which is currently used for various planning, operational and modelling purposes by the Hayabusa2 team.

The Large Synoptic Survey Telescope (LSST) will be a discovery machine for the astronomy and physics communities, revealing astrophysical phenomena from the Solar System to the outer reaches of the observable Universe. While many discoveries will be made using LSST data alone, taking full scientific advantage of LSST will require ground-based optical-infrared (OIR) supporting capabilities, e.g., observing time on telescopes, instrumentation, computing resources, and other infrastructure. This community-based study identifies, from a science-driven perspective, capabilities that are needed to maximize LSST science. Expanding on the initial steps taken in the 2015 OIR System Report, the study takes a detailed, quantitative look at the capabilities needed to accomplish six representative LSST-enabled science programs that connect closely with scientific priorities from the 2010 decadal surveys. The study prioritizes the resources needed to accomplish the science programs and highlights ways that existing, planned, and future resources could be positioned to accomplish the science goals.

Near Earth Objects (NEOs) are small Solar System bodies whose orbits bring them close to the Earth's orbit. We are carrying out a Warm Spitzer Cycle 11 Exploration Science program entitled NEOSurvey --- a fast and efficient flux-limited survey of 597 known NEOs in which we derive diameter and albedo for each target. The vast majority of our targets are too faint to be observed by NEOWISE, though a small sample has been or will be observed by both observatories, which allows for a cross-check of our mutual results. Our primary goal is to create a large and uniform catalog of NEO properties. We present here the first results from this new program: fluxes and derived diameters and albedos for 80 NEOs, together with a description of the overall program and approach, including several updates to our thermal model. The largest source of error in our diameter and albedo solutions, which derive from our single band thermal emission measurements, is uncertainty in eta, the beaming parameter used in our thermal modeling; for albedos, improvements in Solar System absolute magnitudes would also help significantly. All data and derived diameters and albedos from this entire program are being posted on a publicly accessible webpage at nearearthobjects.nau.edu .

The Mission Accessible Near-Earth Objects Survey (MANOS) aims to physically characterize sub-km Near-Earth Objects (NEOs). We report first photometric results from the survey which began in August, 2013. Photometric observations were performed using 1 m to 4 m class telescopes around the world. We present rotational periods and lightcurve amplitudes for 86 sub-km NEOs, though in some cases, only lower limits are provided. Our main goal is to obtain lightcurves for small NEOs (typically, sub-km objects) and estimate their rotational periods, lightcurve amplitudes, and shapes. These properties are used for statistical study to constrain overall properties of the NEO population. A weak correlation seems to indicate that smaller objects are more spherical than the larger ones. We also report 7 NEOs that are fully characterized (lightcurve and visible spectra) as the most suitable candidates for a future human or robotic mission. Viable mission targets are objects fully characterized, with a Delta_v(NHATS) <12 km s^-1, and a rotational period P>1h. Assuming a similar rate of object characterization as reported in this paper, approximately 1,230 NEOs need to be characterized in order to find 100 viable mission targets.

Using the Wide Field Camera for the United Kingdom Infrared Telescope, we measure the near-infrared colors of near-Earth objects (NEOs) in order to put constraints on their taxonomic classifications. The rapid-response character of our observations allows us to observe NEOs when they are close to the Earth and bright. Here we present near-infrared color measurements of 86 NEOs, most of which were observed within a few days of their discovery, allowing us to characterize NEOs with diameters of only a few meters. Using machine-learning methods, we compare our measurements to existing asteroid spectral data and provide probabilistic taxonomic classifications for our targets. Our observations allow us to distinguish between S-complex, C/X-complex, D-type, and V-type asteroids. Our results suggest that the fraction of S-complex asteroids in the whole NEO population is lower than the fraction of ordinary chondrites in the meteorite fall statistics. Future data obtained with UKIRT will be used to investigate the significance of this discrepancy.

Data from the New Horizons mission to Pluto show no craters on Sputnik Planum down to the detection limit (2 km for low resolution data, 625 m for high resolution data). The number of small Kuiper Belt Objects that should be impacting Pluto is known to some degree from various astronomical surveys. We combine these geological and telescopic observations to make an order of magnitude estimate that the surface age of Sputnik Planum must be less than 10 million years. This maximum surface age is surprisingly young and implies that this area of Pluto must be undergoing active resurfacing, presumably through some cryo-geophysical process. We discuss three possible resurfacing mechanisms and the implications of each one for Pluto's physical properties.

The rubble pile spin barrier is an upper limit on the rotation rate of asteroids larger than ~200-300 m. Among thousands of asteroids with diameters larger than ~300 m, only a handful of asteroids are known to rotate faster than 2.0 h, all are in the sub-km range (<=0.6 km). Here we present photometric measurements suggesting that (60716) 2000 GD65, an S-complex, inner-main belt asteroid with a relatively large diameter of 2.3 +0.6-0.7 km, completes one rotation in 1.9529+-0.0002 h. Its unique diameter and rotation period allow us to examine scenarios about asteroid internal structure and evolution: a rubble pile bound only by gravity; a rubble-pile with strong cohesion; a monolithic structure; an asteroid experiencing mass shedding; an asteroid experiencing YORP spin-up/down; and an asteroid with a unique octahedron shape results with a four-peak lightcurve and a 3.9 h period. We find that the most likely scenario includes a lunar-like cohesion that can prevent (60716) 2000 GD65 from disrupting without requiring a monolithic structure or a unique shape. Due to the uniqueness of (60716) 2000 GD65, we suggest that most asteroids typically have smaller cohesion than that of lunar regolith.

We perform a search for dormant comets, asteroidal objects of cometary origin, in the near-Earth asteroid (NEA) population based on dynamical and physical considerations. Our study is based on albedos derived within the ExploreNEOs program and is extended by adding data from NEOWISE and the Akari asteroid catalog. We use a statistical approach to identify asteroids on orbits that resemble those of short-period near-Earth comets using the Tisserand parameter with respect to Jupiter, the aphelion distance, and the minimum orbital intersection distance with respect to Jupiter. From the sample of NEAs on comet-like orbits, we select those with a geometric albedo $p_V \leq 0.064$ as dormant comet candidates, and find that only $\sim$50% of NEAs on comet-like orbits also have comet-like albedos. We identify a total of 23 NEAs from our sample that are likely to be dormant short-period near-Earth comets and, based on a de-biasing procedure applied to the cryogenic NEOWISE survey, estimate both magnitude-limited and size-limited fractions of the NEA population that are dormant short-period comets. We find that 0.3-3.3% of the NEA population with $H \leq 21$, and $9^{+2}_{-5}$% of the population with diameters $d \geq 1$ km, are dormant short-period near-Earth comets.

The surfaces of the large Uranian satellites are characterized by a mixture of H2O ice and a dark, potentially carbon-rich, constituent, along with CO2 ice. At the mean heliocentric distance of the Uranian system, native CO2 ice should be removed on timescales shorter than the age of the Solar System. Consequently, the detected CO2 ice might be actively produced. Analogous to irradiation of icy moons in the Jupiter and Saturn systems, we hypothesize that charged particles caught in Uranus' magnetic field bombard the surfaces of the Uranian satellites, driving a radiolytic CO2 production cycle. To test this hypothesis, we investigated the distribution of CO2 ice by analyzing near-infrared (NIR) spectra of these moons, gathered using the SpeX spectrograph at NASA's Infrared Telescope Facility (IRTF) (2000 - 2013). Additionally, we made spectrophotometric measurements using images gathered by the Infrared Array Camera (IRAC) onboard the Spitzer Space Telescope (2003 - 2005). We find that the detected CO2 ice is primarily on the trailing hemispheres of the satellites closest to Uranus, consistent with other observations of these moons. Our band parameter analysis indicates that the detected CO2 ice is pure and segregated from other constituents. Our spectrophotometric analysis indicates that IRAC is not sensitive to the CO2 ice detected by SpeX, potentially because CO2 is retained beneath a thin surface layer dominated by H2O ice that is opaque to photons over IRAC wavelengths. Thus, our combined SpeX and IRAC analyses suggest that the near-surfaces (i.e., top few 100 microns) of the Uranian satellites are compositionally stratified. We briefly compare the spectral characteristics of the CO2 ice detected on the Uranian moons to icy satellites elsewhere, and we also consider the most likely drivers of the observed distribution of CO2 ice.

Large-area surveys operating at mid-infrared wavelengths have proven to be a valuable means of discovering and characterizing minor planets. Through the use of radiometric models, it is possible to derive physical properties such as diameters, albedos, and thermal inertia for large numbers of objects. Modern detector array technology has resulted in a significant improvement in spatial resolution and sensitivity compared with previous generations of space-based infrared telescopes, giving rise to a commensurate increase in the number of objects that have been observed at these wavelengths. Space-based infrared surveys of asteroids therefore offer an effective means of rapidly gathering information about small body populations' orbital and physical properties. The AKARI, WISE/NEOWISE, Spitzer, and Herschel missions have significantly increased the number of minor planets with well-determined diameters and albedos.

We report on observations of near-Earth asteroid 2011 MD with the Spitzer Space Telescope. We have spent 19.9 h of observing time with channel 2 (4.5 {\mu}m) of the Infrared Array Camera and detected the target within the 2{\sigma} positional uncertainty ellipse. Using an asteroid thermophysical model and a model of nongravitational forces acting upon the object we constrain the physical properties of 2011 MD, based on the measured flux density and available astrometry data. We estimate 2011 MD to be 6 (+4/-2) m in diameter with a geometric albedo of 0.3 (+0.4/-0.2) (uncertainties are 1{\sigma}). We find the asteroid's most probable bulk density to be 1.1 (+0.7/-0.5) g cm^{-3}, which implies a total mass of (50-350) t and a macroporosity of >=65%, assuming a material bulk density typical of non-primitive meteorite materials. A high degree of macroporosity suggests 2011 MD to be a rubble-pile asteroid, the rotation of which is more likely to be retrograde than prograde.

The Deep Ecliptic Survey (DES) discovered hundreds of Kuiper Belt objects from 1998-2005. Follow-up observations yielded 304 objects with good dynamical classifications (Classical, Scattered, Centaur, or 16 mean-motion resonances with Neptune). The DES search fields are well documented, enabling us to calculate the probability of detecting objects with particular orbital parameters and absolute magnitudes at a randomized point in each orbit. Grouping objects together by dynamical class leads, we estimate the orbital element distributions (a, e, i) for the largest three classes (Classical, 3:2, and Scattered) using maximum likelihood. Using H-magnitude as a proxy for the object size, we fit a power law to the number of objects for 8 classes with at least 5 detected members (246 objects). The best Classical slope is alpha=1.02+/-0.01 (observed from 5<=H<=7.2). Six dynamical classes (Scattered plus 5 resonances) are consistent in slope with the Classicals, though the absolute number of objects is scaled. The exception to the power law relation are the Centaurs (non-resonant with perihelia closer than Neptune, and thus detectable at smaller sizes), with alpha=0.42+/-0.02 (7.5<H<11). This is consistent with a knee in the H-distribution around H=7.2 as reported elsewhere (Bernstein et al. 2004, Fraser et al. 2014). Based on the Classical-derived magnitude distribution, the total number of objects (H<=7) in each class are: Classical (2100+/-300 objects), Scattered (2800+/-400), 3:2 (570+/-80), 2:1 (400+/-50), 5:2 (270+/-40), 7:4 (69+/-9), 5:3 (60+/-8). The independent estimate for the number of Centaurs in the same H range is 13+/-5. If instead all objects are divided by inclination into "Hot" and "Cold" populations, following Fraser et al. (2014), we find that alphaHot=0.90+/-0.02, while alphaCold=1.32+/-0.02, in good agreement with that work.

We report on Spitzer Space Telescope IRAC observations of near-Earth object (NEO) 2009 BD that were carried out in support of the NASA Asteroid Robotic Retrieval Mission (ARRM) concept. We did not detect 2009 BD in 25 hrs of integration at 4.5 micron. Based on an upper-limit flux density determination from our data, we present a probabilistic derivation of the physical properties of this object. The analysis is based on the combination of a thermophysical model with an orbital model accounting for the non-gravitational forces acting upon the body. We find two physically possible solutions. The first solution shows 2009 BD as a 2.9+/-0.3 m diameter rocky body (rho = 2.9+/-0.5 g cm-3) with an extremely high albedo of 0.85(+0.20/-0.10) that is covered with regolith-like material, causing it to exhibit a low thermal inertia (Gamma = 30(+20/-10) SI units). The second solution suggests 2009 BD to be a 4+/-1 m diameter asteroid with pV = 0.45(+0.35/-0.15) that consists of a collection of individual bare rock slabs (Gamma = 2000+/-1000 SI units, rho = 1.7(+0.7/-0.4) g cm-3). We are unable to rule out either solution based on physical reasoning. 2009 BD is the smallest asteroid for which physical properties have been constrained, in this case using an indirect method and based on a detection limit, providing unique information on the physical properties of objects in the size range smaller than 10 m.

Centaurs are the transitional population between trans-Neptunian objects (TNOs) and Jupiter-family comets. For this reason it is possible to access the smaller ones, which is more difficult to do with the TNO population. The goal of this work is to characterize a set of 16 Centaurs in terms of their size, albedo, and thermal properties. We study the correlations, for a more extended sample obtained from the literature, of diameter, albedo, orbital parameters, and spectral slopes. We performed three-band photometric observations using Herschel-PACS and used a consistent method for the data reduction and aperture photometry of this sample to obtain monochromatic flux densities at 70, 100, and 160 $\mu$m. Additionally, we used Spitzer-MIPS flux densities at 24 and 70 $\mu$m when available. We also included in our Centaur sample scattered disk objects (SDOs), a dynamical family of TNOS, using results previously published by our team, and some Centaurs observed only with the Spitzer/MIPS instrument. We have determined new radiometric sizes and albedos of 16 Centaurs. The first conclusion is that the albedos of Centaur objects are not correlated with their orbital parameters. Similarly, there is no correlation between diameter and orbital parameters. Most of the objects in our sample are dark (pv < 7%) and most of them are small (D < 120km). However, there is no correlation between albedo and diameter, in particular for the group of the small objects as albedo values are homogeneously distributed between 4 to 16%. The correlation with the color of the objects showed that red objects are all small (mean diameter 65 km), while the gray ones span a wide range of sizes (mean diameter 120 km). Moreover, the gray objects tend to be darker, with a mean albedo of 5.6%, compared with a mean of 8.5% (ranging from 5 to 15%) for the red objects.

The near-Earth object (NEO) population, which mainly consists of fragments from collisions between asteroids in the main asteroid belt, is thought to include contributions from short-period comets as well. One of the most promising NEO candidates for a cometary origin is near-Earth asteroid (3552) Don Quixote, which has never been reported to show activity. Here we present the discovery of cometary activity in Don Quixote based on thermal-infrared observations made with the Spitzer Space Telescope in its 3.6 and 4.5 {\mu}m bands. Our observations clearly show the presence of a coma and a tail in the 4.5 {\mu}m but not in the 3.6 {\mu}m band, which is consistent with molecular band emission from CO2. Thermal modeling of the combined photometric data on Don Quixote reveals a diameter of 18.4 (-0.4/+0.3) km and an albedo of 0.03 (-0.01/+0.02), which confirms Don Quixote to be the third-largest known NEO. We derive an upper limit on the dust production rate of 1.9 kg s^-1 and derive a CO2 gas production rate of (1.1+-0.1)10^26 molecules s^-1. Spitzer IRS spectroscopic observations indicate the presence of fine-grained silicates, perhaps pyroxene rich, on the surface of Don Quixote. Our discovery suggests that CO2 can be present in near-Earth space over a long time. The presence of CO2 might also explain that Don Quixote's cometary nature remained hidden for nearly three decades.

By analyzing brightness variation with ecliptic longitude and using the Lowell Observatory photometric database, we estimate spin-axis longitudes for more than 350 000 asteroids. Hitherto, spin-axis longitude estimates have been made for fewer than 200 asteroids. We investigate longitude distributions in different dynamical groups and asteroid families. We show that asteroid spin-axis longitudes are not isotropically distributed as previously considered. We find that the spin-axis longitude distribution for main-belt asteroids is clearly non-random, with an excess of longitudes from the interval 30{\deg}-110{\deg} and a paucity between 120{\deg}-180{\deg}. The explanation of the non-isotropic distribution is unknown at this point. Further studies have to be conducted to determine if the shape of the distribution can be explained by observational bias, selection effects, a real physical process or other mechanism.

Near-infrared spectroscopy of Near-Earth Objects (NEOs) connects diagnostic spectral features to specific surface mineralogies. The combination of spectroscopy with albedos and diameters derived from thermal infrared observations can increase the scientific return beyond that of the individual datasets. To that end, we have completed a spectroscopic observing campaign to complement the ExploreNEOs Warm Spitzer program that obtained albedos and diameters of nearly 600 NEOs (Trilling et al. 2010). Here we present the results of observations using the low-resolution prism mode (~0.7-2.5 microns) of the SpeX instrument on the NASA Infrared Telescope Facility (IRTF). We also include near-infrared observations of ExploreNEOs targets from the MIT-UH-IRTF Joint Campaign for Spectral Reconnaissance. Our dataset includes near-infrared spectra of 187 ExploreNEOs targets (125 observations of 92 objects from our survey and 213 observations of 154 objects from the MIT survey). We identify a taxonomic class for each spectrum and use band parameter analysis to investigate the mineralogies for the S-, Q-, and V-complex objects. Our analysis suggests that for spectra that contain near-infrared data but lack the visible wavelength region, the Bus-DeMeo system misidentifies some S-types as Q-types. We find no correlation between spectral band parameters and ExploreNEOs albedos and diameters. We find slightly negative Band Area Ratio (BAR) correlations with phase angle for Eros and Ivar, but a positive BAR correlation with phase angle for Ganymed. We find evidence for spectral phase reddening for Eros, Ganymed, and Ivar. We identify the likely ordinary chondrite type analog for a subset of our sample. Our resulting proportions of H, L, and LL ordinary chondrites differ from those calculated for meteorite falls and in previous studies of ordinary chondrite-like NEOs.

We show, by comparing observations with theoretical models, that the observed Kuiper Belt size distribution is well matched by coagulation models, which start from an initial planetesimal population with radii of about 1km, and subsequent collisional evolution. We find that the observed size distribution for R > 30km has not been modified by collisional evolution over the age of the solar system, and that the size distribution below R ~ 30km has been modified by collisions and that its slope is well matched by collisional evolution models that use published strength laws. We investigate in detail the resulting size distribution of bodies ranging from 0.01km to 30km and find that its slope changes several times as a function of radius before approaching the expected value for an equilibrium collisional cascade of material strength dominated bodies for R < 0.1km. Compared to a single power law size distribution that would span the whole range from 0.01km to 30km, we find a strong deficit of bodies around R ~10km and a strong excess of bodies around 2km. This deficit and excess are caused by an excess mass in small planetesimals in the km size range that was left over from the runaway growth phase and that leaves a signature in the size distribution that is not erased after 4.5 Gyrs of collisional evolution. Observations of the small KBO size distribution can therefore test if large KBOs grew as a result of runaway growth and constrain the initial planetesimal sizes. We find that results from recent KBO occultation surveys and the observed KBO size distribution can be best matched by an initial planetesimal population that contained about equal mass per logarithmic mass bin in bodies ranging from 0.4km to 4km in radius. In addition, we find that we cannot match the observed KBO size distribution if most of the planetesimal mass was contained in bodies that were 10km in radius or larger.

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.

We use Spitzer Space Telescope 24 {\mu}m data to search for debris disks among 122 AFGKM stars from the \sim 670 Myr clusters Hyades, Coma Ber, and Praesepe, utilizing a number of advances in data reduction and determining the intrinsic colors of main sequence stars. For our sample, the 1{\sigma} dispersion about the main sequence V-K, K-[24] locus is approximately 3.1%. We identify seven debris disks at 10% or more (\geq 3{\sigma} confidence level) above the expected K-[24] for purely photospheric emission. The incidence of excesses of 10% or greater in our sample at this age is 5.7 +3.1/-1.7%. Combining with results from the literature, the rate is 7.8 +4.2/-2.1% for early- type (B9 - F4) stars and 2.7 +3.3/-1.7% for solar-like (F5 - K9) stars. Our primary sample has strict criteria for inclusion to allow comparison with other work; when we relax these criteria, three additional debris disks are detected. They are all around stars of solar-like type and hence reinforce our conclusion that disks around such stars are still relatively common at 670 Myr and are similar to the rate around early-type stars. The apparently small difference in decay rates between early-type and solar-like stars is inconsistent with the first order theoretical predictions that the later type stellar disks would decay an order of magnitude more quickly than the earlier type ones.

We explore the correlation between an asteroid's taxonomy and photometric phase curve using the H, G12 photometric phase function, with the shape of the phase function described by the single parameter G12. We explore the usability of G12 in taxonomic classification for individual objects, asteroid families, and dynamical groups. We conclude that the mean values of G12 for the considered taxonomic complexes are statistically different, and also discuss the overall shape of the G12 distribution for each taxonomic complex. Based on the values of G12 for about half a million asteroids, we compute the probabilities of C, S, and X complex membership for each asteroid. For an individual asteroid, these probabilities are rather evenly distributed over all of the complexes, thus preventing meaningful classification. We then present and discuss the G12 distributions for asteroid families, and predict the taxonomic complex preponderance for asteroid families given the distribution of G12 in each family. For certain asteroid families, the probabilistic prediction of taxonomic complex preponderance can clearly be made. The Nysa-Polana family shows two distinct regions in the proper element space with different G12 values dominating in each region. We conclude that the G12-based probabilistic distribution of taxonomic complexes through the main belt agrees with the general view of C complex asteroid proportion increasing towards the outer belt. We conclude that the G12 photometric parameter cannot be used in determining taxonomic complex for individual asteroids, but it can be utilized in the statistical treatment of asteroid families and different regions of the main asteroid belt.

We present the faintest mid ecliptic latitude survey in the second part of HST archival search for outer Solar System bodies. We report the discovery of 28 new trans-Neptunian objects and 1 small centaur (R ~ 2km) in the band 5{\circ} - 20{\circ} off the ecliptic. The inclination distribution of these excited ob jects is consistent with the distribution derived from brighter ecliptic surveys. We suggest that the size and inclination distribution should be estimated consistently using suitable surveys with calibrated search algorithms and reliable orbital information.

Space missions to NEOs are being planned at all major space agencies, and recently a manned mission to an NEO was announced as a NASA goal. Efforts to find and select suitable targets (plus backup targets) are severely hampered by our lack of knowledge of the physical properties of dynamically favorable NEOs. In particular, current mission scenarios tend to favor primitive low-albedo objects. For the vast majority of NEOs the albedo is unknown. Here we report new constraints on the size and albedo of 65 NEOs with rendezvous deltaV < 7 km/s. Our results are based on thermal-IR flux data obtained in the framework of our ongoing (2009--2011) ExploreNEOs survey (Trilling et al. 2010) using NASA's "Warm Spitzer" space telescope. As of 2010 July 14, we have results for 293 objects in hand (including the 65 low-deltaV NEOs presented here); before the end of 2011 we expect to have measured the size and albedo of ~ 700 NEOs (including probably ~ 160 low-deltaV NEOs). While there are reasons to believe that primitive volatile-rich materials are universally low in albedo, the converse need not be true: the orbital evolution of some dark objects likely has caused them to lose their volatiles by coming too close to the Sun. For all our targets, we give the closest perihelion distance they are likely to have reached (using orbital integrations from Marchi et al. 2009) and corresponding upper limits on the past surface temperature. Low-deltaV objects for which both albedo and thermal history may suggest a primitive composition include (162998) 2001 SK162, (68372) 2001 PM9, and (100085) 1992 UY4.

We combine nulling interferometry at 10 {\mu}m using the MMT and Keck Telescopes with spectroscopy, imaging, and photometry from 3 to 100 {\mu}m using Spitzer to study the debris disk around {\beta} Leo over a broad range of spatial scales, corresponding to radii of 0.1 to ~100 AU. We have also measured the close binary star o Leo with both Keck and MMT interferometers to verify our procedures with these instruments. The {\beta} Leo debris system has a complex structure: 1.) relatively little material within 1 AU; 2.) an inner component with a color temperature of ~600 K, fitted by a dusty ring from about 2 to 3 AU; and 3.) a second component with a color temperature of ~120 K fitted by a broad dusty emission zone extending from about ~5 AU to ~55 AU. Unlike many other A-type stars with debris disks, {\beta} Leo lacks a dominant outer belt near 100 AU.

Bernstein et al. (2004) found that the population of faint (R>26) trans-Neptunian objects (TNOs) known at that time was dominated by "Classical" objects, which have low inclinations (i<5 degrees) and distances 40--45 AU. Since those observations, the number of faint TNOs whose orbits are sufficiently well known to be classified as "Classical" or "Excited" has grown from seven to 39. We analyze the dynamical classifications of faint TNOs known today and find that this population is dominated by Excited objects. We discuss some implications of this result.

A debris ring around the star HD 207129 (G0V; d = 16.0 pc) has been imaged in scattered visible light with the ACS coronagraph on the Hubble Space Telescope and in thermal emission using MIPS on the Spitzer Space Telescope at 70 microns (resolved) and 160 microns (unresolved). Spitzer IRS (7-35 microns) and MIPS (55-90 microns) spectrographs measured disk emission at >28 microns. In the HST image the disk appears as a ~30 AU wide ring with a mean radius of ~163 AU and is inclined by 60 degrees from pole-on. At 70 microns it appears partially resolved and is elongated in the same direction and with nearly the same size as seen with HST in scattered light. At 0.6 microns the ring shows no significant brightness asymmetry, implying little or no forward scattering by its constituent dust. With a mean surface brightness of V=23.7 mag per square arcsec, it is the faintest disk imaged to date in scattered light.

We introduce a novel search technique that can identify trans-neptunian objects in three to five exposures of a pointing within a single Hubble Space Telescope orbit. The process is fast enough to allow the discovery of candidates soon after the data are available. This allows sufficient time to schedule follow up observations with HST within a month. We report the discovery of 14 slow-moving objects found within 5\circ of the ecliptic in archival data taken with the Wide Field Channel of the Advanced Camera for Surveys. The luminosity function of these objects is consistent with previous ground-based and space-based results. We show evidence that the size distribution of both high and low inclination populations is similar for objects smaller than 100 km, as expected from collisional evolution models, while their size distribution differ for brighter objects. We suggest the two populations formed in different parts of the protoplanetary disk and after being dynamically mixed have collisionally evolved together. Among the objects discovered there is an equal mass binary with an angular separation ~ 0."53.

We have begun the ExploreNEOs project in which we observe some 700 Near Earth Objects (NEOs) at 3.6 and 4.5 microns with the Spitzer Space Telescope in its Warm Spitzer mode. From these measurements and catalog optical photometry we derive albedos and diameters of the observed targets. The overall goal of our ExploreNEOs program is to study the history of near-Earth space by deriving the physical properties of a large number of NEOs. In this paper we describe both the scientific and technical construction of our ExploreNEOs program. We present our observational, photometric, and thermal modeling techniques. We present results from the first 101 targets observed in this program. We find that the distribution of albedos in this first sample is quite broad, probably indicating a wide range of compositions within the NEO population. Many objects smaller than one kilometer have high albedos (>0.35), but few objects larger than one kilometer have high albedos. This result is consistent with the idea that these larger objects are collisionally older, and therefore possess surfaces that are more space weathered and therefore darker, or are not subject to other surface rejuvenating events as frequently as smaller NEOs.

Thermal emission from Kuiper Belt object (136108) Haumea was measured with Herschel-PACS at 100 and 160 micrometers for almost a full rotation period. Observations clearly indicate a 100-micrometer thermal lightcurve with an amplitude of a factor of ~ 2, which is positively correlated with the optical lightcurve. This confirms that both are primarily due to shape effects. A 160-micrometer lightcurve is marginally detected. Radiometric fits of the mean Herschel- and Spitzer- fluxes indicate an equivalent diameter D ~ 1300 km and a geometric albedo p_v ~ 0.70-0.75. These values agree with inferences from the optical lightcurve, supporting the hydrostatic equilibrium hypothesis. The large amplitude of the 100-micrometer lightcurve suggests that the object has a high projected a/b axis ratio (~ 1.3) and a low thermal inertia as well as possible variable infrared beaming. This may point to fine regolith on the surface, with a lunar-type photometric behavior. The quality of the thermal data is not sufficient to clearly detect the effects of a surface dark spot.

The goal of the Herschel Open Time Key programme "TNOs are Cool!" is to derive the physical and thermal properties for a large sample of Centaurs and trans-Neptunian objects (TNOs), including resonant, classical, detached and scattered disk objects. We present results for seven targets either observed in PACS point-source, or in mini scan-map mode. Spitzer-MIPS observations were included for three objects. The sizes of these targets range from 100 km to almost 1000 km, five have low geometric albedos below 10%, (145480) 2005 TB190 has a higher albedo above 15%. Classical thermal models driven by an intermediate beaming factor of $\eta$=1.2 or $\eta$-values adjusted to the observed colour temperature fit the multi-band observations well in most cases. More sophisticated thermophysical models give very similar diameter and albedo values for thermal inertias in the range 0-25 Jm-2s-0.5K-1, consistent with very low heat conductivities at temperatures far away from the Sun. The early experience with observing and model strategies will allow us to derive physical and thermal properties for our complete Herschel TNO sample of 140 targets as a benchmark for understanding the solar system debris disk, and extra-solar ones as well.

Using data from the Deep Ecliptic Survey (DES), we investigate the inclination distributions of objects in the Kuiper Belt. We present a derivation for observational bias removal and use this procedure to generate unbiased inclination distributions for Kuiper Belt objects (KBOs) of different DES dynamical classes, with respect to the Kuiper Belt Plane. Consistent with previous results, we find that the inclination distribution for all DES KBOs is well fit by the sum of two Gaussians, or a Gaussian plus a generalized Lorentzian, multiplied by sin i. Approximately 80% of KBOs are in the high-inclination grouping. We find that Classical object inclinations are well fit by sin i multiplied by the sum of two Gaussians, with roughly even distribution between Gaussians of widths 2.0 -0.5/+0.6 degrees and 8.1 -2.1/+2.6 degrees. Objects in different resonances exhibit different inclination distributions. The inclinations of Scattered objects are best matched by sin i multiplied by a single Gaussian that is centered at 19.1 -3.6/+3.9 degrees with a width of 6.9 -2.7/+4.1 degrees. Centaur inclinations peak just below 20 degrees, with one exceptionally high-inclination object near 80 degrees. The currently observed inclination distribution of the Centaurs is not dissimilar to that of the Scattered Extended KBOs and Jupiter-family comets, but is significantly different from the Classical and Resonant KBOs. While the sample sizes of some dynamical classes are still small, these results should begin to serve as a critical diagnostic for models of Solar System evolution.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

We present 24 micron photometry of the intermediate-age open cluster Praesepe. We assemble a catalog of 193 probable cluster members that are detected in optical databases, the Two Micron All Sky Survey (2MASS), and at 24 micron, within an area of ~ 2.47 square degrees. Mid-IR excesses indicating debris disks are found for one early-type and for three solar-type stars. Corrections for sampling statistics yield a 24 micron excess fraction (debris disk fraction) of 6.5 +- 4.1% for luminous and 1.9 +- 1.2% for solar-type stars. The incidence of excesses is in agreement with the decay trend of debris disks as a function of age observed for other cluster and field stars. The values also agree with those for older stars, indicating that debris generation in the zones that emit at 24 micron falls to the older 1-10 Gyr field star sample value by roughly 750 Myr. We discuss our results in the context of previous observations of excess fractions for early- and solar-type stars. We show that solar-type stars lose their debris disk 24 micron excesses on a shorter timescale than early-type stars. Simplistic Monte Carlo models suggest that, during the first Gyr of their evolution, up to 15-30% of solar-type stars might undergo an orbital realignment of giant planets such as the one thought to have led to the Late Heavy Bombardment, if the length of the bombardment episode is similar to the one thought to have happened in our Solar System. In the Appendix, we determine the cluster's parameters via boostrap Monte Carlo isochrone fitting, yielding an age of 757 Myr (+- 36 Myr at 1 sigma confidence) and a distance of 179 pc (+- 2 pc at 1 sigma confidence), not allowing for systematic errors.

Near Earth Objects (NEOs) are fragments of remnant primitive bodies that date from the era of Solar System formation. At present, the physical properties and origins of NEOs are poorly understood. We have measured thermal emission from three NEOs -- (6037) 1988 EG, 1993 GD, and 2005 GL -- with Spitzer's IRAC instrument at 3.6, 4.5, 5.8, and 8.0 microns (the last object was detected only at 5.8 and 8.0 microns). The diameters of these three objects are 400 m, 180 m, and 160 m, respectively, with uncertainties of around 20% (including both observational and systematic errors). For all three the geometric albedos are around 0.30, in agreement with previous results that most NEOs are S-class asteroids. For the two objects detected at 3.6 and 4.5 microns, diameters and albedos based only on those data agree with the values based on modeling the data in all four bands. This agreement, and the high sensitivity of IRAC, show the promise of the Spitzer Warm Mission for determining the physical parameters for a large number of NEOs.

Spitzer images resolve the debris disk around \gamma Ophiuchi at both 24 and 70 um. The resolved images suggest a disk radius of ~520 AU at 70 um and >=260 AU at 24 um. The images, along with a consistent fit to the spectral energy distribution of the disk from 20 to 350 um, show that the primary disk structure is inclined by ~50 degree from the plane of the sky at a position angle of 55+/-2 degree. Among a group of twelve debris disks that have similar host star spectral types, ages and infrared fractional luminosities, the observed sizes in the infrared and color temperatures indicate that evolution of the debris disks is influenced by multiple parameters in addition to the proto-planetary disk initial mass.

The overarching goal of planetary astronomy is to deduce how the present collection of objects found in our Solar System were formed from the original material present in the proto-solar nebula. As over two hundred exo-planetary systems are now known, and multitudes more are expected, the Solar System represents the closest and best system which we can study, and the only one in which we can clearly resolve individual bodies other than planets. In this White Paper we demonstrate how to use Spitzer Space Telescope InfraRed Array Camera Channels 1 and 2 (3.6 and 4.5 um) imaging photometry with large dedicated surveys to advance our knowledge of Solar System formation and evolution. There are a number of vital, key projects to be pursued using dedicated large programs that have not been pursued during the five years of Spitzer cold operations. We present a number of the largest and most important projects here; more will certainly be proposed once the warm era has begun, including important observations of newly discovered objects.

Detecting heat from minor planets in the outer solar system is challenging, yet it is the most efficient means for constraining the albedos and sizes of Kuiper Belt Objects (KBOs) and their progeny, the Centaur objects. These physical parameters are critical, e.g., for interpreting spectroscopic data, deriving densities from the masses of binary systems, and predicting occultation tracks. Here we summarize Spitzer Space Telescope observations of 47 KBOs and Centaurs at wavelengths near 24 and 70 microns. We interpret the measurements using a variation of the Standard Thermal Model (STM) to derive the physical properties (albedo and diameter) of the targets. We also summarize the results of other efforts to measure the albedos and sizes of KBOs and Centaurs. The three or four largest KBOs appear to constitute a distinct class in terms of their albedos. From our Spitzer results, we find that the geometric albedo of KBOs and Centaurs is correlated with perihelion distance (darker objects having smaller perihelia), and that the albedos of KBOs (but not Centaurs) are correlated with size (larger KBOs having higher albedos). We also find hints that albedo may be correlated with with visible color (for Centaurs). Interestingly, if the color correlation is real, redder Centaurs appear to have higher albedos. Finally, we briefly discuss the prospects for future thermal observations of these primitive outer solar system objects.

We have observed nearly 200 FGK stars at 24 and 70 microns with the Spitzer Space Telescope. We identify excess infrared emission, including a number of cases where the observed flux is more than 10 times brighter than the predicted photospheric flux, and interpret these signatures as evidence of debris disks in those systems. We combine this sample of FGK stars with similar published results to produce a sample of more than 350 main sequence AFGKM stars. The incidence of debris disks is 4.2% (+2.0/-1.1) at 24 microns for a sample of 213 Sun-like (FG) stars and 16.4% (+2.8/-2.9) at 70 microns for 225 Sun-like (FG) stars. We find that the excess rates for A, F, G, and K stars are statistically indistinguishable, but with a suggestion of decreasing excess rate toward the later spectral types; this may be an age effect. The lack of strong trend among FGK stars of comparable ages is surprising, given the factor of 50 change in stellar luminosity across this spectral range. We also find that the incidence of debris disks declines very slowly beyond ages of 1 billion years.

We observed the Mars Trojan asteroids (5261) Eureka and (101429) 1998 VF31 and the candidate Mars Trojan 2001 FR127 at 11.2 and 18.1 microns using Michelle on the Gemini North telescope. We derive diameters of 1.28, 0.78, and <0.52 km, respectively, with corresponding geometric visible albedos of 0.39, 0.32, and >0.14. The albedos for Eureka and 1998 VF31 are consistent with the taxonomic classes and compositions (S(I)/angritic and S(VII)/achrondritic, respectively) and implied histories presented in a companion paper by Rivkin et al. Eureka's surface likely has a relatively high thermal inertia, implying a thin regolith that is consistent with predictions and the small size that we derive.

Mars is the only terrestrial planet known to have Tro jan (co-orbiting) asteroids, with a confirmed population of at least 4 objects. The origin of these objects is not known; while several have orbits that are stable on solar-system timescales, work by Rivkin et al. (2003) showed they have compositions that suggest separate origins from one another. We have obtained infrared (0.8-2.5 micron) spectroscopy of the two largest L5 Mars Tro jans, and confirm and extend the results of Rivkin et al. (2003). We suggest that the differentiated angrite meteorites are good spectral analogs for 5261 Eureka, the largest Mars Trojan. Meteorite analogs for 101429 1998 VF31 are more varied and include primitive achondrites and mesosiderites.

We report the mid- and far-infrared properties of nearby M dwarfs. Spitzer/MIPS measurements were obtained for a sample of 62 stars at 24 um, with subsamples of 41 and 20 stars observed at 70 um and 160 um respectively. We compare the results with current models of M star photospheres and look for indications of circumstellar dust in the form of significant deviations of K-[24 um] colors and 70 um / 24 um flux ratios from the average M star values. At 24 um, all 62 of the targets were detected; 70 um detections were achieved for 20 targets in the subsample observed; and no detections were seen in the 160 um subsample. No clear far-infrared excesses were detected in our sample. The average far infrared excess relative to the photospheric emission of the M stars is at least four times smaller than the similar average for a sample of solar-type stars. However, this limit allows the average fractional infrared luminosity in the M-star sample to be similar to that for more massive stars. We have also set low limits for the maximum mass of dust possible around our stars.

Hubble Space Telescope observations of Uranus- and Neptune-crossing object (65489) Ceto/Phorcys (provisionally designated 2003 FX128) reveal it to be a close binary system. The mutual orbit has a period of 9.554 +/- 0.011 days and a semimajor axis of 1840 +/- 48 km. These values enable computation of a system mass of (5.41 +/- 0.42) 10^18 kg. Spitzer Space Telescope observations of thermal emission at 24 and 70 microns are combined with visible photometry to constrain the system's effective radius (109 +10/-11 km) and geometric albedo (0.084 +0.021/-0.014). We estimate the average bulk density to be 1.37 +0.66/-0.32 g cm^-3, consistent with ice plus rocky and/or carbonaceous materials. This density contrasts with lower densities recently measured with the same technique for three other comparably-sized outer Solar System binaries (617) Patroclus, (26308) 1998 SM165, and (47171) 1999 TC36, and is closer to the density of the saturnian irregular satellite Phoebe. The mutual orbit of Ceto and Phorcys is nearly circular, with an eccentricity <= 0.015. This observation is consistent with calculations suggesting that the system should tidally evolve on a timescale shorter than the age of the solar system.

The dust disks observed around mature stars are evidence that plantesimals are present in these systems on spatial scales that are similar to that of the asteroids and the KBOs in the Solar System. These dust disks (a.k.a. ``debris disks'') present a wide range of sizes, morphologies and properties. It is inferred that their dust mass declines with time as the dust-producing planetesimals get depleted, and that this decline can be punctuated by large spikes that are produced as a result of individual collisional events. The lack of solid state features indicate that, generally, the dust in these disks have sizes larger than approximately 10 microns, but exceptionally, strong silicate features in some disks suggest the presence of large quantities of small grains, thought to be the result of recent collisions. Spatially resolved observations of debris disks show a diversity of structural features, such as inner cavities, warps, offsets, brightness asymmetries, spirals, rings and clumps. There is growing evidence that, in some cases, these structures are the result of the dynamical perturbations of a massive planet. Our Solar System also harbors a debris disk and some of its properties resemble those of extra-solar debris disks. From the cratering record, we can infer that its dust mass has decayed with time, and that there was at least one major ``spike'' in the past during the Late Heavy Bombardment. This offers a unique opportunity to use extra-solar debris disks to shed some light in how the Solar System might have looked in the past. Similarly, our knowledge of the Solar System is influencing our understanding of the types of processes which might be at play in the extra-solar debris disks.

We observed 69 A3-F8 main sequence binary star systems using the Multiband Imaging Photometer for Spitzer onboard the Spitzer Space Telescope. We find emission significantly in excess of predicted photospheric flux levels for 9(+4/-3)% and 40(+7/-6)% of these systems at 24 and 70 microns, respectively. Twenty two systems total have excess emission, including four systems that show excess emission at both wavelengths. A very large fraction (nearly 60%) of observed binary systems with small (<3 AU) separations have excess thermal mission. We interpret the observed infrared excesses as thermal emission from dust produced by collisions in planetesimal belts. The incidence of debris disks around main sequence A3-F8 binaries is marginally higher than that for single old AFGK stars. Whatever combination of nature (birth conditions of binary systems) and nurture (interactions between the two stars) drives the evolution of debris disks in binary systems, it is clear that planetesimal formation is not inhibited to any great degree. We model these dust disks through fitting the spectral energy distributions and derive typical dust temperatures in the range 100--200 K and typical fractional luminosities around 10^-5, with both parameters similar to other Spitzer-discovered debris disks. Our calculated dust temperatures suggest that about half the excesses we observe are derived from circumbinary planetesimal belts and around one third of the excesses clearly suggest circumstellar material. Three systems with excesses have dust in dynamically unstable regions, and we discuss possible scenarios for the origin of this short-lived dust.

We present 24 micron Spitzer/MIPS photometric observations of the ~50 Myr open cluster IC 2391. Thirty-four cluster members ranging in spectral type from B3-M5 were observed in the central square degree of the cluster. Excesses indicative of debris disks were discovered around 1 A star, 6 FGK stars, and possibly 1 M dwarf. For the cluster members observed to their photospheric limit, we find a debris disk frequency of 10 (-3,+17)% for B-A stars and 31 (-9,+13)% for FGK stars using a 15% relative excess threshold. Relative to a model of decaying excess frequency, the frequency of debris disks around A-type stars appears marginally low for the cluster's age while that of FGK stars appears consistent. Scenarios that may qualitatively explain this result are examined. We conclude that planetesimal activity in the terrestrial region of FGK stars is common in the first ~50 Myr and decays on timescales of ~100 Myr. Despite luminosity differences, debris disk evolution does not appear to depend strongly on stellar mass.

We report 24 and/or 70 um measurements of ~160 A-type main-sequence stars using the Multiband Imaging Photometer for Spitzer (MIPS). Their ages range from 5 to 850 Myr based on estimates from the literature (cluster or moving group associations) or from the H-R diagram and isochrones. The thermal infrared excess is identified by comparing the deviation (~3% and ~15% at the 1-\sigma level at 24 and 70 um, respectively) between the measurements and the synthetic Kurucz photospheric predictions. Stars showing excess infrared emission due to strong emission lines or extended nebulosity seen at 24 um are excluded from our sample; therefore, the remaining infrared excesses are likely to arise from circumstellar debris disks. At the 3-sigma confidence level, the excess rate at 24 and 70 um is 32% and >=33% (with an uncertainty of 5%), considerably higher than has been found for old solar analogs and M dwarfs. Our measurements place constraints on the fractional dust luminosities and temperatures in the disks. We find that older stars tend to have lower fractional dust luminosity than younger ones. While the fractional luminosity from the excess infrared emission follows a general 1/t relationship, the values at a given stellar age vary by at least two orders of magnitude. We also find that (1) older stars possess a narrow range of temperature distribution peaking at colder temperatures, and (2) the disk emission at 70 um persists longer than that at 24 um. Both results suggest that the debris-disk clearing process is more effective in the inner regions.

With the MIPS camera on Spitzer, we have searched for far-infrared emission from dust in the planetary system orbiting pulsar PSR 1257+12. With accuracies of 0.05 mJy at 24 um and 1.5 mJy at 70 um, photometric measurements find no evidence for emission at these wavelengths. These observations place new upper limits on the luminosity of dust with temperatures between 20 and 1000 K. They are particularly sensitive to dust temperatures of 100-200 K, for which they limit the dust luminosity to below $3 \times 10^{-5}$ of the pulsar's spin-down luminosity, three orders of magnitude better than previous limits. Despite these improved constraints on dust emission, an asteroid belt similar to the Solar System's cannot be ruled out.

We present observations of six late-type members of the young cluster IC 348 detected at 24 microns with the Multiband Imaging Photometer for Spitzer(MIPS). At least four of the objects are probably substellar. Combining these data with ground-based optical and near-infrared photometry and complementary observations with the Infrared Array Camera (IRAC), we have modeled the spectral energy distributions using detailed models of irradiated accretion disks. We are able to fit the observations with models using a range of maximum grain sizes from ISM-type dust to grains as large as 1 millimeter. Two objects show a lack of excess emission at wavelengths shortward of 5.8-8 microns but significant excess at longer wavelengths, indicative of large optically thin or evacuated inner holes. Our models indicate a inner hole of radius ~ 0.5-0.9 AU for the brown dwarf L316; this is the first brown dwarf with evidence for an AU-scale inner disk hole. We examine several possible mechanisms for the inner disk clearing in this case, including photoevaporation and planet formation.

We measured the system-integrated thermal emission of the binary Kuiper Belt Object 1999 TC36 at wavelengths near 24 and 70 microns using the Spitzer space telescope. We fit these data and the visual magnitude using both the Standard Thermal Model and thermophysical models. We find that the effective diameter of the binary is 405 km, with a range of 350 -- 470 km, and the effective visible geometric albedo for the system is 0.079 with a range of 0.055 -- 0.11. The binary orbit, magnitude contrast between the components, and system mass have been determined from HST data (Margot et al., 2004; 2005a; 2005b). Our effective diameter, combined with that system mass, indicate an average density for the objects of 0.5 g/cm3, with a range 0.3 -- 0.8 g/cm3. This density is low compared to that of materials expected to be abundant in solid bodies in the trans-Neptunian region, requiring 50 -- 75% of the interior of 1999 TC36 be taken up by void space. This conclusion is not greatly affected if 1999 TC36 is ``differentiated'' (in the sense of having either a rocky or just a non-porous core). If the primary is itself a binary, the average density of that (hypothetical) triple system would be in the range 0.4 -- 1.1 g/cm3, with a porosity in the range 15 -- 70%.

We report the results of a spectroscopic search for debris disks surrounding 41 nearby solar type stars, including 8 planet-bearing stars, using the {\it Spitzer Space Telescope}. With accurate relative photometry using the Infrared Spectrometer (IRS) between 7-34 $\micron$ we are able to look for excesses as small as $\sim$2% of photospheric levels with particular sensitivity to weak spectral features. For stars with no excess, the $3\sigma$ upper limit in a band at 30-34 $\mu$m corresponds to $\sim$ 75 times the brightness of our zodiacal dust cloud. Comparable limits at 8.5-13 $\mu$m correspond to $\sim$ 1,400 times the brightness of our zodiacal dust cloud. These limits correspond to material located within the $<$1 to $\sim$5 AU region that, in our solar system, originates from debris associated with the asteroid belt. We find excess emission longward of $\sim$25 $\mu$m from five stars of which four also show excess emission at 70 $\mu$m. This emitting dust must be located around 5-10 AU. One star has 70 micron emission but no IRS excess. In this case, the emitting region must begin outside 10 AU; this star has a known radial velocity planet. Only two stars of the five show emission shortward of 25 $\micron$ where spectral features reveal the presence of a population of small, hot dust grains emitting in the 7-20 $\mu$m band. The data presented here strengthen the results of previous studies to show that excesses at 25 $\micron$ and shorter are rare: only 1 star out of 40 stars older than 1 Gyr or $\sim 2.5$% shows an excess. Asteroid belts 10-30 times more massive than our own appear are rare among mature, solar-type stars.

We report the results of a spectroscopic search for debris disks surrounding 41 nearby solar type stars, including 8 planet-bearing stars, using the {\it Spitzer Space Telescope}. With accurate relative photometry using the Infrared Spectrometer (IRS) between 7-34 $\micron$ we are able to look for excesses as small as $\sim$2% of photospheric levels with particular sensitivity to weak spectral features. For stars with no excess, the $3\sigma$ upper limit in a band at 30-34 $\mu$m corresponds to $\sim$ 75 times the brightness of our zodiacal dust cloud. Comparable limits at 8.5-13 $\mu$m correspond to $\sim$ 1,400 times the brightness of our zodiacal dust cloud. These limits correspond to material located within the $<$1 to $\sim$5 AU region that, in our solar system, originates from debris associated with the asteroid belt. We find excess emission longward of $\sim$25 $\mu$m from five stars of which four also show excess emission at 70 $\mu$m. This emitting dust must be located around 5-10 AU. One star has 70 micron emission but no IRS excess. In this case, the emitting region must begin outside 10 AU; this star has a known radial velocity planet. Only two stars of the five show emission shortward of 25 $\micron$ where spectral features reveal the presence of a population of small, hot dust grains emitting in the 7-20 $\mu$m band. The data presented here strengthen the results of previous studies to show that excesses at 25 $\micron$ and shorter are rare: only 1 star out of 40 stars older than 1 Gyr or $\sim 2.5$% shows an excess. Asteroid belts 10-30 times more massive than our own appear are rare among mature, solar-type stars.

We report high precision photometry of three small and one larger Kuiper Belt Objects (KBOs) obtained with the Advanced Camera for Surveys onboard the Hubble Space Telescope (ACS/HST). The three small bodies are the smallest KBOs for which lightcurve measurements are available. 2003 BF91 has a diameter of 20 kilometers (assuming 10% albedo) and a 1.09 magnitude, 9.1-hour lightcurve that is feasibly explained by the rotation of an elongated, coherent body that is supported by material strength and best imagined as an icy outer Solar System analog to asteroid (243) Ida. Two other small KBOs, 2003 BG91 and 2003 BH91 (diameters 31 and 18 km, with albedo 10%), exhibit an unremarkable lightcurve and no detectable photometric variation, respectively. For the larger KBO 2000 FV53 (116 km diameter, assuming 10% albedo) we strongly detect a non-sinusoidal periodic (7.5 hours) brightness variation with a very small amplitude (0.07 mag). This KBO may be nearly spherical, a result that might not be unusual in the Kuiper Belt but would be remarkable among outer Solar System satellites of similar size. We carry out a study of possible physical states and bulk densities under the assumptions of both fluid equilibrium and finite, non-zero internal friction. The densities for the these KBOs are likely to be in the range 1--2 g/cm3, and a plausible solution for 2000 FV53 is a rubble pile of this density that is held slightly out of the minimum-energy shape by internal friction among constituent blocks that are relatively small. Our interpretation of 2000 FV53 as a pulverized but essentially primordial object and 2003 BF91 as a collisional fragment is consistent with models of collisional timescales in the outer Solar System. We compile all published KBO lightcurve data and compare our results to the larger population. [abridged]

We have searched for infrared excesses around a well defined sample of 69 FGK main-sequence field stars. These stars were selected without regard to their age, metallicity, or any previous detection of IR excess; they have a median age of ~4 Gyr. We have detected 70 um excesses around 7 stars at the 3-sigma confidence level. This extra emission is produced by cool material (< 100 K) located beyond 10 AU, well outside the ``habitable zones'' of these systems and consistent with the presence of Kuiper Belt analogs with ~100 times more emitting surface area than in our own planetary system. Only one star, HD 69830, shows excess emission at 24 um, corresponding to dust with temperatures > 300 K located inside of 1 AU. While debris disks with Ld/L* > 10^-3 are rare around old FGK stars, we find that the disk frequency increases from 2+-2% for Ld/L* > 10^-4 to 12+-5% for Ld/L* > 10^-5. This trend in the disk luminosity distribution is consistent with the estimated dust in our solar system being within an order of magnitude, greater or less, than the typical level around similar nearby stars.

Spitzer photometry and spectroscopy of the star HD69830 reveal an excess of emission relative to the stellar photosphere between 8 and 35 $\mu$m dominated by strong features attributable to crystalline silicates with an emitting surface area more than 1,000 times that of our zodiacal cloud. The spectrum closely resembles that of the comet C/1995 O1 (Hale-Bopp). Since no excess is detected at 70 $\mu$m, the emitting material must be quite warm, be confined within a few AU of the star, and originate in grains with low long-wavelength emissivity, i.e. grains much smaller than $70 \mu$m$/2\pi\sim 10 $ $\mu$m. The strong mineralogical features are evidence for even smaller, possibly sub-micron sized grains. This small grain size is in direct contrast to the 10-100 $\mu$m grains that dominate the relatively featureless spectra of our zodiacal dust cloud and most other main sequence stars with excesses. The upper limit at 70 $\mu$m also implies that any Kuiper Belt analog must either be very cold or be less massive than $\sim$ 5 times our own Kuiper Belt. With collisional and Poynting-Robertson drag times of about a thousand years for small grains, the emitting material must either: a) be created through continual grinding down of material in a dense asteroid belt; or b) originate in cometary debris arising from either a single ``super-comet'' or a very large number of individual comets arriving from a distant reservoir. In the case of a cometary origin for the emission, the mass requirements for continuous generation by many individual comets are unreasonable and we favor the capture of a single ``super comet'' into a 0.5-1 AU orbit where it can evolve a large number of small grains over a 2 Myr period.

We present high spatial resolution mid- and far-infrared images of the Vega debris disk obtained with the Multiband Imaging Photometer for Spitzer (MIPS). The disk is well resolved and its angular size is much larger than found previously. The radius of the disk is at least 43" (330 AU), 70"(543 AU), and 105" (815 AU) in extent at 24, 70 and 160 um, respectively. The disk images are circular, smooth and without clumpiness at all three wavelengths. The radial surface brightness profiles imply an inner boundary at a radius of 11"+/-2" (86 AU). Assuming an amalgam of amorphous silicate and carbonaceous grains, the disk can be modeled as an axially symmetric and geometrically thin disk, viewed face-on, with the surface particle number density following an r^-1 power law. The disk radiometric properties are consistent with a range of models using grains of sizes ~1 to ~50 um. We find that a ring, containing grains larger than 180 um and at radii of 86-200 AU from the star, can reproduce the observed 850 um flux, while its emission does not violate the observed MIPS profiles. This ring could be associated with a population of larger asteroidal bodies analogous to our own Kuiper Belt. Cascades of collisions starting with encounters amongthese large bodies in the ring produce the small debris that is blown outward by radiation pressure to much larger distances where we detect its thermal emission. The dust production rate is >~10^15 g/s based on the MIPS results. This rate would require a very massive asteroidal reservoir for the dust to be produced in a steady state throughout Vega's life. Instead, we suggest that the disk we imaged is ephemeral and that we are witnessing the aftermath of a large and relatively recent collisional event, and subsequent collisional cascade.

We report our efforts to constrain the form of the low-mass star and brown dwarf mass function via Bayesian inference. Recent surveys of M, L, and T dwarfs in the local solar neighborhood are an essential component of our study. Uncertainties in the age distribution of local field stars make reliable inference complicated. We adopt a wide range of plausible assumptions about the rate of galactic star formation and show that their deviations from a uniform rate produce little effect on the resulting luminosity function for a given mass function. We use a Bayesian statistical formalism to evaluate the probability of commonly used mass functions in light of recent discoveries. We consider three functional forms of the mass function, include a two-segment power law, a single power law with a low-mass cutoff, and a log-normal distribution. Our results show that, at a 60% confidence level, the power-law index, $\alpha$, for the low-mass arm of a two-segment power law has a value between -0.5 and 0.5 for objects with masses between $0.04 M_{\odot}$ and $0.10 M_{\odot}$. The best-fit index is $\alpha = 0.3\pm0.6$ at the 60% confidence level for a single-segment mass function. Current data require this function extend to at least $0.05 M_{\odot}$ with no restrictions placed on a lower mass cutoff. Inferences of the parameter values for a log-normal mass function are virtually unaffected by recent estimates of the local space density of L and T dwarfs. We find no preference among these three forms using this method. We discuss current and future capabilities that may eventually discriminate between mass-function models and refine estimates of their associated parameter values.

As part of a large Spitzer/MIPS GTO program, we have searched for infrared excesses due to debris disks toward 26 FGK field stars known from radial velocity (RV) studies to have one or more planets. While none of these stars show excesses at 24um, we have detected 70um excesses around 6 stars at the 3-sigma confidence level. The excesses are produced by cool material (< 100 K) located beyond 10 AU, well outside the "habitable zones" of these systems and consistent with the presence of Kuiper Belt analogues with ~100 times more emitting surface area than in our own planetary system. These planet-bearing stars are, by selection for RV studies, typically older than 1 Gyr, and the stars identified here with excesses have a median age of 4 Gyr. We find a preliminary correlation of both the frequency and the magnitude of dust emission with the presence of known planets. These are the first stars outside the solar system identified as having both well-confirmed planetary systems and well-confirmed IR excesses.

[Condensed] We search 0.02 deg^2 for trans-Neptunian objects (TNOs) with m<=29.2 (diameter ~15 km) using the ACS on HST. Three new objects are discovered, roughly 25 times fewer than expected from extrapolation of the differential sky density Sigma(m) of brighter objects. The ACS and other recent TNO surveys show departures from a power law size distribution. Division of the TNO sample into ``classical Kuiper belt'' (CKB) and ``Excited'' samples reveals that Sigma(m) differs for the two populations at 96% confidence. A double power law adequately fits all data. Implications include: The total mass of the CKB is ~0.010 M_Earth, only a few times Pluto's mass, and is predominately in the form of ~100 km bodies. The mass of Excited objects is perhaps a few times larger. The Excited class has a shallower bright-end size distribution; the largest objects, including Pluto, comprise tens of percent of the total mass whereas the largest CKBOs are only ~2% of its mass. The predicted mass of the largest Excited body is close to the Pluto mass; the largest CKBO is ~60 times less massive. The deficit of small TNOs occurs for sizes subject to disruption by present-day collisions, suggesting extensive depletion by collisions. Both accretion and erosion appearing to have proceeded to more advanced stages in the Excited class than the CKB. The absence of distant TNOs implies that any distant (60 AU) population must have less than the CKB mass in the form of objects 40 km or larger. The CKB population is sparser than theoretical estimates of the required precursor population for short period comets, but the Excited population could be a viable precursor population.

We numerically investigate the possibility that a close stellar encounter could account for the high inclinations of the Kuiper belt, as originally proposed by Ida, Larwood and Burkert, however we consider encounters with pericenters within those explored by previous works. A star of mass 0.2 Msol entering the solar system on a nearly parabolic, low inclination, retrograde orbit with perihelion of 50 AU can scatter ~30% of the Kuiper belt into moderate inclination (i>10 degrees) and eccentricity orbits, while leaving the rest of the belt at low eccentricities and inclinations. This scenario provides a possible explanation for the dichotomy of the Kuiper belt, accounting for the scattered and classical populations, including objects with high eccentricities and inclinations. High eccentricity objects that were not placed into resonance with Neptune are more likely to be removed from the system by subsequent long timescale evolution. Consequently this scenario suggests that some Plutinos could have originally been objects in circular orbits which were scattered into resonant high eccentricity and high inclination orbits by the stellar encounter. Objects could have been placed into stable resonant regions with Neptune by the encounter. However long lived regions of dynamical stability should differ from those selected by resonance capture. This scenario predicts the presence of objects in resonances with Neptune that are unlikely to have been captured via Neptune's migration.

The Deep Ecliptic Survey is a project whose goal is to survey a large area of the near-ecliptic region to a faint limiting magnitude (R ~ 24) in search of objects in the outer solar system. We are collecting a large homogeneous data sample from the Kitt Peak Mayall 4-m and Cerro Tololo Blanco 4-m telescopes with the Mosaic prime-focus CCD cameras. Our goal is to collect a sample of 500 objects with good orbits to further our understanding of the dynamical structure of the outer solar system. This survey has been in progress since 1998 and is responsible for 272 designated discoveries as of March 2003. We summarize our techniques, highlight recent results, and describe publically available resources.

We predict near-infrared luminosity functions of young (5 Myr to 1 Gyr) star clusters by combining evolutionary models of very low-mass ($1 M_J$ to $0.15 M_{\odot}$) dwarfs with empirical bolometric corrections. We identify several characteristic features in our results. These can be attributed to three causes: (1) deuterium burning in the most massive substellar objects; (2) methane absorption in bodies with $T_{eff}$ less than 1300 K, the temperature of the L/T transition; and (3) the formation of dust clouds and the rainout of dust at roughly the same effective temperature as methane formation. Accurate reconstruction of the substellar mass function from luminosity function observations requires that these phenomena are taken into account. At present, few observational studies extend to sufficient sensitivities to allow detection of these effects. However, the luminosity function of the young open cluster IC 2391 shows a clear peak at $M_I \sim 14$ which we attribute to the result of deuterium burning in substellar objects. The location of this feature is a strong function of age, and we estimate an age of 35 Myr for IC 2391. This is significantly younger than the 53 Myr derived from the location of the lithium depletion boundary but agrees with the main sequence turnoff age. We consider the implications of this result and our multi-band luminosity functions for future observational studies. All predicted luminosity function features are, or will be, accessible to observations using new wide-field IR imagers and the Space Infrared Telescope Facility.

As part of our ongoing Deep Ecliptic Survey (DES) of the Kuiper belt, we report on the occupation of the 1:1 (Trojan), 4:3, 3:2, 7:4, 2:1, and 5:2 Neptunian mean-motion resonances (MMRs). The occupation of the 1:1 and 5:2 MMRs is not easily understood within the standard model of resonance sweeping by a migratory Neptune over an initially dynamically cold belt. Our dynamically hot, 5:2 resonant objects can librate with modest amplitudes of 90 deg within the resonance for at least 1 Gyr. Their trajectories cannot be explained by close encounters with Neptune alone, given the latter's current orbit. The hot orbits of such 5:2 resonant KBOs, unlike hot orbits of previously known resonant KBOs, may imply that these objects were pre-heated to large inclination and large eccentricity prior to resonance capture by a migratory Neptune. Our first discovered Neptunian Trojan, 2001QR322, may not owe its existence to Neptune's migration at all. The trajectory of 2001QR322 is remarkably stable; the object can undergo tadpole-type libration about Neptune's leading Lagrange (L4) point for at least 1 Gyr with a libration amplitude of 24 deg. Trojan capture probably occurred while Neptune accreted the bulk of its mass. For an assumed albedo of 12--4%, our Trojan is 130--230 km in diameter. Model-dependent estimates place the total number of Neptune Trojans resembling 2001QR322 at 20--60. Their existence might rule out violent orbital histories for Neptune.

We present a statistical study of the post-formation migration of giant planets in a range of initial disk conditions. For given initial conditions we model the evolution of giant planet orbits under the influence of disk, stellar, and mass loss torques. We determine the mass and semi-major axis distribution of surviving planets after disk dissipation, for various disk masses, lifetimes, viscosities, and initial planet masses. The majority of planets migrate too fast and are destroyed via mass transfer onto the central star. Most surviving planets have relatively large orbital semi-major axes of several AU or larger. We conclude that the extrasolar planets observed to date, particularly those with small semi-major axes, represent only a small fraction (~25% to 33%) of a larger cohort of giant planets around solar-type stars, and many undetected giant planets must exist at large (>1-2 AU) distances from their parent stars. As sensitivity and completion of the observed sample increases with time, this distant majority population of giant planets should be revealed. We find that the current distribution of extrasolar giant planet masses implies that high mass (more than 1-2 Jupiter masses) giant planet formation must be relatively rare. Finally, our simulations imply that the efficiency of giant planet formation must be high: at least 10% and perhaps as many as 80% of solar-type stars possess giant planets during their pre-main sequence phase. These predictions, including those for pre-main sequence stars, are testable with the next generation of ground- and space-based planet detection techniques

We present ground-based near-infrared (H-band) imaging of the circumstellar disk around the nearby classical T Tauri star TW Hydrae. The scattered light image shows a face-on disk with radius 4 arcseconds (corresponding to 225 AU) and a morphology which agrees with recent images from the Hubble Space Telescope and the Very Large Array. The best fit power-law for the disk's radial surface brightness profile obeys the law r^{-3.3\pm0.3}. We use our image and published continuum flux densities to derive properties of the disk with a simple model of emission from an optically thick, flat disk. The best-fit values for disk mass and inner radius are 0.03 M_{\odot} and 0.3 AU; best-fit values for temperature, density, and grain opacity power law exponents (q, p, and \beta) are 0.7, 1.3, and 0.9, respectively. These properties are similar to those of disks around classical T Tauri stars located in more distant molecular clouds. Because of TW Hydrae's nearby location and pole-on orientation, it is a uniquely favorable object for future studies of radial disk structure at the classical T Tauri stage.

We discuss the possibility that gaseous giant planets drive strong outflows during early phases of their formation. We consider the range of parameters appropriate for magneto-centrifugally driven stellar and disk outflow models and find that if the proto-Jovian planet or accretion disk had a magnetic field of >~ 10 Gauss and moderate mass inflow rates through the disk of less than 10^-7 M_J/yr that it is possible to drive an outflow. Estimates based both on scaling from empirical laws observed in proto-stellar outflows and the magneto-centrigugal disk and stellar+disk wind models suggest that winds with mass outflow rates of 10^-8 M_J/yr and velocities of order ~ 20 km/s could be driven from proto-Jovian planets. Prospects for detection and some implications for the formation of the solar system are briefly discussed.

Giant planets in circumstellar disks can migrate inward from their initial (formation) positions. Radial migration is caused by inward torques between the planet and the disk; by outward torques between the planet and the spinning star; and by outward torques due to Roche lobe overflow and consequent mass loss from the planet. We present self-consistent numerical considerations of the problem of migrating giant planets. Summing torques on planets for various physical parameters, we find that Jupiter-mass planets can stably arrive and survive at small heliocentric distances, thus reproducing observed properties of some of the recently discovered extra-solar planets. Inward migration timescales can be approximately equal to or less than disk lifetimes and star spindown timescales. Therefore, the range of fates of massive planets is broad, and generally comprises three classes: (I) planets which migrate inward too rapidly and lose all their mass; (II) planets which migrate inward, lose some but not all of their mass, and survive in very small orbits; and (III) planets which do not lose any mass. Some planets in Class III do not migrate very far from their formation locations. Our results show that there is a wide range of possible fates for Jupiter-mass planets for both final heliocentric distance and final mass.