Search

We report a statistically significant detection of nongravitational acceleration on the sub-kilometer near-Earth asteroid (523599) 2003 RM. Due to its orbit, 2003 RM experiences favorable observing apparitions every 5 years. Thus, since its discovery, 2003 RM has been extensively tracked with ground-based optical facilities in 2003, 2008, 2013, and 2018. We find that the observed plane-of-sky positions cannot be explained with a purely gravity-driven trajectory. Including a transverse nongravitational acceleration allows us to match all observational data, but its magnitude is inconsistent with perturbations typical of asteroids such as the Yarkovsky effect or solar radiation pressure. After ruling out that the orbital deviations are due to a close approach or collision with another asteroid, we hypothesize that this anomalous acceleration is caused by unseen cometary outgassing. A detailed search for evidence of cometary activity with archival and deep observations from Pan-STARRS and the VLT does not reveal any detectable dust production. However, the best-fitting H$_2$O sublimation model allows for brightening due to activity consistent with the scatter of the data. We estimate the production rate required for H$_2$O outgassing to power the acceleration, and find that, assuming a diameter of 300 m, 2003 RM would require Q(H$_2$O)$\sim10^{23}$ molec s$^{-1}$ at perihelion. We investigate the recent dynamical history of 2003 RM and find that the object most likely originated in the mid-to-outer main belt ($\sim86\%$) as opposed to from the Jupiter-family comet region ($\sim11\%$). Further observations, especially in the infrared, could shed light on the nature of this anomalous acceleration., Comment: 14 pages, 9 figures, 2 tables, accepted by the Planetary Science Journal

2009 MS9 is a trans-Neptunian object (TNO) whose perihelion brings it close to the distance where some long period comets are seen to become active. Knowing this, and the fact that this object appears to brighten in excess of it's predicted nucleus brightness suggests that 2009 MS9 has a delayed onset of activity brought on by the sublimation of a species more volatile than water. In this paper we characterize 2009 MS9's physical properties and investigate potential outgassing through composite images, sublimation models, and measurements of spectral reflectivity. We find that deep composite images of the object at various epochs along its orbit show no evidence of dust yet place sensitive limits to the dust production. We estimate the nucleus radius to be 11.5 km $\pm 3.5$ km using thermal IR modeling from NEOWISE data and use this and data pre-perihelion to estimate a geometric albedo of 0.25. We compare a CO-sublimation activity model to its post perihelion heliocentric light curve and find this data supports an active fractional area of $5 \times 10^{-6}$ assuming 2 $\mu$m sized grains and other typical comet parameters. The spectral reflectivity of the surface materials obtained with the Gemini Observatory and CFHT at different epochs shows a reddening spectral slope. We compare the physical properties of 2009 MS9 to both TNO and comet populations, and speculate that 2009 MS9's reddening may be due to the buildup of a dust mantle on the surface and could be an explanation of why TNOs exhibit a color bimodality., Comment: 25 pages, 8 figures, Accepted to PSJ

Manx objects approach the inner solar system on long-period comet (LPC) orbits with the consequent high inbound velocities, but unlike comets, Manxes display very little to no activity even near perihelion. This suggests that they may have formed in circumstances different from typical LPCs; moreover, this lack of significant activity also renders them difficult to detect at large distances. Thus, analyzing their physical properties can help constrain models of solar system formation as well as sharpen detection methods for those classified as NEOs. Here, we focus on the Manx candidate A/2018 V3 as part of a larger effort to characterize Manxes as a whole. This particular object was observed to be inactive even at its perihelion at $q$ = 1.34 au in 2019 September. Its spectral reflectivity is consistent with typical organic-rich comet surfaces with colors of $g'-r'= 0.67\pm0.02$, $r'-i' = 0.26\pm0.02$, and $r'-z' = 0.45\pm0.02$, corresponding to a spectral reflectivity slope of $10.6\pm 0.9$ %/100nm. A least-squares fit of our constructed light curve to the observational data yields an average nucleus radius of $\approx$2 km assuming an albedo of 0.04. This is consistent with the value measured from NEOWISE. A surface brightness analysis for data taken 2020 July 13 indicated possible low activity ($\lesssim0.68$ g $\rm s^{-1}$), but not enough to lift optically significant amounts of dust. Finally, we discuss Manxes as a constraint on solar system dynamical models as well as their implications for planetary defense., Comment: 13 pages, 7 figures; updated to add NEOWISE acknowledgment

Comet 49P/Arend-Rigaux is a well known low-activity Jupiter Family comet. Despite the low activity, we have witnessed outgassing activity in 1992, 2004, and 2012. In 2012 a broad tail-like feature (PA$\sim270^\circ, \sim2.3\times10^5$ km) and a narrow jet-like feature (PA$\sim180^\circ, \sim9.3\times10^4$ km) were seen simultaneously. Using Finson-Probstein (FP) dust dynamical models we determine: grain sizes released in each event; duration of activity; when activity peaked; and velocity of the dust particles, allowing us to make comparisons between the events. We find that the tail feature in 2012 is similar to the tail in 1992 with large grains (40-4000 $\mu$m) peaking in activity near perihelion with a long outgassing duration greater than 150 days. The jet feature from 2012, however, is more similar to the 2004 event which we model with small grains (1-8 $\mu$m) with a short duration of activity ($\sim$1 month). The main difference between these two features is that the 2004 event occurs prior to perihelion, while the 2012 event is post-perihelion. We use the grain sizes from the FP models to constrain ice sublimation models. Between 1985 and 2018 we cover 6 apparitions with 26 nights of our own observations plus data from the literature and the Minor Planet Center, which together, allow us to model the heliocentric light curve. We find that the models are consistent with H$_2$O ice sublimation as the volatile responsible for driving activity over most of the active phases and a combination of H$_2$O and CO$_2$ ices are responsible for driving activity near perihelion. We measure the fractional active area over time for H$_2$O and discover that the activity decreases from an average active area of $\sim3\%$ to $\sim0.2\%$. This secular decrease in activity implies that the comet is becoming depleted of volatiles and is in the process of transitioning to a dormant or dead state., Comment: 18 pages, 8 figures, Accepted for publication by Icarus

