In the outer solar system, the Kuiper Belt contains dynamical sub-populations sculpted by a combination of planet formation and migration and gravitational perturbations from the present-day giant planet configuration. The subdivision of observed Kuiper Belt objects (KBOs) into different dynamical classes is based on their current orbital evolution in numerical integrations of their orbits. Here we demonstrate that machine learning algorithms are a promising tool for reducing both the computational time and human effort required for this classification. Using a Gradient Boosting Classifier, a type of machine learning regression tree classifier trained on features derived from short numerical simulations, we sort observed KBOs into four broad, dynamically distinct populations - classical, resonant, detached, and scattering - with a >97 per cent accuracy for the testing set of 542 securely classified KBOs. Over 80 per cent of these objects have a $>3\sigma$ probability of class membership, indicating that the machine learning method is classifying based on the fundamental dynamical features of each population. We also demonstrate how, by using computational savings over traditional methods, we can quickly derive a distribution of class membership by examining an ensemble of object clones drawn from the observational errors. We find two major reasons for misclassification: inherent ambiguity in the orbit of the object - for instance, an object that is on the edge of resonance - and a lack of representative examples in the training set. This work provides a promising avenue to explore for fast and accurate classification of the thousands of new KBOs expected to be found by surveys in the coming decade.

We investigated the dynamical stability of high-multiplicity Kepler and K2 planetary systems. Our numerical simulations find instabilities in $\sim20\%$ of the cases on a wide range of timescales (up to $5\times10^9$ orbits) and over an unexpectedly wide range of initial dynamical spacings. To identify the triggers of long-term instability in multi-planet systems, we investigated in detail the five-planet Kepler-102 system. Despite having several near-resonant period ratios, we find that mean motion resonances are unlikely to directly cause instability for plausible planet masses in this system. Instead, we find strong evidence that slow inward transfer of angular momentum deficit (AMD) via secular chaos excites the eccentricity of the innermost planet, Kepler-102 b, eventually leading to planet-planet collisions in $\sim80\%$ of Kepler-102 simulations. Kepler-102 b likely has a mass $>\sim0.1M_{\oplus}$, hence a bulk density exceeding about half Earth's, in order to avoid dynamical instability. To investigate the role of secular chaos in our wider set of simulations, we characterize each planetary system's AMD evolution with a "spectral fraction" calculated from the power spectrum of short integrations ($\sim5\times10^6$ orbits). We find that small spectral fractions ($\lesssim0.01$) are strongly associated with dynamical stability on long timescales ($5\times10^9$ orbits) and that the median time to instability decreases with increasing spectral fraction. Our results support the hypothesis that secular chaos is the driver of instabilities in many non-resonant multi-planet systems, and also demonstrate that the spectral analysis method is an efficient numerical tool to diagnose long term (in)stability of multi-planet systems from short simulations.

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

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

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

The composition of comets in the Solar System come in multiple groups thought to encode information about their formation in different regions of the outer protosolar disk. The recent discovery of the second interstellar object, 2I/Borisov, allows for spectroscopic investigations into its gas content and a preliminary classification of it within the Solar System comet taxonomies to test the applicability of planetesimal formation models to other stellar systems. We present spectroscopic and imaging observations from 2019 September 20th to October 26th at the Bok, MMT, and LBT telescopes. We identify CN in the comet's spectrum and set precise upper limits on the abundance of C2 on all dates. We use a Haser model to convert our integrated fluxes to production rates and find Q(CN) = 5.0 +/- 2.0 * 10^24 mol/s on September 20th and Q(CN) = 1.1 - 1.9 * 10^24 mol/s on later dates, both consistent with contemporaneous observations. We set our lowest upper limit on a C2 production rate, Q(C2) < 1.6 * 10^23 mol/s, on October 10th. The measured ratio upper limit for that date, Q(C2)/Q(CN) < 0.095 indicates that 2I/Borisov is strongly in the (carbon chain) 'depleted' taxonomic group. The only comparable Solar System comets have detected ratios near this limit, making 2I/Borisov statistically likely to be more depleted than any known comet. Most 'depleted' comets are Jupiter Family Comets, perhaps indicating a similiarity in formation conditions between the most depleted of the JFCs and 2I/Borisov. More work is needed to understand the applicability of our knowledge of Solar System comet taxonomies onto interstellar objects, and we discuss future work that could help clarify the usefulness of the approach.

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

The Centaurs are the small solar system bodies intermediate between the active inner solar system Jupiter Family Comets and their inactive progenitors in the trans-Neptunian region. Among the fraction of Centaurs which show comet-like activity, 174P/Echeclus is best known for its massive 2005 outburst in which a large apparently active fragment was ejected above the escape velocity from the primary nucleus. We present visible imaging and near-infrared spectroscopy of Echeclus during the first week after its December 2017 outburst taken at the Faulkes North & South Telescopes and the NASA IRTF, the largest outburst since 2005. The coma was seen to be highly asymmetric. A secondary peak was seen in the near-infrared 2D spectra, which is strongly hinted at in the visible images, moving hyperbolically with respect to the nucleus. The retrieved reflectance spectrum of Echelcus is consistent with the unobscured nucleus but becomes bluer when a wider extraction aperture is used. We find that Echeclus's coma is best explained as dominated by large blue dust grains, which agrees with previous work. We also conducted a high-resolution orbital integration of Echeclus's recent evolution and found no large orbital changes that could drive its modern evolution. We interpret the second peak in the visible and near-infrared datasets as a large cloud of larger-than-dust debris ejected at the time of outburst. If Echeclus is typical of the Centaurs, there may be several debris ejection or fragmentation events per year on other Centaurs that are going unnoticed.

Jupiter-family comets (JFCs) are the evolutionary products of trans-Neptunian objects (TNOs) that evolve through the giant planet region as Centaurs and into the inner solar system. Through numerical orbital evolution calculations following a large number of TNO test particles that enter the Centaur population, we have identified a short-lived dynamical Gateway, a temporary low-eccentricity region exterior to Jupiter through which the majority of JFCs pass. We apply an observationally based size distribution function to the known Centaur population and obtain an estimated Gateway region population. We then apply an empirical fading law to the rate of incoming JFCs implied by the the Gateway region residence times. Our derived estimates are consistent with observed population numbers for the JFC and Gateway populations. Currently, the most notable occupant of the Gateway region is 29P/Schwassmann-Wachmann 1 (SW1), a highly active, regularly outbursting Centaur. SW1's present-day, very-low-eccentricity orbit was established after a 1975 Jupiter conjunction and will persist until a 2038 Jupiter conjunction doubles its eccentricity and pushes its semi-major axis out to its current aphelion. Subsequent evolution will likely drive SW1's orbit out of the Gateway region, perhaps becoming one of the largest JFCs in recorded history. The JFC Gateway region coincides with a heliocentric distance range where the activity of observed cometary bodies increases significantly. SW1's activity may be typical of the early evolutionary processing experienced by most JFCs. Thus, the Gateway region, and its most notable occupant SW1, are critical to both the dynamical and physical transition between Centaurs and JFCs.

While collisional families are common in the asteroid belt, only one is known in the Kuiper belt, linked to the dwarf planet Haumea. The characterization of Haumea's family helps to constrain its origin and, more generally, the collisional history of the Kuiper belt. However, the size distribution of the Haumea family is difficult to constrain from the known sample, which is affected by discovery biases. Here, we use the Outer Solar System Origins Survey (OSSOS) Ensemble to look for Haumea family members. In this OSSOS XVI study we report the detection of three candidates with small ejection velocities relative to the family formation centre. The largest discovery, 2013 UQ15, is conclusively a Haumea family member, with a low ejection velocity and neutral surface colours. Although the OSSOS Ensemble is sensitive to Haumea family members to a limiting absolute magnitude (Hr) of 9.5 (inferred diameter of ~90 km), the smallest candidate is significantly larger, Hr = 7.9. The Haumea family members larger than approximately 20 km in diameter must be characterized by a shallow H-distribution slope in order to produce only these three large detections. This shallow size distribution suggests that the family formed in a graze-and-merge scenario, not a catastrophic collision.

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