The detection of Interstellar Objects passing through the Solar System offers the promise of constraining the physical and chemical processes involved in planetary formation in other extrasolar systems. While the effect of outgassing by 1I/2017 U1 ('Oumuamua) was dynamically observed, no direct detection of the ejected material was made. The discovery of the active interstellar comet 2I/Borisov means spectroscopic investigations of the sublimated ices is possible for this object. We report the first detection of gas emitted by an interstellar comet via the near-UV emission of CN from 2I/Borisov at a heliocentric distance of $r$ = 2.7 au on 2019 September 20. The production rate was found to be Q(CN) = $(3.7\pm0.4)\times10^{24}$ s$^{-1}$, using a simple Haser model with an outflow velocity of 0.5 km s$^{-1}$. No other emission was detected, with an upper limit to the production rate of C$_2$ of $4\times10^{24}$ s$^{-1}$. The spectral reflectance slope of the dust coma over $3900$ \AA\ $< \lambda< 6000$ \AA \ is steeper than at longer wavelengths, as found for other comets. Broad band $R_c$ photometry on 2019 September 19 gave a dust production rate of $Af\rho=143\pm10$ cm. Modelling of the observed gas and dust production rates constrains the nuclear radius to $0.7-3.3$ km assuming reasonable nuclear properties. Overall, we find the gas, dust and nuclear properties for the first active Interstellar Object are similar to normal Solar System comets., Comment: 9 pages, 4 figures, accepted for publication in ApJL

We report on the catastrophic disintegration of P/2016 G1 (PANSTARRS), an active asteroid, in April 2016. Deep images over three months show that object is constituted by a central concentration of fragments surrounded by an elongated coma, and presents previously unreported sharp arc-like and a narrow linear features. The morphology and evolution of these characteristics independently point toward a brief event on 2016 March 6. The arc and the linear feature can be reproduced by large particles on a ring, moving at ~2.5 m/s. The expansion of the ring defines a cone with a ~40deg half-opening. We propose that G1 was hit by a small object which caused its (partial or total) disruption, and that the ring corresponds to large fragments ejected during the final stages of the crater formation.}

On 2019 January 5 a streamer associated with the 4--10 km main-belt asteroid (6478)~Gault was detected by the ATLAS sky survey, a rare discovery of activity around a main-belt asteroid. Archival data from ATLAS and Pan-STARRS1 show the trail in early December 2018, but not between 2010 and January 2018. The feature has significantly changed over one month, perfectly matching predictions of pure dust dynamical evolution and changes in observing geometry for a short release of dust around 2018 October 28. Follow-up observations with HST show a second narrow trail corresponding to a brief release of dust on 2018 December 30. Both releases occurred with negligible velocity. We find the dust grains to be fairly large, with power-law size distributions in the $10^{-5} - 10^{-3}$~m range and power-law indices of $\sim -1.5$. Three runs of ground-based data find a signature of $\sim 2\,\rm h$ rotation, close to the rotational limit, suggesting that the activity is the result of landslides or reconfigurations after YORP spin-up., Comment: 9 pages, 4 figures, 1 table

We show that `Oumuamua's excited spin could be in a high energy LAM state, which implies that its shape could be far from the highly elongated shape found in previous studies. CLEAN and ANOVA algorithms are used to analyze `Oumuamua's lightcurve using 818 observations over 29.3~days. Two fundamental periodicities are found at frequencies (2.77$\pm$0.11) and (6.42$\pm$0.18)~cycles/day, corresponding to (8.67$\pm$0.34)~h and (3.74$\pm$0.11)~h, respectively. The phased data show that the lightcurve does not repeat in a simple manner, but approximately shows a double minimum at 2.77~cycles/day and a single minimum at 6.42~cycles/day. This is characteristic of an excited spin state. `Oumuamua could be spinning in either the long (LAM) or short (SAM) axis mode. For both, the long axis precesses around the total angular momentum vector with an average period of (8.67$\pm$0.34)~h. For the three LAMs we have found, the possible rotation periods around the long axis are 6.58, 13.15, or 54.48~h, with 54.48~h being the most likely. `Oumuamua may also be nutating with respective periods of half of these values. We have also found two possible SAM states where `Oumuamua oscillates around the long axis with possible periods at 13.15 and 54.48~h, the latter as the most likely. In this case any nutation will occur with the same periods. Determination of the spin state, the amplitude of the nutation, the direction of the TAMV, and the average total spin period may be possible with a direct model fit to the lightcurve. We find that `Oumuamua is "cigar-shaped"', if close to its lowest rotational energy, and an extremely oblate spheroid if close to its highest energy state for its total angular momentum., Comment: 8 pages, 3 figures

Comet C/2017 K2 (PANSTARRS) was discovered by the Pan-STARRS1 (PS1) Survey on 2017 May 21 at a distance 16.09 au from the Sun, the second most distant discovery of an active comet. Pre-discovery images in the PS1 archive back to 2014 and additional deep CFHT images between 2013 May 10-13 showed the comet to be active at 23.75 au. We derive an upper limit to the nucleus radius of $R_N$=80 km, assuming a 4\% albedo. The spectral reflectivity of the comet surface is similar to "fresh" regions seen on comet 67P/Churyumov-Gerasimenko using the $Rosetta$ OSIRIS camera. Pre-discovery photometry combined with new data obtained with Megacam on the CFHT show that the activity is consistent with CO-ice sublimation and inconsistent with CO$_2$-ice sublimation. The ice sublimation models were run out to perihelion in 2022 at 1.8 au to predict the CO production rates, assuming that the outgassing area does not change. Assuming a canonical 4\% active surface area for water-ice sublimation, we present production rate ratios, $Q_{\rm CO}$/$Q_{\rm H2O}$, for a range of nucleus sizes. Comparing these results with other CO-rich comets we derive a lower limit to the nucleus radius of $\sim$14 km. We present predictions for $Q_{\rm CO}$ at a range of distances that will be useful for planning observations with JWST and large ground-based facilities., Comment: 5 pages, 4 figures

We present initial time-resolved observations of the split comet 332P/Ikeya-Murakami taken using the Hubble Space Telescope. Our images reveal a dust-bathed cluster of fragments receding from their parent nucleus at projected speeds in the range 0.06 to 3.5 m s$^{-1}$ from which we estimate ejection times from October to December 2015. The number of fragments with effective radii $\gtrsim$20 m follows a differential power law with index $\gamma$ = -3.6$\pm$0.6, while smaller fragments are less abundant than expected from an extrapolation of this power-law. We argue that, in addition to losses due to observational selection, torques from anisotropic outgassing are capable of destroying the small fragments by driving them quickly to rotational instability. Specifically, the spin-up times of fragments $\lesssim$20 m in radius are shorter than the time elapsed since ejection from the parent nucleus. The effective radius of the parent nucleus is $r_e \le$ 275 m (geometric albedo 0.04 assumed). This is about seven times smaller than previous estimates and results in a nucleus mass at least 300 times smaller than previously thought. The mass in solid pieces, $2\times10^9$ kg, is about 4% of the mass of the parent nucleus. As a result of its small size, the parent nucleus also has a short spin-up time. Brightness variations in time-resolved nucleus photometry are consistent with rotational instability playing a role in the release of fragments., Comment: 19 pages, 1 table, 4 figures, To be published on ApJL

We describe 2016 January to April observations of the fragments of 332P/Ikeya-Murakami, a comet earlier observed in a 2010 October outburst (Ishiguro et al 2014). We present photometry of the fragments, and perform simulations to infer the time of breakup. We argue that the eastern-most rapidly brightening fragment ($F4$) best corresponds to the original nucleus, rather than the initial bright fragment $F1$. We compute radial and tangential non-gravitational parameters, $A_1 = (1.5 \pm 0.4) \times 10^{-8}$ AU day$^{-2}$ and $(7.2 \pm 1.9) \times 10^{-9}$ AU day$^{-2}$; both are consistent with zero at the $4\sigma$ level. Monte Carlo simulations indicate that the fragments were emitted on the outbound journey well after the 2010 outburst, with bright fragment $F1$ splitting in mid--2013 and the fainter fragments within months of the 2016 January recovery. Western fragment $F7$ is the oldest, dating from 2011. We suggest that the delayed onset of the splitting is consistent with a self-propagating crystallization of water ice., Comment: 8 pages, 4 figures, published in ApJ Letters

We present detailed chemical element abundance ratios of 17 elements in three metal poor stars in the Ursa Minor dwarf spheroidal galaxy, which we combine with extant data from the literature to assess the predictions of a novel suite of galaxy chemical evolution models. The spectroscopic data were obtained with the Keck/HIRES instrument and revealed low metallicities of [Fe/H]=-2.12, -2.13 and -2.67 dex. While the most metal poor star in our sample shows an overabundance of [Mn/Fe] and other Fe-peak elements, our overall findings are in agreement with previous studies of this galaxy: elevated values of the [alpha/Fe] ratios that are similar to, or only slightly lower than, the halo values but with SN Ia enrichment at very low metallicity, as well as an enhancement of the ratio of first to second peak neutron capture elements [Y/Ba] with decreasing metallicity. The chemical evolution models which were tailored to reproduce the metallicity distribution function of the dSph, indicate that UMi had an extended star formation which lasted nearly 5 Gyr with low efficiency and are able to explain the [Y/Ba] enhancement at low metallicity for the first time. In particular, we show that the present day lack of gas is probably due to continuous loss of gas from the system, which we model as winds., Comment: 10 pages, 7 figures, tables

We present an observational and dynamical study of newly discovered main-belt comet 313P/Gibbs. We find that the object is clearly active both in observations obtained in 2014 and in precovery observations obtained in 2003 by the Sloan Digital Sky Survey, strongly suggesting that its activity is sublimation-driven. This conclusion is supported by a photometric analysis showing an increase in the total brightness of the comet over the 2014 observing period, and dust modeling results showing that the dust emission persists over at least three months during both active periods, where we find start dates for emission no later than 2003 July 24+/-10 for the 2003 active period and 2014 July 28+/-10 for the 2014 active period. From serendipitous observations by the Subaru Telescope in 2004 when the object was apparently inactive, we estimate that the nucleus has an absolute R-band magnitude of H_R=17.1+/-0.3, corresponding to an effective nucleus radius of r_e~1.00+/-0.15 km. The object's faintness at that time means we cannot rule out the presence of activity, and so this computed radius should be considered an upper limit. We find that 313P's orbit is intrinsically chaotic, having a Lyapunov time of T_l=12000 yr and being located near two 3-body mean-motion resonances with Jupiter and Saturn, 11J-1S-5A and 10J+12S-7A, yet appears stable over >50 Myr in an apparent example of stable chaos. We furthermore find that 313P is the second main-belt comet, after P/2012 T1 (PANSTARRS), to belong to the ~155 Myr old Lixiaohua asteroid family., Comment: 8 pages, 4 figures; accepted for publication in ApJ Letters