We present numerical simulations of giant planet migration in our solar system and examine how the speed of planetary migration affects inclinations in the resulting population of small bodies (test particles) scattered outward and subsequently captured into Neptune's 3:2 mean motion resonance (the Plutinos) as well as the hot classical Kuiper belt population. We do not find a consistent relationship between the degree of test particle inclination excitation and e-folding planet migration timescales in the range 5-50 Myr. Our results present a counter-example to Nesvorny 2015's finding that the Plutino and hot classical inclinations showed a marked increase with increasing e-folding timescales for Neptune's migration. We argue that these differing results are likely due to differing secular architectures of the giant planets during and after migration. Small changes in the planets' initial conditions and differences in the numerical implementation of planet migration can result in different amplitudes of the planets' inclination secular modes, and this can lead to different final inclination distributions for test particles in the simulations. We conclude that the observed large inclination dispersion of Kuiper belt objects does not require Neptune's migration to be slow; planetary migration with e-folding timescales of 5, 10, 30, and 50 Myr can all yield inclination dispersions similar to the observed Plutino and hot classical populations, with no correlation between the degree of inclination excitation and migration speed.

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

We present an analysis and comparison of the 30 micron dust features seen in the Spitzer Space Telescope spectra of 207 carbon-rich asymptotic giant branch (AGB) stars, post-AGB objects, and planetary nebulae located in the Milky Way, the Magellanic Clouds (MCs), or the Sagittarius dwarf spheroidal galaxy (Sgr dSph), which are characterised by different average metallicities. We investigated whether the formation of the 30 micron feature carrier may be a function of the metallicity. Through this study we expect to better understand the late stages of stellar evolution of carbon-rich stars in these galaxies. Our analysis uses the `Manchester method' as a basis for estimating the temperature of dust for the carbon-rich AGB stars and the planetary nebulae in our sample. We used a black-body function with a single temperature deduced from the Manchester method or its modification to approximate the continuum under the 30 micron feature. The most important conclusion of our work is the fact that the formation of the 30 micron feature is affected by metallicity. Specifically that, as opposed to more metal-poor samples of AGB stars in the MCs, the feature is seen at lower mass-loss rates, higher temperatures, and has seen to be more prominent in Galactic carbon stars. The averaged feature (profile) in the AGB, post-AGB objects, and PNe seems unaffected by metallicity at least between a fifth and solar metallicity, but in the case of PNe it is shifted to significantly longer wavelengths.

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

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

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

Making an inventory of the Solar System is one of the four fundamental science requirements for the Large Synoptic Survey Telescope (LSST). The current baseline footprint for LSST's main Wide-Fast-Deep (WFD) Survey observes the sky below 0$^\circ$ declination, which includes only half of the ecliptic plane. Critically, key Solar System populations are asymmetrically distributed on the sky: they will be entirely missed, or only partially mapped, if only the WFD occurs. We propose a Northern Ecliptic Spur (NES) mini survey, observing the northern sky up to +10$^\circ$ ecliptic latitude, to maximize Solar System science with LSST. The mini survey comprises a total area of $\sim$5800 deg$^2$/604 fields, with 255 observations/field over the decade, split between g,r, and z bands. Our proposed survey will 1) obtain a census of main-belt comets; 2) probe Neptune's past migration history, by exploring the resonant structure of the Kuiper belt and the Neptune Trojan population; 3) explore the origin of Inner Oort cloud objects and place significant constraints on the existence of a hypothesized planet beyond Neptune; and 4) enable precise predictions of KBO stellar occultations. These high-ranked science goals of the Solar System Science Collaboration are only achievable with this proposed northern survey.

Making an inventory of the Solar System is one of the four pillars that the requirements for the Large Synoptic Survey Telescope (LSST) are built upon. The choice between same-filter nightly pairs or different-filter nightly pairs in the Wide-Fast-Deep (WFD) Survey will have a dramatic effect on the ability of the Moving Object Pipeline System (MOPS) to detect certain classes of Solar System objects; many of the possible filter pairings would result in significant ($\sim50\%$ or more) loss of Solar System object detections. In particular, outer Solar System populations can be significantly redder than those in the inner Solar System, and nightly pairs in $r$-band will result in the deepest survey for the outer Solar System. To maximize the potential for outer Solar System science, we thus advocate for ensuring that the WFD survey contains a sufficient number of $r$-$r$ nightly pairs for each field during a discovery season to ensure detection and linking using MOPS. We also advocate for adding additional spectral energy distributions (SEDs) that more accurately model outer Solar System populations to the pipeline for evaluating the outputs of the LSST operations simulator. This will enable a better estimate of how many Solar System population detections are lost or gained for different filter choices in the WFD survey.

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

Observations of Kuiper belt objects (KBOs) in Neptune's 5:2 resonance present two puzzles: this third order resonance hosts a surprisingly large population, comparable to the prominent populations of Plutinos and Twotinos in the first order 3:2 and 2:1 resonances, respectively; secondly, their eccentricities are concentrated near $0.4$. To shed light on these puzzles, we investigate the phase space near this resonance with use of Poincar\'e sections of the circular planar restricted three body model. We find several transitions in the phase space structure with increasing eccentricity, which we explain with the properties of the resonant orbit relative to Neptune's. The resonance width is narrow for very small eccentricities, but widens dramatically for $e\gtrsim0.2$, reaching a maximum near $e\approx0.4$, where it is similar to the maximum widths of the 2:1 and 3:2 resonances. We confirm these results with N-body numerical simulations, including the effects of all four giant planets and a wide range of orbital inclinations of the KBOs. We find that the boundaries of the stable resonance zone are not strongly sensitive to inclination and remain very similar to those found with the simplified three body model, with the caveat that orbits of eccentricity above $\sim0.53$ are unstable; higher eccentricity orbits are phase-protected from destabilizing encounters with Neptune but not with Uranus. These results show that the 5:2 resonant KBOs are not more puzzling than the Plutinos and Twotinos; however, detailed understanding of the origins of eccentric, inclined resonant KBOs remains a challenge.

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

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

A substantial fraction of our solar system's trans-Neptunian objects (TNOs) are in mean motion resonance with Neptune. Many of these objects were likely caught into resonances by planetary migration---either smooth or stochastic---approximately 4 Gyr ago. Some, however, gravitationally scattered off of Neptune and became transiently stuck in more recent events. Here, we use numerical simulations to predict the number of transiently-stuck objects, captured from the current actively scattering population, that occupy 111 resonances at semimajor axes $a=$30--100 au. Our source population is an observationally constrained model of the currently-scattering TNOs. We predict that, integrated across all resonances at these distances, the current transient sticking population comprises 40\% of total transiently-stuck+scattering TNOs, suggesting that these objects should be treated as a single population. We compute the relative distribution of transiently-stuck objects across all $p$:$q$ resonances with $1/6 \le q/p < 1$, $p<40$, and $q<20$, providing predictions for the population of transient objects with $H_r < 8.66$ in each resonance. We find that the relative populations are approximately proportional to each resonance's libration period and confirm that the importance of transient sticking increases with semimajor axis in the studied range. We calculate the expected distribution of libration amplitudes for stuck objects and demonstrate that observational constraints indicate that both the total number and the amplitude-distribution of 5:2 resonant TNOs are inconsistent with a population dominated by transient sticking from the current scattering disk. The 5:2 resonance hence poses a challenge for leading theories of Kuiper belt sculpting.

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

We measured the mean plane of the Kuiper belt as a function of semi-major axis. For the classical Kuiper belt as a whole (the non-resonant objects in the semi-major axis range 42--48~au), we find a mean plane of inclination $i_m=1.8^{\circ}$$^{+0.7^{\circ}}_{-0.4^{\circ}}$ and longitude of ascending node $\Omega_m=77^{\circ}$$^{+18^{\circ}}_{-14^{\circ}}$ (in the J2000 ecliptic-equinox coordinate system), in accord with theoretical expectations of the secular effects of the known planets. With finer semi-major axis bins, we detect a statistically significant warp in the mean plane near semi-major axes 40--42~au. Linear secular theory predicts a warp near this location due to the $\nu_{18}$ nodal secular resonance, however the measured mean plane for the 40.3-42~au semi-major axis bin (just outside the $\nu_{18}$) is inclined $\sim13^{\circ}$ to the predicted plane, a nearly 3-$\sigma$ discrepancy. For the more distant Kuiper belt objects of semi-major axes in the range 50--80~au, the expected mean plane is close to the invariable plane of the solar system, but the measured mean plane deviates greatly from this: it has inclination $i_m=9.1^{\circ}$$^{+6.6^{\circ}}_{-3.8^{\circ}}$ and longitude of ascending node $\Omega_m=227^{\circ}$$^{+18^{\circ}}_{-44^{\circ}}$. We estimate this deviation from the expected mean plane to be statistically significant at the $\sim97-99\%$ confidence level. We discuss several possible explanations for this deviation, including the possibility that a relatively close-in ($a\lesssim100$~au), unseen small planetary-mass object in the outer solar system is responsible for the warping.

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

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

The Infrared Spectrograph (IRS) on the {\em Spitzer Space Telescope} observed nearly 800 point sources in the Large Magellanic Cloud (LMC), taking over 1,000 spectra. 197 of these targets were observed as part of the Sage-Spec Spitzer Legacy program; the remainder are from a variety of different calibration, guaranteed time and open time projects. We classify these point sources into types according to their infrared spectral features, continuum and spectral energy distribution shape, bolometric luminosity, cluster membership, and variability information, using a decision-tree classification method. We then refine the classification using supplementary information from the astrophysical literature. We find that our IRS sample is comprised substantially of YSO and H\,{\sc ii} regions, post-Main Sequence low-mass stars: (post-)AGB stars and planetary nebulae and massive stars including several rare evolutionary types. Two supernova remnants, a nova and several background galaxies were also observed. We use these classifications to improve our understanding of the stellar populations in the Large Magellanic Cloud, study the composition and characteristics of dust species in a variety of LMC objects, and to verify the photometric classification methods used by mid-IR surveys. We discover that some widely-used catalogues of objects contain considerable contamination and others are missing sources in our sample.

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

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

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

The four longest period Kuiper belt objects have orbital periods close to integer ratios with each other. A hypothetical planet with orbital period $\sim$17,117 years, semimajor axis $\sim$665 AU, would have N/1 and N/2 period ratios with these four objects. The orbital geometries and dynamics of resonant orbits constrain the orbital plane, the orbital eccentricity and the mass of such a planet, as well as its current location in its orbital path.

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

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

The Infrared Spectrograph on the Spitzer Space Telescope observed 184 carbon stars in the Magellanic Clouds. This sample reveals that the dust-production rate (DPR) from carbon stars generally increases with the pulsation period of the star. The composition of the dust grains follows two condensation sequences, with more SiC condensing before amorphous carbon in metal-rich stars, and the order reversed in metal-poor stars. MgS dust condenses in optically thicker dust shells, and its condensation is delayed in more metal-poor stars. Metal-poor carbon stars also tend to have stronger absorption from C2H2 at 7.5 um. The relation between DPR and pulsation period shows significant apparent scatter, which results from the initial mass of the star, with more massive stars occupying a sequence parallel to lower-mass stars, but shifted to longer periods. Accounting for differences in the mass distribution between the carbon stars observed in the Small and Large Magellanic Clouds reveals a hint of a subtle decrease in the DPR at lower metallicities, but it is not statistically significant. The most deeply embedded carbon stars have lower variability amplitudes and show SiC in absorption. In some cases they have bluer colors at shorter wavelengths, suggesting that the central star is becoming visible. These deeply embedded stars may be evolving off of the AGB and/or they may have non-spherical dust geometries.

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

The Magellanic clouds are uniquely placed to study the stellar contribution to dust emission. Individual stars can be resolved in these systems even in the mid-infrared, and they are close enough to allow detection of infrared excess caused by dust.We have searched the Spitzer Space Telescope data archive for all Infrared Spectrograph (IRS) staring-mode observations of the Small Magellanic Cloud (SMC) and found that 209 Infrared Array Camera (IRAC) point sources within the footprint of the Surveying the Agents of Galaxy Evolution in the Small Magellanic Cloud (SAGE-SMC) Spitzer Legacy programme were targeted, within a total of 311 staring mode observations. We classify these point sources using a decision tree method of object classification, based on infrared spectral features, continuum and spectral energy distribution shape, bolometric luminosity, cluster membership and variability information. We find 58 asymptotic giant branch (AGB) stars, 51 young stellar objects (YSOs), 4 post-AGB objects, 22 Red Supergiants (RSGs), 27 stars (of which 23 are dusty OB stars), 24 planetary nebulae (PNe), 10Wolf-Rayet (WR) stars, 3 Hii regions, 3 R Coronae Borealis (R CrB) stars, 1 Blue Supergiant and 6 other objects, including 2 foreground AGB stars. We use these classifications to evaluate the success of photometric classification methods reported in the literature.

The Kepler mission results indicate that systems of tightly-packed inner planets (STIPs) are present around of order 5% of FGK field stars (whose median age is ~5 Gyr). We propose that STIPs initially surrounded nearly all such stars and those observed are the final survivors of a process in which long-term metastability eventually ceases and the systems proceed to collisional consolidation or destruction, losing roughly equal fractions of systems every decade in time. In this context, we also propose that our Solar System initially contained additional large planets interior to the current orbit of Venus, which survived in a metastable dynamical configuration for 1-10% of the Solar System's age. Long-term gravitational perturbations caused the system to orbit cross, leading to a cataclysmic event which left Mercury as the sole surviving relic.

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

We have investigated the light variability in a sample of 22 carbon-rich post-AGB stars in the Large Magellanic Cloud (LMC) and Small Magellanic Cloud (SMC), based primarily on photometric data from the OGLE survey. All are found to vary. Dominant periods are found in eight of them; these periods range from 49 to 157 days, and most of these stars have F spectral types. These eight are found to be similar to the Milky Way Galaxy (MWG) carbon-rich proto-planetary nebulae (PPNs) in several ways: (a) they are in the same period range of ~38 to ~160 days, (b) they have similar spectral types, (c) they are (all but one) redder when fainter, (d) they have multiple periods, closely spaced in time, with a average ratio of secondary to primary period of ~1.0, and as an ensemble, (e) they show a trend of decreasing period with increasing temperature, and (f) they show a trend of decreasing amplitude with decreasing period. However, they possibly differ in that the decreasing trend of period with temperature may be slightly offset from that of the MWG. These eight are classified as PPNs. The other 14 all show evidence of variability on shorter timescales. They are likely hotter PPNs or young planetary nebulae. However, in the MWG the numbers of PPNs peak in the F-G spectral types, while it appears that in the LMC they peak at a hotter B spectral type. One of the periodic ones shows a small, R Coronae Borealis-type light curve drop.

When modelling an ionised plasma, all spectral synthesis codes need the thermally averaged free-free Gaunt factor defined over a very wide range of parameter space in order to produce an accurate prediction for the spectrum. Until now no data set exists that would meet these needs completely. We have therefore produced a table of relativistic Gaunt factors over a much wider range of parameter space than has ever been produced before. We present tables of the thermally averaged Gaunt factor covering the range log10(gamma^2) = -6 to 10 and log10(u) = -16 to 13 for all atomic numbers Z = 1 through 36. The data were calculated using the relativistic Bethe-Heitler-Elwert (BHE) approximation and were subsequently merged with accurate non-relativistic results in those parts of the parameter space where the BHE approximation is not valid. These data will be incorporated in the next major release of the spectral synthesis code Cloudy. We also produced tables of the frequency integrated Gaunt factor covering the parameter space log10(gamma^2) = -6 to 10 for all values of Z between 1 and 36. All the data presented in this paper are available online.

Modern spectral synthesis codes need the thermally averaged free-free Gaunt factor defined over a very wide range of parameter space in order to produce an accurate prediction for the spectrum emitted by an ionized plasma. Until now no set of data exists that would meet this need in a fully satisfactory way. We have therefore undertaken to produce a table of very accurate non-relativistic Gaunt factors over a much wider range of parameters than has ever been produced before. We first produced a table of non-averaged Gaunt factors, covering the parameter space log10(epsilon_i) = -20 to +10 and log10(w) = -30 to +25. We then continued to produce a table of thermally averaged Gaunt factors covering the parameter space log10(gamma^2) = -6 to +10 and log10(u) = -16 to +13. Finally we produced a table of the frequency integrated Gaunt factor covering the parameter space log10(gamma^2) = -6 to +10. All the data presented in this paper are available online.

K-band spectra are presented for a sample of 39 Spitzer IRS SAGE-Spec sources in the Large Magellanic Cloud. The spectra exhibit characteristics in very good agreement with their positions in the near infrared - Spitzer color-magnitude diagrams and their properties as deduced from the Spitzer IRS spectra. Specifically, the near infrared spectra show strong atomic and molecular features representative of oxygen-rich and carbon-rich asymptotic giant branch stars, respectively. A small subset of stars were chosen from the luminous and red extreme "tip" of the color magnitude diagram. These objects have properties consistent with dusty envelopes but also cool, carbon-rich "stellar" cores. Modest amounts of dust mass loss combine with the stellar spectral energy distribution to make these objects appear extreme in their near infrared and mid infrared colors. One object in our sample, HV 915, a known post asymptotic giant branch star of the RV Tau type exhibits CO 2.3 micron band head emission consistent with previous work that demonstrates the object has a circumstellar disk.

Infrared spectra of carbon-rich objects which have evolved off the asymptotic giant branch reveal a range of dust properties, including fullerenes, polycyclic aromatic hydrocarbons (PAHs), aliphatic hydrocarbons, and several unidentified features, including the 21 um emission feature. To test for the presence of fullerenes, we used the position and width of the feature at 18.7-18.9 um and examined other features at 17.4 and 6-9 um. This method adds three new fullerene sources to the known sample, but it also calls into question three previous identifications. We confirm that the strong 11 um features seen in some sources arise primarily from SiC, which may exist as a coating around carbonaceous cores and result from photo-processing. Spectra showing the 21 um feature usually show the newly defined Class D PAH profile at 7-9 um. These spectra exhibit unusual PAH profiles at 11-14 um, with weak contributions at 12.7 um, which we define as Class D1, or show features shifted to ~11.4, 12.4, and 13.2 um, which we define as Class D2. Alkyne hydrocarbons match the 15.8 um feature associated with 21 um emission. Sources showing fullerene emission but no PAHs have blue colors in the optical, suggesting a clear line of sight to the central source. Spectra with 21 um features and Class D2 PAH emission also show photometric evidence for a relatively clear line of sight to the central source. The multiple associations of the 21 um feature to aliphatic hydrocarbons suggest that the carrier is related to this material in some way.

JWST's Near-Infrared Imager and Slitless Spectrograph (NIRISS) includes an Aperture Masking Interferometry (AMI) mode designed to be used between 2.7{\mu}m and 4.8{\mu}m. At these wavelengths, it will have the highest angular resolution of any mode on JWST, and, for faint targets, of any existing or planned infrastructure. NIRISS AMI is uniquely suited to detect thermal emission of young massive planets and will permit the characterization of the mid-IR flux of exoplanets discovered by the GPI and SPHERE adaptive optics surveys. It will also directly detect massive planets found by GAIA through astrometric accelerations, providing the first opportunity ever to get both a mass and a flux measurement for non-transiting giant planets. NIRISS AMI will also enable the study of the nuclear environment of AGNs.

This paper reports variations of polycyclic aromatic hydrocarbons (PAHs) features that were found in Spitzer Space Telescope spectra of carbon-rich post-asymptotic giant branch (post-AGB) stars in the Large Magellanic Cloud (LMC). The paper consists of two parts. The first part describes our Spitzer spectral observing programme of 24 stars including post-AGB candidates. The latter half of this paper presents the analysis of PAH features in 20 carbon-rich post-AGB stars in the LMC, assembled from the Spitzer archive as well as from our own programme. We found that five post-AGB stars showed a broad feature with a peak at 7.7 micron, that had not been classified before. Further, the 10--13 micron PAH spectra were classified into four classes, one of which has three broad peaks at 11.3, 12.3 and 13.3 micron rather than two distinct sharp peaks at 11.3 and 12.7 micron, as commonly found in HII regions. Our studies suggest that PAHs are gradually processed while the central stars evolve from post-AGB phase to PNe, changing their composition before PAHs are incorporated into the interstellar medium. Although some metallicity dependence of PAH spectra exists, the evolutionary state of an object is more significant than its metallicity in determining the spectral characteristics of PAHs for LMC and Galactic post-AGB stars.

The Centaurs are a population of small, planet-crossing objects in the outer solar system. They are dynamically short-lived and represent the transition population between the Kuiper belt and the Jupiter family short-period comets. Dynamical models and observations of the physical properties of the Centaurs indicate that they may have multiple source populations in the trans-Neptunian region. It has been suggested that the inclination distribution of the Centaurs may be useful in distinguishing amongst these source regions. The Centaurs, however, undergo many close encounters with the giant planets during their orbital evolution; here we show that these encounters can substantially determine the inclination distribution of the Centaurs. Almost any plausible initial inclination distribution of a Kuiper belt source results in Centaurs having inclinations peaked near 10-20 degrees. Our studies also find that the Kuiper belt is an extremely unlikely source of the retrograde Centaur that has been observed.

The Haumea family is currently the only identified collisional family in the Kuiper belt. We numerically simulate the long-term dynamical evolution of the family to estimate a lower limit of the family's age and to assess how the population of the family and its dynamical clustering are preserved over Gyr timescales. We find that the family is not younger than 100 Myr, and its age is at least 1 Gyr with 95% confidence. We find that for initial velocity dispersions of 50-400 m/s, approximately 20-45% of the family members are lost to close encounters with Neptune after 3.5 Gyr of orbital evolution. We apply these loss rates to two proposed models for the formation of the Haumea family, a graze-and-merge type collision between two similarly sized, differentiated KBOs or the collisional disruption of a satellite orbiting Haumea. For the graze-and-merge collision model, we calculate that >85% of the expected mass in surviving family members within 150 m/s of the collision has been identified, but that one to two times the mass of the known family members remains to be identified at larger velocities. For the satellite-break-up model, we estimate that the currently identified family members account for ~50% of the expected mass of the family. Taking observational incompleteness into account, the observed number of Haumea family members is consistent with either formation scenario at the 1 sigma level, however both models predict more objects at larger relative velocities (>150 m/s) than have been identified.

Two types of dust disks around white dwarfs (WDs) have been reported: small dust disks around cool metal-rich WDs consisting of tidally disrupted asteroids, and a large dust disk around the hot central WD of the Helix planetary nebula (PN) possibly produced by collisions among Kuiper Belt-like objects. To search for more dust disks of the latter type, we have conducted a Spitzer MIPS 24 um survey of 71 hot WDs or pre-WDs, among which 35 are central stars of PNe (CSPNs). Nine of these evolved stars are detected and their 24 um flux densities are at least two orders of magnitude higher than their expected photospheric emission. Considering the bias against detection of distant objects, the 24 um detection rate for the sample is >~15%. It is striking that seven, or ~20%, of the WD and pre-WDs in known PNe exhibit 24 um excesses, while two, or 5-6%, of the WDs not in PNe show 24 um excesses and they have the lowest 24 um flux densities. We have obtained follow-up Spitzer IRS spectra for five objects. Four show clear continuum emission at 24 um, and one is overwhelmed by a bright neighboring star but still show a hint of continuum emission. In the cases of WD 0950+139 and CSPN K1-22, a late-type companion is present, making it difficult to determine whether the excess 24 um emission is associated with the WD or its red companion. High-resolution images in the mid-IR are needed to establish unambiguously the stars responsible for the 24 um excesses.

We investigate the long-term evolution of the inclinations of the known classical and resonant Kuiper belt objects (KBOs). This is partially motivated by the observed bimodal inclination distribution and by the putative physical differences between the low- and high-inclination populations. We find that some classical KBOs undergo large changes in inclination over gigayear timescales, which means that a current member of the low-inclination population may have been in the high-inclination population in the past, and vice versa. The dynamical mechanisms responsible for the time-variability of inclinations are predominantly distant encounters with Neptune and chaotic diffusion near the boundaries of mean motion resonances. We reassess the correlations between inclination and physical properties including inclination time-variability. We find that the size-inclination and color-inclination correlations are less statistically significant than previously reported (mostly due to the increased size of the data set since previous works with some contribution from inclination variability). The time-variability of inclinations does not change the previous finding that binary classical KBOs have lower inclinations than non-binary objects. Our study of resonant objects in the classical Kuiper belt region includes objects in the 3:2, 7:4, 2:1, and eight higher-order mean motion resonances. We find that these objects (some of which were previously classified as non-resonant) undergo larger changes in inclination compared to the non-resonant population, indicating that their current inclinations are not generally representative of their original inclinations. They are also less stable on gigayear timescales.

Low- and intermediate-mass stars are one of the important dust sources in the interstellar medium (ISM) of galaxies. The compositions of dust ejected from these stars are likely to affect those in the ISM. We investigate dust in post-Asymptotic Giant Branch (AGB) stars, which are in a late evolutionary phase for low- and intermediate-mass stars, and which produce a wide variety of dust grains. We are particularly targeting post-AGB stars in the Large Magellanic Cloud (LMC), which has about half of the solar metallicity, to investigate the effects of sub-solar metallicity on dust compositions. Using the Spitzer Space Telescope, we obtained 5-30 micron spectra of 24 post-AGB candidates in the LMC. Five are C-rich post-AGB stars, and this presentation focuses on spectra of these stars. We found that rare dust features in the Milky Way, such as a 21 micron unidentified feature are commonly found in LMC post-AGB stars. The 6-8 micron spectra are compared with those of Galactic objects. Four spectra match the Galactic templates of polycyclic aromatic hydrocarbon (PAH) features. However, we found the three objects show 7.85 micron feature which have not found in Galactic post-AGB stars. Low metallicity conditions definitely affect the dust formation process and compositions.

High-resolution HST imaging of the compact planetary nebula NGC 6644 has revealed two pairs of bipolar lobes and a central ring lying close to the plane of the sky. From mid-infrared imaging obtained with the Gemini Telescope, we have found a dust torus which is oriented nearly perpendicular to one pair of the lobes. We suggest that NGC 6644 is a multipolar nebula and have constructed a 3-D model which allows the visualization of the object from different lines of sight. These results suggest that NGC 6644 may have similar intrinsic structures as other multipolar nebulae and the phenomenon of multipolar nebulosity may be more common than previously believed.

We present the classification of 197 point sources observed with the Infrared Spectrograph in the SAGE-Spec Legacy program on the Spitzer Space Telescope. We introduce a decision-tree method of object classification based on infrared spectral features, continuum and spectral energy distribution shape, bolometric luminosity, cluster membership, and variability information, which is used to classify the SAGE-Spec sample of point sources. The decision tree has a broad application to mid-infrared spectroscopic surveys, where supporting photometry and variability information are available. We use these classifications to make deductions about the stellar populations of the Large Magellanic Cloud and the success of photometric classification methods. We find 90 asymptotic giant branch (AGB) stars, 29 young stellar objects, 23 post-AGB objects, 19 red supergiants, eight stellar photospheres, seven background galaxies, seven planetary nebulae, two HII regions and 12 other objects, seven of which remain unclassified.

We present a 2Dust model for the dust shell around a LMC long-period variable (LPV) previously studied as part of the OGLE survey. OGLE LMC LPV 28579 (SAGE J051306.40-690946.3) is a carbon-rich asymptotic giant branch (AGB) star for which we have photometry and spectra from the Spitzer SAGE and SAGE-Spec programs along with UBVIJHK_s photometry. By modeling this source, we obtain a baseline set of dust properties to be used in the construction of a grid of models for carbon stars. We reproduce its spectral energy distribution using a mixture of AmC and SiC (15% by mass). The grain sizes are distributed according to the KMH model. The best-fit model has an optical depth of 0.28 for the shell at the peak of the SiC feature, with R_in~1430 R_sun or 4.4 R_star. The temperature at this inner radius is 1310 K. Assuming an expansion velocity of 10 km s^-1, we obtain a dust mass-loss rate of 2.5x10^-9 M_sun yr-1. We calculate a 15% variation in this rate by testing the fit sensitivity against variation in input parameters. We also present a simple model for the molecular gas in the extended atmosphere that could give rise to the 13.7 \mu m feature seen in the spectrum. We find that a combination of CO and C_2H_2 gas at an excitation temperature of about 1000 K and column densities of 3x10^21 cm^-2 and 10^19 cm^-2 respectively are able to reproduce the observations. Given that the excitation temperature is close to T_dust(R_in), most of the molecular contribution probably arises from the inner shell region. The luminosity corresponding to the first epoch of SAGE observations is 6580 L_sun. For an effective temperature of about 3000 K, this implies a stellar mass of 1.5-2 M_sun and an age of 1-2.5 Gyr. For a gas:dust ratio of 200, we obtain a gas mass-loss rate of 5.0x10^-7 M_sun yr^-1, consistent with the gas mass-loss rates estimated from the period, color and 8 \mu m flux of the source.

The SAGE-Spec Spitzer Legacy program is a spectroscopic follow-up to the SAGE-LMC photometric survey of the Large Magellanic Cloud carried out with the Spitzer Space Telescope. We present an overview of SAGE-Spec and some of its first results. The SAGE-Spec program aims to study the life cycle of gas and dust in the Large Magellanic Cloud, and to provide information essential to the classification of the point sources observed in the earlier SAGE-LMC photometric survey. We acquired 224.6 hours of observations using the InfraRed Spectrograph and the SED mode of the Multiband Imaging Photometer for Spitzer. The SAGE-Spec data, along with archival Spitzer spectroscopy of objects in the Large Magellanic Cloud, are reduced and delivered to the community. We discuss the observing strategy, the specific data reduction pipelines applied and the dissemination of data products to the scientific community. Initial science results include the first detection of an extragalactic "21 um" feature towards an evolved star and elucidation of the nature of disks around RV Tauri stars in the Large Magellanic Cloud. Towards some young stars, ice features are observed in absorption. We also serendipitously observed a background quasar, at a redshift of z~0.14, which appears to be host-less.

The "Frankfurter Neutronenquelle am Stern-Gerlach-Zentrum" (FRANZ), which is currently under development, will be the strongest neutron source in the astrophysically interesting energy region in the world. It will be about three orders of magnitude more intense than the well-established neutron source at the Research Center Karlsruhe (FZK).

Spitzer MIPS 24 um images were obtained for 36 Galactic planetary nebulae (PNe) whose central stars are hot white dwarfs (WDs) or pre-WDs with effective temperatures of ~100,000 K or higher. Diffuse 24 um emission is detected in 28 of these PNe. The eight non-detections are angularly large PNe with very low H-alpha surface brightnesses. We find three types of correspondence between the 24 um emission and H-alpha line emission of these PNe: six show 24 um emission more extended than H-alpha emission, nine have a similar extent at 24 um and H-alpha, and 13 show diffuse 24 um emission near the center of the H-alpha shell. The sizes and surface brightnesses of these three groups of PNe and the non-detections suggest an evolutionary sequence, with the youngest ones being brightest and the most evolved ones undetected. The 24 um band emission from these PNe is attributed to [O IV] 25.9 um and [Ne V] 24.3 um line emission and dust continuum emission, but the relative contributions of these three components depend on the temperature of the central star and the distribution of gas and dust in the nebula.

Using the data obtained with the Spitzer Space telescope as part of the Surveying the Agents of a Galaxy's Evolution (SAGE) legacy survey, we have studied the variations of the dust composition and abundance across the Large Magellanic Cloud (LMC). Such variations are expected, as the explosive events which have lead to the formation of the many HI shells observed should have affected the dust properties. Using a model and comparing with a reference spectral energy distribution from our Galaxy, we deduce the relative abundance variations of small dust grains across the LMC. We examined the infrared color ratios as well as the relative abundances of very small grains (VSGs) and polycyclic aromatic hydrocarbons (PAHs) relative to the big grain (BG) abundance. Results show that each dust component could have different origins or evolution in the interstellar medium (ISM). The VSG abundance traces the star formation activity and could result from shattering of larger grains, whereas the PAH abundance increases around molecular clouds as well as in the stellar bar, where they could have been injected into the ISM during mass loss from old stars.

We report on an analysis of the gas and dust budget in the the interstellar medium (ISM) of the Large Magellanic Cloud (LMC). Recent observations from the Spitzer Space Telescope enable us to study the mid-infrared dust excess of asymptotic giant branch (AGB) stars in the LMC. This is the first time we can quantitatively assess the gas and dust input from AGB stars over a complete galaxy, fully based on observations. The integrated mass-loss rate over all intermediate and high mass-loss rate carbon-rich AGB candidates in the LMC is 8.5x10^-3 solar mass per year, up to 2.1x10^-2 solar mass per year. This number could be increased up to 2.7x10^-2 solar mass per year, if oxygen-rich stars are included. This is overall consistent with theoretical expectations, considering the star formation rate when these low- and intermediate-mass stars where formed, and the initial mass functions. AGB stars are one of the most important gas sources in the LMC, with supernovae (SNe), which produces about 2-4x10^-2 solar mass per year. At the moment, the star formation rate exceeds the gas feedback from AGB stars and SNe in the LMC, and the current star formation depends on gas already present in the ISM. This suggests that as the gas in the ISM is exhausted, the star formation rate will eventually decline in the LMC, unless gas is supplied externally. Our estimates suggest `a missing dust-mass problem' in the LMC, which is similarly found in high-z galaxies: the accumulated dust mass from AGB stars and possibly SNe over the dust life time (400--800 Myrs) is significant less than the dust mass in the ISM. Another dust source is required, possibly related to star-forming regions.