We analyze 760475 observations of 333026 main-belt objects obtained by the Pan-STARRS1 (PS1) survey telescope between 2012 May 20 and 2013 November 9, a period during which PS1 discovered two main-belt comets, P/2012 T1 (PANSTARRS) and P/2013 R3 (Catalina-PANSTARRS). PS1 comet detection procedures currently consist of the comparison of the point spread functions (PSFs) of moving objects to those of reference stars, and the flagging of objects that show anomalously large radial PSF widths. Based on the number of missed discovery opportunities among comets discovered by other observers, we estimate an upper limit comet discovery efficiency rate of ~70% for PS1. Additional analyses that could improve comet discovery yields in future surveys include linear PSF analysis, modeling of trailed stellar PSFs for comparison to trailed moving object PSFs, searches for azimuthally localized activity, comparison of point-source-optimized photometry to extended-source-optimized photometry, searches for photometric excesses in objects with known absolute magnitudes, and crowd-sourcing. Analysis of PS1 survey statistics indicates an expected fraction of 59 MBCs per 10^6 outer main-belt asteroids, and a 95% confidence upper limit of 96 MBCs per 10^6 outer main-belt asteroids. We note that more sensitive future surveys could detect many more MBCs than estimated here. We find an excess of high eccentricities (0.1 < e < 0.3) among all known MBCs relative to the background asteroid population. Theoretical calculations show that, given these eccentricities, the sublimation rate for a typical MBC is orders of magnitude larger at perihelion than at aphelion, providing a plausible physical explanation for the observed behavior of MBCs peaking in observed activity strength near perihelion., Comment: 38 pages, 14 figures; accepted for publication in Icarus

We have calculated 90% confidence limits on the steady-state rate of catastrophic disruptions of main belt asteroids in terms of the absolute magnitude at which one catastrophic disruption occurs per year (HCL) as a function of the post-disruption increase in brightness (delta m) and subsequent brightness decay rate (tau). The confidence limits were calculated using the brightest unknown main belt asteroid (V = 18.5) detected with the Pan-STARRS1 (Pan-STARRS1) telescope. We measured the Pan-STARRS1's catastrophic disruption detection efficiency over a 453-day interval using the Pan-STARRS moving object processing system (MOPS) and a simple model for the catastrophic disruption event's photometric behavior in a small aperture centered on the catastrophic disruption event. Our simplistic catastrophic disruption model suggests that delta m = 20 mag and 0.01 mag d-1 < tau < 0.1 mag d-1 which would imply that H0 = 28 -- strongly inconsistent with H0,B2005 = 23.26 +/- 0.02 predicted by Bottke et al. (2005) using purely collisional models. We postulate that the solution to the discrepancy is that > 99% of main belt catastrophic disruptions in the size range to which this study was sensitive (100 m) are not impact-generated, but are instead due to fainter rotational breakups, of which the recent discoveries of disrupted asteroids P/2013 P5 and P/2013 R3 are probable examples. We estimate that current and upcoming asteroid surveys may discover up to 10 catastrophic disruptions/year brighter than V = 18.5., Comment: 61 Pages, 10 Figures, 3 Tables

We present the results of a search for the reactivation of active asteroid 176P/LINEAR during its 2011 perihelion passage using deep optical observations obtained before, during, and after that perihelion passage. Deep composite images of 176P constructed from data obtained between June 2011 and December 2011 show no visible signs of activity, while photometric measurements of the object during this period also show no significant brightness enhancements similar to that observed for 176P between November 2005 and December 2005 when it was previously observed to be active. An azimuthal search for dust emission likewise reveals no evidence for directed emission (i.e., a tail, as was previously observed for 176P), while a one-dimensional surface brightness profile analysis shows no indication of a spherically symmetric coma at any time in 2011. We conclude that 176P did not in fact exhibit activity in 2011, at least not on the level on which it exhibited activity in 2005, and suggest that this could be due to the devolatization or mantling of the active site responsible for its activity in 2005., Comment: 17 pages, 9 figures; published in AJ

We report photometric observations for comet C/2012 S1 (ISON) obtained during the time period immediately after discovery (r=6.28 AU) until it moved into solar conjunction in mid-2013 June using the UH2.2m, and Gemini North 8-m telescopes on Mauna Kea, the Lowell 1.8m in Flagstaff, the Calar Alto 1.2m telescope in Spain, the VYSOS-5 telescopes on Mauna Loa Hawaii and data from the CARA network. Additional pre-discovery data from the Pan STARRS1 survey extends the light curve back to 2011 September 30 (r=9.4 AU). The images showed a similar tail morphology due to small micron sized particles throughout 2013. Observations at sub-mm wavelengths using the JCMT on 15 nights between 2013 March 9 (r=4.52 AU) and June 16 (r=3.35 AU) were used to search for CO and HCN rotation lines. No gas was detected, with upper limits for CO ranging between (3.5-4.5)E27 molec/s. Combined with published water production rate estimates we have generated ice sublimation models consistent with the photometric light curve. The inbound light curve is likely controlled by sublimation of CO2. At these distances water is not a strong contributor to the outgassing. We also infer that there was a long slow outburst of activity beginning in late 2011 peaking in mid-2013 January (r~5 AU) at which point the activity decreased again through 2013 June. We suggest that this outburst was driven by CO injecting large water ice grains into the coma. Observations as the comet came out of solar conjunction seem to confirm our models., Comment: 8 pages, 2 figures, 3 tables

We present initial results from observations and numerical analyses aimed at characterizing main-belt comet P/2012 T1 (PANSTARRS). Optical monitoring observations were made between October 2012 and February 2013 using the University of Hawaii 2.2 m telescope, the Keck I telescope, the Baade and Clay Magellan telescopes, Faulkes Telescope South, the Perkins Telescope at Lowell Observatory, and the Southern Astrophysical Research (SOAR) telescope. The object's intrinsic brightness approximately doubles from the time of its discovery in early October until mid-November and then decreases by ~60% between late December and early February, similar to photometric behavior exhibited by several other main-belt comets and unlike that exhibited by disrupted asteroid (596) Scheila. We also used Keck to conduct spectroscopic searches for CN emission as well as absorption at 0.7 microns that could indicate the presence of hydrated minerals, finding an upper limit CN production rate of QCN<1.5x10^23 mol/s, from which we infer a water production rate of QH2O<5x10^25 mol/s, and no evidence of the presence of hydrated minerals. Numerical simulations indicate that P/2012 T1 is largely dynamically stable for >100 Myr and is unlikely to be a recently implanted interloper from the outer solar system, while a search for potential asteroid family associations reveal that it is dynamically linked to the ~155 Myr-old Lixiaohua asteroid family., Comment: 15 pages, 4 figures, accepted for publication in ApJ Letters

We describe the Pan-STARRS Moving Object Processing System (MOPS), a modern software package that produces automatic asteroid discoveries and identifications from catalogs of transient detections from next-generation astronomical survey telescopes. MOPS achieves > 99.5% efficiency in producing orbits from a synthetic but realistic population of asteroids whose measurements were simulated for a Pan-STARRS4-class telescope. Additionally, using a non-physical grid population, we demonstrate that MOPS can detect populations of currently unknown objects such as interstellar asteroids. MOPS has been adapted successfully to the prototype Pan-STARRS1 telescope despite differences in expected false detection rates, fill-factor loss and relatively sparse observing cadence compared to a hypothetical Pan-STARRS4 telescope and survey. MOPS remains >99.5% efficient at detecting objects on a single night but drops to 80% efficiency at producing orbits for objects detected on multiple nights. This loss is primarily due to configurable MOPS processing limits that are not yet tuned for the Pan-STARRS1 mission. The core MOPS software package is the product of more than 15 person-years of software development and incorporates countless additional years of effort in third-party software to perform lower-level functions such as spatial searching or orbit determination. We describe the high-level design of MOPS and essential subcomponents, the suitability of MOPS for other survey programs, and suggest a road map for future MOPS development., Comment: 57 Pages, 26 Figures, 13 Tables

We describe the recovery of faint Main Belt comet P/2008 R1 Garradd using several telescopes, culminating in a successful low $S/N$ recovery with the Gemini North telescope with GMOS. This recovery was a time-critical effort for a mission proposal, and had to be performed in a crowded field. We describe techniques and software tools for eliminating systematic noise artifacts and stellar residuals, bringing the final detection image statistics close to the Gaussian ideal for a median image stack, and achieving a detection sensitivity close to this theoretical optimum. The magnitude of $R_c$=26.1$\pm$0.2 with an assumed geometric albedo of 0.05 corresponds to a radius of 0.3 km. For ice to have survived in this object over the age of the solar system, it implies that the object is a more recent collisional fragment. We discuss the implications of the unexpectedly faint magnitude and nuclear size of P/2008 R1 on the survival of ice inside very small bodies., Comment: 15 pages, 4 figures; accepted for publications in PASP

P/2010 A2 is an object on an asteroidal orbit that was observed to have an extended tail or debris trail in January 2010. In this work, we fit the outburst of P/2010 A2 with a conical burst model, and verify previous suspicions that this was a one--time collisional event rather than an sustained cometary outburst, implying that P/2010 A2 is not a new Main Belt Comet driven by ice sublimation. We find that the best--fit cone opening angle is about 40 to 50 degrees, in agreement with numerical and laboratory simulations of cratering events. Mapping debris orbits to sky positions suggests that the distinctive arc features in the debris correspond to the same debris cone inferred from the extended dust. From the velocity of the debris, and from the presence of a velocity maximum at around 15 cm/s, we infer that the surface of A2 probably has a very low strength (<1 kPa), comparable to lunar regolith., Comment: 14 pages, 25 figures; accepted by Astronomy and Astrophysics

Main belt asteroid (300163) 2006 VW139 (later designated P/2006 VW139) was discovered to exhibit comet-like activity by the Pan-STARRS1 survey telescope using automated point-spread-function analyses performed by PS1's Moving Object Processing System. Deep follow-up observations show both a short (\sim 10") antisolar dust tail and a longer (\sim 60") dust trail aligned with the object's orbit plane, similar to the morphology observed for another main-belt comet, P/2010 R2 (La Sagra), and other well-established comets, implying the action of a long-lived, sublimation-driven emission event. Photometry showing the brightness of the near-nucleus coma remaining constant over \sim 30 days provides further evidence for this object's cometary nature, suggesting it is in fact a main-belt comet, and not a disrupted asteroid. A spectroscopic search for CN emission was unsuccessful, though we find an upper limit CN production rate of Q_CN < 1.3x10^24 mol/s, from which we infer a water production rate of Q_H2O < 10^26 mol/s. We also find an approximately linear optical spectral slope of 7.2%/1000A, similar to other cometary dust comae. Numerical simulations indicate that P/2006 VW139 is dynamically stable for > 100 Myr, while a search for a potential asteroid family around the object reveals a cluster of 24 asteroids within a cutoff distance of 68 m/s. At 70 m/s, this cluster merges with the Themis family, suggesting that it could be similar to the Beagle family to which another main-belt comet, 133P/Elst-Pizarro, belongs., Comment: 10 pages, 4 figures, accepted for publication in ApJ Letters