We present empirical relations describing excess emission from evolved stars in the Large Magellanic Cloud (LMC) using data from the SAGE (Surveying the Agents of a Galaxy's Evolution) survey which includes the IRAC 3.6, 4.5, 5.8 and 8.0 \mu m and MIPS 24, 70 and 160 \mu m bands. We combine the SAGE data with the Two Micron All Sky Survey (2MASS; J, H and Ks) and the optical Magellanic Cloud Photometric Survey (MCPS; U, B, V and I) point source catalogs to create complete spectral energy distributions (SEDs) of the asymptotic giant branch (AGB) star candidates in the LMC. AGB star outflows are among the main producers of dust in a galaxy, and this mass loss results in an excess in the fluxes observed in the 8 and 24 \mic m bands. We identify oxygen-rich, carbon-rich and extreme AGB star populations in our sample based on their 2MASS and IRAC colors. We calculate excesses fluxes in the mid-IR bands by comparison of the SEDs with model photospheres. We find about 16,000 O-rich, 6300 C-rich and 1000 extreme sources with reliable 8 \mu m excesses, and about 4500 O-rich, 5300 C-rich and 960 extreme sources with reliable 24 \mic m excesses. The excesses are in the range 0.1 mJy--5 Jy. The 8 and 24 \mic m excesses for all three types of AGB candidates show a general increasing trend with luminosity. The dust color temperature derived from the ratio of the 8 and 24 \mic m excesses decreases with an increase in excess, while the 24 \mic m optical depth increases with excess. The extreme AGB candidates are the major contributors to the mass loss, and we estimate the total AGB mass-loss return to the LMC to be (5.9--13)\times 10^{-3} Msun yr^{-1}.

We have carried out mid-infrared spectroscopy of seven Galactic proto-planetary nebulae (PPNs) using the Spitzer Space Telescope. They were observed from 10-36 microns at relatively high spectral resolution, R~600. The sample was chosen because they all gave some evidence in the visible of a carbon-rich chemistry. All seven of the sources show the broad, unidentified 21 micron emission feature; three of them are new detections (IRAS 06530-0213, 07430+1115, and 19477+2401) and the others are observed at higher S/N than in previous spectra. These have the same shape and central wavelength (20.1 microns) as found in the ISO spectra of the brighter PPNs. The 30 micron feature was seen in all seven objects. However, it is not resolved into two separate features (26 and 33 microns) as was claimed on the basis of ISO spectra, which presumably suffered from the noisy detector bands in this region. All showed the infrared aromatic bands (AIB) at 11.3, 12.4, and 13.3 microns. Five of these also appear to have the C2H2 molecular band at 13.7 microns, one in absorption and four in emission. This is extremely rare, with only one other evolved star, IRC+10216, in which C2H2 emission has been observed. Four also possessed a broad, unidentified emission feature at 15.8 microns that may possibly be related to the 21 micron feature. Model fits were made to the spectral energy distributions for these PPNs to determine properties of the detached circumstellar envelopes. The 21 micron feature has been seen in all Galactic carbon-rich PPNs observed, and thus its carrier appears to be a common component of the outflow around these objects.

The short period Jupiter family comets (JFCs) are thought to originate in the Kuiper Belt; specifically, a dynamical subclass of the Kuiper Belt known as the `scattered disk' is argued to be the dominant source of JFCs. However, the best estimates from observational surveys indicate that this source may fall short by more than two orders of magnitude the estimates obtained from theoretical models of the dynamical evolution of Kuiper belt objects into JFCs. We re-examine the scattered disk as a source of the JFCs and make a rigorous estimate of the discrepancy. We find that the uncertainties in the dynamical models combined with a change in the size distribution function of the scattered disk at faint magnitudes (small sizes) beyond the current observational limit offer a possible but problematic resolution to the discrepancy. We discuss several other possibilities: that the present population of JFCs is a large fluctuation above their long term average, that larger scattered disk objects tidally break-up into multiple fragments during close planetary encounters as their orbits evolve from the trans-Neptune zone to near Jupiter, or that there are alternative source populations that contribute significantly to the JFCs. Well-characterized observational investigations of the Centaurs, objects that are transitioning between the trans-Neptune Kuiper belt region and the inner solar system, can test the predictions of the non-steady state and the tidal break-up hypotheses. The classical and resonant classes of the Kuiper belt are worth re-consideration as significant additional or alternate sources of the JFCs.

We present diffraction limited, 10um imaging polarimetry data for the central regions of the archetypal Seyfert AGN, NGC1068. The position angle of polarization is consistent with three dominant polarizing mechanisms. We identify three distinct regions of polarization: (a) north of the nucleus, arising from aligned dust in the NLR, (b) south, east and west of the nucleus, consistent with dust being channeled toward the central engine and (c) a central minimum of polarization consistent with a compact (<22pc) torus. These observations provide continuity between the geometrically and optically thick torus and the host galaxy's nuclear environments. These images represent the first published mid-IR polarimetry from an 8-m class telescope and illustrate the potential of such observations.

Molecular oxygen, O2 has been expected historically to be an abundant component of the chemical species in molecular clouds and, as such, an important coolant of the dense interstellar medium. However, a number of attempts from both ground and from space have failed to detect O2 emission. The work described here uses heterodyne spectroscopy from space to search for molecular oxygen in the interstellar medium. The Odin satellite carries a 1.1 m sub-millimeter dish and a dedicated 119 GHz receiver for the ground state line of O2. Starting in 2002, the star forming molecular cloud core rho Oph A was observed with Odin for 34 days during several observing runs. We detect a spectral line at v(LSR) = 3.5 km/s with dv(FWHM) = 1.5 km/s, parameters which are also common to other species associated with rho Ohp A. This feature is identified as the O2 (N_J = 1_1 - 1_0) transition at 118 750.343 MHz. The abundance of molecular oxygen, relative to H2,, is 5E-8 averaged over the Odin beam. This abundance is consistently lower than previously reported upper limits.

We have mapped [Ne II] (12.8um) and [S IV] (10.5um) emission from W51 IRS 2 with TEXES on Gemini North, and we compare these data to VLA free-free observations and VLT near-infrared images. With 0.5" spatial and 4 km/s spectral resolution we are able to separate the ionized gas into several components: an extended H II region on the front surface of the molecular cloud, several embedded compact H II regions, and a streamer of high velocity gas. We interpret the high velocity streamer as a precessing or fan-like jet, which has emerged from the molecular cloud into an OB star cluster where it is being ionized.

Excess emission from a point-like source coincident with the central star of the Helix Nebula is detected with Spitzer at 8, 24, and 70 um. At 24 um, the central source is superposed on an extended diffuse emission region. While the [OIV] 25.89 um line contributes to the diffuse emission, a 10-35 um spectrum of the central source shows a strong thermal continuum. The excess emission from the star most likely originates from a dust disk with blackbody temperatures of 90--130 K. Assuming a simple optically thin debris disk model, the dust is distributed in a ring between ~35 and ~150 AU from the central star, possibly arising from collisions of Kuiper-Belt-like Objects or the break-up of comets from an Oort-like cloud that have survived from the post-main-sequence evolution of the central star.

Color-magnitude diagrams (CMDs) are presented for the Spitzer SAGE (Surveying the Agents of a Galaxy's Evolution) survey of the Large Magellanic Cloud (LMC). IRAC and MIPS 24 um epoch one data are presented. These data represent the deepest, widest mid-infrared CMDs of their kind ever produced in the LMC. Combined with the 2MASS survey, the diagrams are used to delineate the evolved stellar populations in the Large Magellanic Cloud as well as Galactic foreground and extragalactic background populations. Some 32000 evolved stars brighter than the tip of the red giant branch are identified. Of these, approximately 17500 are classified as oxygen-rich, 7000 carbon-rich, and another 1200 as ``extreme'' asymptotic giant branch (AGB) stars. Brighter members of the latter group have been called ``obscured'' AGB stars in the literature owing to their dusty circumstellar envelopes. A large number (1200) of luminous oxygen--rich AGB stars/M supergiants are also identified. Finally, there is strong evidence from the 24 um MIPS channel that previously unexplored, lower luminosity oxygen-rich AGB stars contribute significantly to the mass loss budget of the LMC (1200 such sources are identified).

High-resolution mid-infrared images have been obtained in N-band and Q-band for the proto-planetary nebula IRAS 16594-4656. A bright equatorial torus and a pair of bipolar lobes can clearly be seen in the infrared images. The torus appears thinner at the center than at the edges, suggesting that it is viewed nearly edge-on. The infrared lobes correspond to the brightest lobes of the reflection nebula seen in the Hubble Space Telescope (HST) optical image, but with no sign of the point-symmetric structure seen in the visible image. The lobe structure shows a close correspondence with a molecular hydrogen map obtained with HST, suggesting that the dust emission in the lobes traces the distribution of the shocked gas. The shape of the bipolar lobes shows clearly that the fast outflow is still confined by the remnant circumstellar envelope of the progenitor asymptotic giant branch (AGB) star. However, the non-detection of the dust outside of the lobes suggests that the temperature of the dust in the AGB envelope is too low for it to be detected at 20 microns.

We are performing a uniform and unbiased, ~7x7 degrees imaging survey of the Large Magellanic Cloud (LMC), using the IRAC and MIPS instruments on board the Spitzer Space Telescope in order to survey the agents of a galaxy's evolution (SAGE), the interstellar medium (ISM) and stars in the LMC. The detection of diffuse ISM with column densities >1.2x10^21 H cm^-2 permits detailed studies of dust processes in the ISM. SAGE's point source sensitivity enables a complete census of newly formed stars with masses >3 solar masses that will determine the current star formation rate in the LMC. SAGE's detection of evolved stars with mass loss rates >1x10^-8 solar masses per year will quantify the rate at which evolved stars inject mass into the ISM of the LMC. The observing strategy includes two epochs in 2005, separated by three months, that both mitigate instrumental artifacts and constrain source variability. The SAGE data are non-proprietary. The data processing includes IRAC and MIPS pipelines and a database for mining the point source catalogs, which will be released to the community in support of Spitzer proposal cycles 4 and 5. We present initial results on the epoch 1 data with a special focus on the N79 and N83 region. The SAGE epoch 1 point source catalog has ~4 million sources. The point source counts are highest for the IRAC 3.6 microns band and decrease dramatically towards longer wavelengths consistent with the fact that stars dominate the point source catalogs and that the dusty objects, e.g. young stellar objects and dusty evolved stars that detected at the longer wavelengths, are rare in comparison. We outline a strategy for identifying foreground MW stars, that may comprise as much as 18% of the source list, and background galaxies, that may comprise ~12% of the source list.

The bright GRB 050408 was localized by HETE-II near local midnight, enabling an impressive ground-based followup effort as well as space-based followup from Swift. The Swift data from the X-Ray Telescope (XRT) and our own optical photometry and spectrum of the afterglow provide the cornerstone for our analysis. Under the traditional assumption that the visible waveband was above the peak synchrotron frequency and below the cooling frequency, the optical photometry from 0.03 to 5.03 days show an afterglow decay corresponding to an electron energy index of p_lc = 2.05 +/- 0.04, without a jet break as suggested by others. A break is seen in the X-ray data at early times (at ~12600 sec after the GRB). The spectral slope of the optical spectrum is consistent with p_lc assuming a host-galaxy extinction of A_V = 1.18 mag. The optical-NIR broadband spectrum is also consistent with p = 2.05, but prefers A_V = 0.57 mag. The X-ray afterglow shows a break at 1.26 x 10^4 sec, which may be the result of a refreshed shock. This burst stands out in that the optical and X-ray data suggest a large H I column density of N_HI ~ 10^22 cm^-2; it is very likely a damped Lyman alpha system and so the faintness of the host galaxy (M_V > -18 mag) is noteworthy. Moreover, we detect extraordinarily strong Ti II absorption lines with a column density through the GRB host that exceeds the largest values observed for the Milky Way by an order of magnitude. Furthermore, the Ti II equivalent width is in the top 1% of Mg II absorption-selected QSOs. This suggests that the large-scale environment of GRB 050408 has significantly lower Ti depletion than the Milky Way and a large velocity width (delta v > 200 km/s).

We present Odin observations of the AGB star W Hya in the ground-state transition of ortho-H2O, 1(10)-1(01), at 557GHz. The line is clearly of circumstellar origin. Radiative transfer modelling of the water lines observed by Odin and ISO results in a mass-loss rate of (2.5 +/- 0.5)E-7 Msol/yr, and a circumstellar H2O abundance of (2.0 +/- 1.0)E-3. The inferred mass-loss rate is consistent with that obtained from modelling the circumstellar CO radio line emission, and also with that obtained from the dust emission modelling combined with a dynamical model for the outflow. The very high water abundance, higher than the cosmic oxygen abundance, can be explained by invoking an injection of excess water from evaporating icy bodies in the system. The required extra mass of water is quite small, on the order of ~ 0.1 earth mass.

Ever since the pioneering study of Spitzer, it has been widely recognized that grains play an important role in the heating and cooling of photo-ionized environments. A detailed code is necessary to model grains in a photo-ionized medium. In this paper we will use the spectral synthesis code Cloudy for this purpose. The grain model of this code has recently undergone a comprehensive upgrade. Part of this is the newly developed hybrid grain charge model, which will be described in detail. This model allows discrete charge states of very small grains to be modelled accurately while simultaneously avoiding the overhead of fully resolving the charge distribution of large grains, thus making the model both accurate and computationally efficient. A comprehensive comparison with the fully resolved charge state models of Weingartner & Draine (2001) shows that the agreement is very satisfactory for realistic size distributions. The effect of the grain size distribution on the line emission from photo-ionized regions is studied by taking standard models for an H II region and a planetary nebula and adding a dust component to the models with varying grain size distributions. A comparison of the models shows that varying the size distribution has a dramatic effect on the emitted spectrum, and affects the ionization balance as well. These results clearly demonstrate that the grain size distribution is an important parameter in photo-ionization models.

An Artificial Neural Network (ANN) has been employed using a supervised back-propagation scheme to classify 2000 bright sources from the Calgary database of IRAS (Infrared Astronomy Satellite) spectra in the wavelength region of 8-23 microns. The data base has been classified into 17 pre-determined classes based on spectral morphology. We have been able to classify more than 80 percent of the 2000 sources correctly at the first instance. The speed and robustness of the scheme will allow us to classify the whole of LRS database, containing more than 50,000 sources in the future.

In this paper we present results obtained with the new grain code in Cloudy which underline the strong effect of photo-electric heating by grains in photo-ionized regions. We study the effect that the distribution of grain sizes has on the magnitude of the effect, and show that this effect is nothing short of dramatic. This makes the grain size distribution an important parameter in modeling of photo-ionized regions such as H II regions and planetary nebulae.

We have mapped the submillimeter emission from the 1(10)-1(01) transition of ortho-water in the W3 star-forming region. A 5'x5' map of the W3 IRS4 and W3 IRS5 region reveals strong water lines at half the positions in the map. The relative strength of the Odin lines compared to previous observations by SWAS suggests that we are seeing water emission from an extended region. Across much of the map the lines are double-peaked, with an absorption feature at -39 km/s; however, some positions in the map show a single strong line at -43 km/s. We interpret the double-peaked lines as arising from optically thick, self-absorbed water emission near the W3 IRS5, while the narrower blue-shifted lines originate in emission near W3 IRS4. In this model, the unusual appearance of the spectral lines across the map results from a coincidental agreement in velocity between the emission near W3 IRS4 and the blue peak of the more complex lines near W3 IRS5. The strength of the water lines near W3 IRS4 suggests we may be seeing water emission enhanced in a photon-dominated region.

Odin has successfully observed the molecular core rho Oph A in the 572.5 GHz rotational ground state line of ammonia, NH3 (J,K = 1,0 - 0,0). The interpretation of this result makes use of complementary molecular line data obtained from the ground (C17O and CH3OH) as part of the Odin preparatory work. Comparison of these observations with theoretical model calculations of line excitation and transfer yields a quite ordinary abundance of methanol, X(CH3OH) = 3e-9. Unless NH3 is not entirely segregated from C17O and CH3OH, ammonia is found to be significantly underabundant with respect to typical dense core values, viz. X(NH3) = 8e-10.

Odin has successfully observed three regions in the Orion A cloud, i.e. Ori KL, Ori S and the Orion Bar, in the 572.5 GHz rotational ground state line of ammonia, ortho-NH3 (J,K) = (1,0) -> (0,0), and the result for the Orion Bar represents the first detection in an ammonia line. Several velocity components are present in the data. Specifically, the observed line profile from the Orion Bar can be decomposed into two components, which are in agreement with observations in high-J CO lines by Wilson et al. 2001. Using the source model for the Orion Bar by these authors, our Odin observation implies a total ammonia abundance of NH3/H2 = 5E-9.

Ever since the pioneering study of Spitzer, it has been widely recognized that grains play an important role in the heating and cooling of photo-ionized environments. This includes the diffuse ISM, as well as H II regions, planetary nebulae, and photo-dissociation regions. A detailed code is necessary to model grains in a photo-ionized medium since the interactions of grains with their environment include a host of microphysical processes, and their importance can only be judged by performing a complete simulation. In this paper we will use the spectral synthesis code Cloudy for this purpose. A comprehensive upgrade of the grain model has been recently incorporated in Cloudy, and certain aspects of this upgrade will be discussed. Special emphasis will be on the new grain charge model. We will consider in detail the physics of grains in both ionized and neutral environments, and will present a calculation of photo-electric heating and collisional cooling rates for a range of physical conditions and grain materials and for a range of grain sizes (including a realistic size distribution). We conclude with a brief discussion of the problems currently hampering progress in this field. The new grain model will be used to model the silicate emission in the Ney-Allen nebula, and will help us better understand the nature of the grains in that part of the Orion complex.

We have fitted the HST images and the spectral energy distributions (SEDs) of three proto-planetary nebulae (PPN) with a 2-D radiation transfer model. The geometric and mass loss properties of these PPN are also derived.

The optical spectrum of the infrared source IRAS 04296+3429 (optical counterpart-G0 Ia star, V=14.2) was obtained with the echelle spectrometer PFES at the prime focus of the 6 m telescope. We discover emission bands (0,0) and (0,1) of the Swan system of the C2 molecule in the optical spectrum of IRAS 04296+3429. Comparison with the spectrum of the Hale-Bopp comet leads us to propose that in both cases the same mechanism (resonance fluorescence) is responsible for the emission in the C2 molecular bands. Several strong absorption features whose positions coincide with known diffuse interstellar bands are revealed in the spectrum of IRAS 04296+3429. The infrared spectrum of IRAS 04296+3429 shows the famous 21 um feature, but this object has not been observed by KAO. However, like IRAS 05113+1347, IRAS 05341+0852 and IRAS 22223+4327, our detailed modelling of its spectral energy distribution suggested that this source also should show the 30 um band. In fact, ISO discovered a broad, relatively strong feature around 30 um for IRAS 04296+3429. The surface chemical composition of the source IRAS 04296+3429 is metal-deficient (the averaged value of the abundances of the iron group elements Ti, V, Cr and Fe relative to the solar values is [M/H]=-0.9 and has been considerably altered during the evolution: carbon, nitrogen and s-process elements are overabundant relative to the metallicity. The totality of physical and chemical parameters derived for IRAS 04296+3429 confirms a relation between presence of the feature at 21 um in the spectrum of a carbon rich star and an excess of the s-process elements.