We present observations of comet-like main-belt object P/2010 R2 (La Sagra) obtained by Pan-STARRS 1 and the Faulkes Telescope-North on Haleakala in Hawaii, the University of Hawaii 2.2 m, Gemini-North, and Keck I telescopes on Mauna Kea, the Danish 1.54 m telescope at La Silla, and the Isaac Newton Telescope on La Palma. An antisolar dust tail is observed from August 2010 through February 2011, while a dust trail aligned with the object's orbit plane is also observed from December 2010 through August 2011. Assuming typical phase darkening behavior, P/La Sagra is seen to increase in brightness by >1 mag between August 2010 and December 2010, suggesting that dust production is ongoing over this period. These results strongly suggest that the observed activity is cometary in nature (i.e., driven by the sublimation of volatile material), and that P/La Sagra is therefore the most recent main-belt comet to be discovered. We find an approximate absolute magnitude for the nucleus of H_R=17.9+/-0.2 mag, corresponding to a nucleus radius of ~0.7 km, assuming an albedo of p=0.05. Using optical spectroscopy, we find no evidence of sublimation products (i.e., gas emission), finding an upper limit CN production rate of Q_CN<6x10^23 mol/s, from which we infer an H2O production rate of Q_H2O<10^26 mol/s. Numerical simulations indicate that P/La Sagra is dynamically stable for >100 Myr, suggesting that it is likely native to its current location and that its composition is likely representative of other objects in the same region of the main belt, though the relatively close proximity of the 13:6 mean-motion resonance with Jupiter and the (3,-2,-1) three-body mean-motion resonance with Jupiter and Saturn mean that dynamical instability on larger timescales cannot be ruled out., Comment: 23 pages, 13 figures, accepted for publication in AJ

We present results from a wide-field imaging campaign at the Canada-France-Hawaii Telescope to study the spectacular outburst of comet 17P/Holmes in late 2007. Using image-processing techniques we probe inside the spherical dust coma and find sixteen fragments having both spatial distribution and kinematics consistent with isotropic ejection from the nucleus. Photometry of the fragments is inconsistent with scattering from monolithic, inert bodies. Instead, each detected fragment appears to be an active cometesimal producing its own dust coma. By scaling from the coma of the primary nucleus of 17P/Holmes, assumed to be 1.7 km in radius, we infer that the sixteen fragments have maximum effective radii between ~ 10 m and ~ 100 m on UT 2007 Nov. 6. The fragments subsequently fade at a common rate of ~ 0.2 mag/day, consistent with steady depletion of ices from these bodies in the heat of the Sun. Our characterization of the fragments supports the hypothesis that a large piece of material broke away from the nucleus and crumbled, expelling smaller, icy shards into the larger dust coma around the nucleus., Comment: 41 pages, 12 figures. Accepted for publication by the Astronomical Journal.

There has been a vast recent improvement in photometric and kinematic data for star clusters, Ultra Compact dwarfs, galactic nuclei, and local dSph galaxies, with Subaru contributing substantially to the photometric studies in particular. These data show that there is a bimodal distribution in half-light radii, with stable star clusters always being smaller than 35pc, while stable galaxies are always larger than 120pc. We extend the previously known observational relationships and interpret them in terms of a more fundamental pair of intrinsic properties of dark matter itself: dark matter forms cored mass distributions, with a core scale length of greater than about 100pc, and always has a maximum central mass density with a narrow range. The dark matter in dSph galaxies appears to be clustered such that there is a mean volume mass density within the stellar distribution which has the very low value of about 0.1$M_{\odot}$ pc$^{-3}$. None of the dSphs displays kinematics which require the presence of an inner cusp, while in two dSphs there is evidence that the density profile is shallow (cored) in the inner regions. The maximum central dark matter density derived is model dependent, but is likely to have a mean value (averaged over a volume of radius 10pc) of about 0.1$M_{\odot}$ pc$^{-3}$, which is 5GeV/c$^2$cm$^{-3}$). Galaxies are embedded in dark matter halos with these properties; smaller systems containing dark matter are not observed., Comment: proceedings of The 1st Subaru International Conference, "Panoramic Views of Galaxy Formation and Evolution" (Hayama, Japan, 11-16 Dec 2007) ASP Conf ser

(ABRIDGED) We present a new Schwarzschild orbit-superposition code designed to model discrete datasets composed of velocities of individual kinematic tracers in a dynamical system. This constitutes an extension of previous implementations that can only address continuous data in the form of (the moments of) velocity distributions, thus avoiding potentially important losses of information due to data binning. Furthermore, the code can handle any combination of available velocity components, i.e., only line-of-sight velocities, only proper motions, or a combination of both. It can also handle a combination of discrete and continuous data. The code finds the distribution function (DF, a function of the three integrals of motion E, Lz, and I3) that best reproduces the available kinematic and photometric observations in a given axisymmetric gravitational potential. The fully numerical approach ensures considerable freedom on the form of the DF f(E,Lz,I3). This allows a very general modeling of the orbital structure, thus avoiding restrictive assumptions about the degree of (an)isotropy of the orbits. We describe the implementation of the discrete code and present a series of tests of its performance based on the modeling of simulated datasets generated from a known DF. We find that the discrete Schwarzschild code recovers the original orbital structure, M/L ratios, and inclination of the input datasets to satisfactory accuracy, as quantified by various statistics. The code will be valuable, e.g., for modeling stellar motions in Galactic globular clusters, and those of individual stars, planetary nebulae, or globular clusters in nearby galaxies. This can shed new light on the total mass distributions of these systems, with central black holes and dark matter halos being of particular interest., Comment: ApJ, in press; 51 pages, 11 figures; manuscript revised following comments by referee

We present a synthesis of recent photometric and kinematic data for several of the most dark-matter dominated galaxies. There is a bimodal distribution in half-light radii, with stable star clusters always being smaller than $\sim30$pc, while stable galaxies are always larger than $\sim120$pc. We extend the previously known observational relationships and interpret them in terms of a more fundamental pair of intrinsic properties of dark matter itself: dark matter forms cored mass distributions, with a core scale length of greater than about 100pc, and always has a maximum central massdensity with a narrow range. The dark matter in dSph galaxies appears to be clustered such that there is a mean volume mass density within the stellar distribution which has the very low value of about 0.1$\Msun$ pc$^{-3}$ (about 5GeV/c$^2$ cm$^{-3}$). All dSphs have velocity dispersions equivalent to circular velocities at the edge of their light distributions of $\sim 15$km s$^{-1}$. In two dSphs there is evidence that the density profile is shallow (cored) in the inner regions, and so far none of the dSphs display kinematics which require the presence of an inner cusp. The maximum central dark matter density derived is model dependent, but is likely to have a mean value (averaged over a volume of radius 10pc) of $\sim0.1\Msun$ pc$^{-3}$ (about 5GeV/c$^2$ cm$^{-3}$) for our proposed cored dark mass distributions (where it is similar to the mean value), or $\sim60\Msun$ pc$^{-3}$ (about 2TeV/c$^2$ cm$^{-3}$) if the dark matter density distribution is cusped. Galaxies are embedded in dark matter halos with these properties; smaller systems containing dark matter are not observed., Comment: Astrophys J in press, 35pp, 3 tables, 5 figures

We surveyed 1.75 square degrees of sky near Neptune to an R-band 50% detection efficiency of 25.8 mags (corresponding to radii of about 17 km for an assumed albedo of 0.04). We discovered one new outer satellite, Psamathe (S/2003 N1), about 20 km in radius and having a distant retrograde orbit and moderate eccentricity. Until 2003 Neptune was only known to have two satellites which exhibited orbital signatures indicative of capture. Both of these, Triton and Nereid, are unusual when compared to the irregular satellites of other giant planets. With recent discoveries of four additional satellites by Holman et al. (2004) it is now apparent that Neptune has a distant ``normal'' irregular satellite system in which the satellites have radii and orbital properties similar to those of the satellites of other giant planets. We find that the satellite size distribution at Neptune is not well determined given the few objects known to date, being especially sensitive to the inclusion of Triton and Nereid in the sample. Finally, we note that Psamathe and S/2002 N4 have similar semi-major axes, inclinations and eccentricities. They may be fragments of a once larger satellite., Comment: 10 pages, 9 figures (6 color), 3 tables; Accepted for the July 2006 Astronomical Journal

The history of the Milky Way Galaxy is written in the properties of its stellar populations. Here we analyse stars observed as part of surveys of local dwarf spheroidal galaxies, but which from their kinematics are highly probable to be non-members. The selection function -- designed to target metal-poor giants in the dwarf galaxies, at distances of ~100kpc -- includes F-M dwarfs in the Milky Way, at distances of up to several kpc. Thestars whose motions are analysed here lie in the cardinal directions of Galactic longitude l ~ 270 and l ~ 90, where the radial velocity is sensitive to the orbital rotational velocity. We demonstrate that the faint F/G stars contain a significant population with V_phi ~ 100km/s, similar to that found by a targeted, but limited in areal coverage, survey of thick-disk/halo stars by Gilmore, Wyse & Norris (2002). This value of mean orbital rotation does not match either the canonical thick disk or the stellar halo. We argue that this population, detected at both l ~ 270 and l ~ 90, has the expected properties of `satellite debris' in the thick-disk/halo interface, which we interpret as remnants of the merger that heated a pre-existing thin disk to form the thick disk., Comment: Accepted, Astrophysical Journal Letters

There are two main tidal effects which can act on the Local Group dwarf spheroidals (dSphs): tidal stripping and tidal shocking. Using N-body simulations, we show that tidal stripping always leads to flat or rising projected velocity dispersions beyond a critical radius; it is $\sim 5$ times more likely, when averaging over all possible projection angles, that the cylindrically averaged projected dispersion will rise, rather than be flat. In contrast, the Local Group dSphs, as a class, show flat or falling projected velocity dispersions interior to $\sim 1$ kpc. This argues for tidal stripping being unimportant interior to $\sim 1$ kpc for most of the Local Group dSphs observed so far. These results have important implications for the formation of the dSphs and for cosmology. A tidal origin for the formation of these Local Group dSphs (in which they contain no dark matter) is strongly disfavoured. In the cosmological context, a naive solution to the missing satellites problem is to allow only the most massive substructure dark matter halos around the Milky Way to form stars. It is possible for dSphs to reside within these halos ($\sim 10^{10}$M$_\odot$) and have their velocity dispersions lowered through the action of tidal shocks, but only if they have a central density core in their dark matter, rather than a cusp. dSphs can reside within cuspy dark matter halos if their halos are less massive ($\sim 10^{9}$M$_\odot$) and therefore have smaller central velocity dispersions initially (abridged)., Comment: 15 pages, 9 figures. Accepted for publication in MNRAS. No change from previous version

We present an improved analytic calculation for the tidal radius of satellites and test our results against N-body simulations. The tidal radius in general depends upon four factors: the potential of the host galaxy, the potential of the satellite, the orbit of the satellite and {\it the orbit of the star within the satellite}. We demonstrate that this last point is critical and suggest using {\it three tidal radii} to cover the range of orbits of stars within the satellite. In this way we show explicitly that prograde star orbits will be more easily stripped than radial orbits; while radial orbits are more easily stripped than retrograde ones. This result has previously been established by several authors numerically, but can now be understood analytically. For point mass, power-law (which includes the isothermal sphere), and a restricted class of split power law potentials our solution is fully analytic. For more general potentials, we provide an equation which may be rapidly solved numerically. Over short times ($\simlt 1-2$ Gyrs $\sim 1$ satellite orbit), we find excellent agreement between our analytic and numerical models. Over longer times, star orbits within the satellite are transformed by the tidal field of the host galaxy. In a Hubble time, this causes a convergence of the three limiting tidal radii towards the prograde stripping radius. Beyond the prograde stripping radius, the velocity dispersion will be tangentially anisotropic., Comment: 10 pages, 5 figures. Final version accepted for publication in MNRAS. Some new fully analytic tidal radii have been added for power law density profiles (including the isothermal sphere) and some split power laws

We present the results from a suite of N-body simulations of the tidal stripping of two-component dwarf galaxies comprising some stars and dark matter. We show that recent kinematic data from the local group dwarf spheroidal (dSph) galaxies suggests that dSph galaxies must be sufficiently massive ($10^9 - 10^{10}$M$_\odot$) that tidal stripping is of little importance for the stars. We discuss the implications of these massive dSph galaxies for cosmology and galaxy formation., Comment: 4 pages, 1 figure, to appear in the proceedings of the IAUC198 "Near-Field Cosmology with Dwarf Elliptical Galaxies", H. Jerjen & B. Binggeli (eds.). Comments welcome

We present a deep optical survey of Uranus' Hill sphere for small satellites. The Subaru 8-m telescope was used to survey about 3.5 square degrees with a 50% detection efficiency at limiting red magnitude m = 26.1 mag. This magnitude corresponds to objects that are about 7 km in radius (assuming an albedo of 0.04). We detected (without prior knowledge of their positions) all previously known outer satellites and discovered two new irregular satellites (S/2001 U2 and S/2003 U3). The two inner satellites Titania and Oberon were also detected. One of the newly discovered satellites (S/2003 U3) is the first known irregular prograde of the planet. The population, size distribution and orbital parameters of Uranus' irregular satellites are remarkably similar to the irregular satellites of gas giant Jupiter. Both have shallow size distributions (power law indices q~2 for radii > 7 km) with no correlation between the sizes of the satellites and their orbital parameters. However, unlike those of Jupiter, Uranus' irregular satellites do not appear to occupy tight distinct dynamical groups in semi-major axis versus inclination phase space. Two groupings in semi-major axis versus eccentricity phase space appear to be statistically significant., Comment: 10 Figures (2 color), 3 Tables, to appear in AJ January 2005

We surveyed the Hill sphere of Mars for irregular satellites. Our search covered nearly the entire Hill Sphere, but scattered light from Mars excluded the inner few arcminutes where the satellites Phobos and Deimos reside. No new satellites were found to an apparent limiting red magnitude of 23.5, which corresponds to radii of about 0.09 km using an albedo of 0.07., Comment: 5 figures (1 color), 2 Tables, to appear in AJ Nov. 2004

The nearby Ursa Minor dwarf spheroidal (UMi dSph) is one of the most dark matter dominated galaxies known, with a central mass to light ratio roughly equal to 70. Somewhat anomalously, it appears to contain morphological substructure in the form of a second peak in the stellar number density. It is often argued that this substructure must be transient because it could not survive for the > 10 Gyr age of the system, given the crossing time implied by UMi's 8.8 km/s internal velocity dispersion. In this paper, however, we present evidence that the substructure has a cold kinematical signature, and argue that UMi's clumpiness could indeed be a primordial artefact. Using numerical simulations, we demonstrate that substructure is incompatible with the cusped dark matter haloes predicted by the prevailing Cold Dark Matter (CDM) paradigm, but is consistent with an unbound stellar cluster sloshing back and forth within the nearly harmonic potential of a cored dark matter halo. Thus CDM appears to disagree with observation at the least massive, most dark matter dominated end of the galaxy mass spectrum., Comment: Astrophysical Journal (Letters), in press

We present the first clear evidence for an extended dark matter halo in the Draco dwarf spheroidal galaxy based on a sample of new radial velocities for 159 giant stars out to large projected radii. Using a two parameter family of halo models spanning a range of density profiles and velocity anisotropies, we are able to rule out (at about the 2.5 sigma confidence level) haloes in which mass follows light. The data strongly favor models in which the dark matter is significantly more extended than the visible dwarf. However, haloes with harmonic cores larger than the light distribution are also excluded. When combined with existing measurements of the proper motion of Draco, our data strongly suggest that Draco has not been tidally truncated within ~1 kpc. We also show that the rising velocity dispersion at large radii represents a serious problem for modified gravity (MOND)., Comment: to be published in ApJL; 5 pages, 4 figures

We report submillimeter 450 and 850 microns dust continuum observations for comet C/2012 S1 (ISON) obtained at heliocentric distances 0.31-0.08 au prior to perihelion on 2013 November 28 (rh=0.0125 au). These observations reveal a rapidly varying dust environment in which the dust emission was initially point-like. As ISON approached perihelion, the continuum emission became an elongated dust column spread out over as much as 60″ (greater than 10(exp 5) km in the anti-solar direction. Deconvolution of the November 28.04 850 microns image reveals numerous distinct clumps consistent with the catastrophic disruption of comet ISON, producing approximately 5.2×10(exp 10) kg of submillimeter-sized dust. Orbital computations suggest that the SCUBA-2 emission peak coincides with the cometʼs residual nucleus.